'Theoretical and experimental studies in pulsed nuclear 

QUADRUPOLE RESONANCE SPECTROSCOPY : APPLICATION OF 
TENSOR OPERATOR FORMALISM 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


1 


by 

RAVINDER REDDY 


m the 

DEPARTMENT OF CHEMISTRY 
INDIAN INSTITUTE OF TECHNOLOGY. KANPUR 

MAY. ISSt 



- Rth-T 


f'- 8 NOV 1909 

'"’^^ARY 

lovk4l^. 


T/, 

5^3' 0B77 



II>K:D I CATE 3D 
TO 

M-V" F> A “F? W KI *T* ^ 



i 


DEF/vRTHEKT OF CHEMISTPY 
TUTE OF TECHKOLCGY KAi:: 


I A 


CERTIFICATE OF COURSE WORK 


This is to certify that Mr. Ravinder Reddy has satis- 
factorily completed all the courses required for the Fh.D. 
degree progranmo. These courses include: 


Chm 505 Principles of Organic Chemistry 
Chm 521 Chemical Binding 

Chm 524 Modern Physical Methods in Chemistry 
Chm' 525 Principles of Physical Chemistry 
Chm 5 34 Electronic?;*, for Chemists 
Chm 545 Principles of Inorganic Chemistry 
Chm 634 Symmetry and Molecular Structure 


Chm 800 General Seminar 
Chm 801 Special Seminar 
Chm 900 Ph.D. Thesis 


Mr, Ravinder Roddy v;as admitted to the >'<indidacy of tho 
fh.D. degree in August 1983 after he succes.vfully completed the 
•rfrltton and oral qualifying examinations. 


M- 

(N. Sathyamurthy) ^ 


Professor and Acting Head, 
Department of Chemistry, 


I.I.T., KANPUR 



(S. Sarkar) 
Convener, 

Departmental Post- 
Graduate Committee, 
Dopartm'nd’ oi- s try , 

KANPUR 



ii 


CERTIFICATE II 

Certified that the work contained in this thesis entitled: 
"Theoretical and E;:perimental Studies in Pulsed Kuclear Quadru- 
pole Resonance Spectroscopy: Application of Tensor Operator 
Formalism" has been carried out by Mr. Ravinder Reddy under my 
supervision and the same has not been submitted elsewhere for 
a degree. 



(P. T. Narasimhan) 
Professor 

Departraent of Chemistry, 
I.I.T. Kanpur 

Thesis Supervisor 


Kanpur, 
May 1988. 



iii 


STATEMENT 

I hereby declare that the matter embodied in this thesis 
is the result of investigations carried out by me in the Depart- 
ment of Chemistry/ Indian Institute of Technology, Kanpur, India, 
under the supervision of Professor P.T. Narasimhan. 

In keeping with the general practice of reporting scienti- 
fic observations, due acknowledgements have been made wherever 
the work described is based on the findings of other investiga- 
tors . 


Ravinder Reddy 

Kanpur 
May 1988. 



iv 


ACKNOWLEDGEMENTS 


It is a great pleasure to express my gratitude to Professor 
P.T. Narasirnhan who introduced me to the area of magnetic reson- 
ance in general and NQR in particular. The past several years 
of my association with him have helped shape my attitudes towards 
academic as well as . non- academic matters. His inspiring guidance 
and constant encouragement have been the main factors responsible 
for the completion of this work. 

I also thank. 

Dr. (Mrs.) Amrita Sheikh, Messrs Arun L. Bhavsar, A.K. Dubey and 
A. Ramamoorthy for providing an ideal working atmosphere in the 
laboratory. 

Dr. Narsimha Reddy for his invaluable help during the course of 
this work. 

Professor P. Raghunathan for his encouragement and moral support. 

Professor S. Chandrasekaran and K. Sathyamurthy for their help, 
in particular at the initial stages of my programme. 

Dr. N. Chandrakumar for useful discussions on MQNQR. 

All my friends who helped me in many ways and made my stay at 
IIT-K a pleasant one, 

Mr. R.D. Singh for efficient typing of the manuscript. 

Messrs Umesh and Ram Singh Chauhan for their help in the 
laboratory. 

Ravinder Reddy 



SYNOPSIS 


The thesis deals with theoretical and experimental studies 
in pulsed nuclear quadrupole resonance (NQR) spectroscopy. Tensor 
operator formalism has been employed for the first time to ana- 
lyse pulsed NQR experiments in a variety of situations. 

The thesis consists of five chapters. Chapter I deals with 
the basic principles, methods of detection and applications of 
NQR spectroscopy. A description of the pulsed NQR spectrometer 
employed in the present investigations is also presented in this 
chapter. A brief outline of the various theoretical approaches 
available for the description of pulsed NQR is given. Scope of 
the present work is presented at the end of this chapter. 

Chapter II is concerned with the application of tensor 
operator formalism for pure NQR. Since quadrupolar nuclei may 
be present in axially symmetric or non-axially symmetric elec- 
tric field gradients (efg) the study of the responses of quadru- 
polar nuclei to radio frequency (r.f.) perturbations in these 
situations is of considerable importcince. Of special interest 
is the study of response of spin 3/2 nuclei to r.f. pulse excita- 
tions in the presence of a weak Zeeman field in order to obtain 
both the quadrupole coupling constant (e qQ) and the asymmetry 
parameter n . 

Using the tensor operator formalism the response of non- 
interacting nuclear spin systems with 1=1, 3/2, 5/2 in single 
crystals to single- and two- r.f. pulse excitations have been 
evaluated. The following cases have been considered: 



i) Spin I = 3/2, 5/2 nuclei in .axially symmetric efg's. 

ii) Spin 1=1, 3/2 nuclei in non-axially symmetric efg's. 

iii) Spin I = 3/2 in axially symmetric field gradients in thfe 
presence of a weak Zeeman field applied at an angle 'O' 
w.r.t. the principal Z-axis of the efg. 

iv) Spin I = 3/2 in non-axial field gradients, in the presence 
of Zeeman field applied along principal X-, Y-, or Z-axis 
of the efg. 

In all these calculations relaxation effects are not considered- 
The time evolution of complete set of tensor operators for 
I = 3/2 under the action of r.f. Hamiltonian (including the phase 
of the r.f.) has been investigated in detail. For the purpose of 
illustration experimental recordings of the pulse responses in 
some cases are also presented. 

Wherever possible the present results have been compared 
with available literature reports and agreement has been found. 
The present set of results are, however, more exhaustive than 
those given in previous literature reports. 

Chapter III deals with the application of tensor operator 
formalism to the study of Zeeman-perturbed spin echo envelope 
modulations (ZSEEM) in spin I = 3/2 nuclei. In the first section 
of this chapter a brief discussion of the two-pulse ZSEEM experi- 
ment is given and illustrated using experimental recordings of 

35 

ZSEEM -patterns obtained from cl in polycrystalline specimen of 
AgClOj • In the second section tensor operator formalism has been 
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used for the calculation of the expression for the ZSEEM pattern. 
The present results agree with those of Rarnacliandran and Nara- 
simhan who employed the matrix transformations approach. 

Chapter IV is concerned with the experimental and theore- 
tical investigations of responses of quadrupolar nuclei to multi- 

3 5 

pie-pulse sequences. The nuclei investigated are: Cl (I = 3/2) 

in (i) single crystal of NaClO^/ Cii) polycrystalline samples of 
NaClO^. KCIO 3 , SbCl 3 and ^®^Re (both with I = 5/2) in 

polycrystalline sample of KReO^. The pulse sequences studied 
are: (a) spin-locking sequence (SLS) : (t-© 

(b) phase-alternated pulse sequence (PAPS) : (rc/2) -(T-(e) 

DC DC 

2t:-(9) -t) and (c) WAHUHA multiple pulse sequence: (^/2) - 

DC jn DC 

T(T:-(n/2) ,-t-(tt/2)„-2t- (Tr/2) -x-(n/ 2 ) -x-) Here, (0) , 

—X y — y Xu 0 

represents a jjulse of flip; angle Q and phase / with t correspond- 
ing to the pulse separation. 


In Section IV. A a brief review of multiple-pulse studies 

in NQR is given. Section IV. B presents the experimental results 

with the aforesaid sequences. Under the influence of SLS and 

PAPS the establishment of a quasi-equilibriurn state of the 

35 

magnetization has been observed both in Cl (in single crystal 
and powder samples) and Re, Re in KReO^. In all the cases 
studied the decay constant, ^2^ quasi-equilibrium magne- 

tization (QEM) has been found to be much larger than the spin- 


spin relaxatio: 


time T 2 . The dependence of the decay of the 


QEM on the oarlineters such as ©, and resonance off-set. Aw 
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has been investigated. From these studies it is inferred that 
the spin-locking phenomenon can be observed in single crystals 
as well as in powder samples with arbitrary ©, JZT values of the 
sequence pulses contrary to the theoretical predictions in the 
literature . 

23 35 

Section IV. C presents the results of a Na- Cl double 
resonance experiment carried out under the multiple pulse spin- 
locking condition in a polycrystalline sample of NaClO^- Section 
IV. D outlines the theoretical analysis of PAPS/ WAHUHA and Mans- 
field-Rim-Elleman-Vaughan (MREV-8) sequences in the case of 
spin I = 3/2 nuclei. Using tensor oioerator formalism and 
average Hamiltonian theory it has been shown that in the case 
of single crystals under the action of PAPS the internal inter- 
actions due to (i) efg inhomogeneities ( ' (ii) torsional 

oscillations ( snd (iii) heteronuclear dipole-dipole inter- 
actions ( , ) which contribute to the NQR line width are aver- 

het 

aged out. In the case of polycrystalline samples it has been 
possible to show that under PAPS '^tors averaged 

out. For both the WAHUHA and MREV-8 sequences it has been possi- 
ble to show that in the case of single crystals ^j_nhom "^tors 

are not averaged out. 

Results of theoretical investigation of double quantum co- 
herence in pure NQR of non-interacting nuclear sjjins with I = 5/2 
for the case of axially symmetric .efg in a single crystals are 
presented in Chapter V. Employing the tensor operator formalism 
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it has been shown that the creation and detection of double quan- 
tum coherence can be done by the application of a pu2se scheme of 
the form (n:/2) ( 2ou) (t) (co,n) - t- (t) . The effect of resonance 

off-set on the double quantum coherence has also been considered. 

The thesis concludes with some critical comments on the 
merits and demerits of the tensor operator formalism as applied 
to pulsed NQR spectroscopy. 
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In this chapter a brief "discussion on the principles, 
applications and detection of Nuclear Quadrupole Resonance 
(NQR) and some theoretical aspects of pulsed KQR are presented. 

Section I. A outlines the basic principles of NQR while 
some applications of NQR are presented in Section I.B. A 
brief outline of the raethods of detection of NQR and descrip- 
tion of the pulsed NQR spectrometer employed in the present 
work is given in Section I.C. Section I.B pi-esents the basic 
theoretical outlines concerning pulse responses in NQR. In 
Section I.E the scope of the present work is described. 


I. A INTRODUCTION TO NQR 

Nuclear quadrupole resonance spectroscopy deals with 
interactions between radio frequency (r.f.) field and nuclear 
spins with I >1 in crystalline (solids) samples. I is the 
nuclear spin angular momentum measured in units of "h. A nucl- 
eus with spin I > 1 exhibits a non-spherical charge distribu- 
tion and possess a magnetic dipole moment as well as an elec- 
tric quadrupole moment and hence it can interact with surround- 
ing magnetic fields and also with electric field gradients 
(efg) . Such nuclei probe magnetic, electrostatic and electro- 
dynamic properties of chemical systems and thus yield valuable 


information 
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The scalar quadrupole moment (eQ) is a measure of the 
deviation of the nuclear charge distribution from spherical 

lAr 

symmetry [l#2] . when a nucleus with a quadrupole moment is 
placed in an inhomogeneous electric field arising from, exter- 
nal charges, it orients itself in different directions with 
different energies.' The observed NQR spectra correspond to 
the dipole transitions/ between such energy levels induced by 
an oscillating r.f. field. Like in NMR these transitions are 
magnetic dipole transitions. 

The theoretical treatment of the interaction of a qua- 
drupolar nucleus v/ith the surrounding electronic charge cloud 
in atoms v;as originally given by Casimir [3j. From a study 
of atomic spectra Schiilar and Schmidt [4j have proved the 
existence of nuclear quadrupole moments, huclear quadrupole 
interaction in m.olecules was first observed by Kellogg et al . 
[ 5 ] in 1939 using the molecular beam technique. The first 
experimental observation of pure NQR was made by Dehmelt and 
Kruger [6J. 

In the principal axes system of efg, the Hamiltonian 
corresponding to the interaction between the electric quadru- 
pole moment and electric field gradient at the nucleus can be 
written as [?] 


^References appear at the end of each chapter. 
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% = 4iT5??iy - r 


(i.i) 


where rj the asymmetry parameter is defined as 


V 


^yD~^YY 


V 


zz 


( 1 . 2 ) 


r) is a measure of the deviation of the efg from cylindrical 

2 

symmetry. Conventionally, e qQ is referred to as the nuclear 
quadrupole coupling constant (QCC) . The principal axes sys- 
tem of efg tensor has been chosen such that the diagonal compo- 
nents of the efg, '^YY '^ZZ condition 



> 


V. 


YY 



If we assume that the electric field at the nucleus is produ- 
ced entirely by charges completely external to the nucleus 
then the Laplace equation +- V^y + ^zz ~ ^ satisfied 


and we have only two independent parameters V 


ZZ 


eq and t). 


The two conditions, viz., V 


zz 


> ^yYl> 


XX 


and V^, + Vyy 


+ V„ = 0 lead to the result that the value of 7) can vary 


between 0 and 1. h is zero for the axially asymmetric case. 


In order to arrive at the energy levels of the quadru- 
polar nucleus, one needs to evaluate the matrix elements of 
the type = <^i,m| '^Q|l,m'^ . For the case of axially 

symmetric field gradient ( 7)= 0) the Hamiltonian matrix is 
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diagonal in the representation and the energy levels are 

Zj 

given by 


E 


m 


[3 - I(I+l)j 

41(21-1) 


.. (1.3) 


where m is the Z-component of the nuclear spin I. Since the 
energies of the states with .are the same the energy levels 
are degenerate in m. 

The presence of asymmetry par ante ter 77 / mixes the states 
having Am = + 2 and hence the Hamiltonian matrix is not dia- 
gonal in the 'm' representation. In this case one has to 
construct linear combinations of spin functions |l/m)>' , [l,m'^ 
having Am = + 2 , and calculate the coefficients in these 
linear combinations from the energies obtained as the eigen- 
values of the matrix. For various values of nuclear 

spin I, the secular equations have been tabulated in the book 
by Das and Hahn [s]. 

For the case of spin 1=1 the energies of states with 
m = 0, +1 are given respectively by 


E = 
o 


e^q-Q 


.. (1.4) 



(1 + 77 ) 


e!aQ 

4 


(1.5) 


Hence the frequencies of transition are given by 
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V 


4* 


3e^qQ 

“4 


(1 + 77/3) 


.. ( 1 . 6 ) 


except for nuclear spin I = 3/2 in all other cases the nuclear 

spin energy levels will be established in such a way that more 

than one transition of Am = +1 between these levels is possible. 

Thus by observing these transitions the quadrupole coupling 
2 

constant e qQ and 77 can be determined. 


In the case of spin I = 3/2 we still get two doubly- 
degenerate levels which are given by [8,9] 



(1 + 77 ^ 3 )^^^ 
^ (1 t 7)2/3) 


.. (1.7) 

.. ( 1 . 8 ) 


where A^, B^, A^, are appropriate linear combinations of pure 

m functions. As there are only two energy levels there is only 

2 

one transition and it is not possible to obtain both e qQ and 77 
by studying the pure NQR of spin I = 3/2 nuclei. In order to 
obtain both these parameters one has to do a Zeeman NQR study. 
Zeeman NQR studies have been reported both on single crystals [ 9 ] 
as well as on powder sample [lo] to obtain these parameters. 

Until ^recently the eigenvalues of matrix in the case 
of spin I >3/2 were not available in the analytical form. 

Bersohn [H] dealt with the secular equations for I,> 3/2 by 
a perturbation procedure which is convergent for 7)^o.25. , 
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Cohen [l2] has solved the secular equations for I = 5/2/ 7/2 
and 9/2 numerically and tabulated eigenvalues for values of 77 
ranging from 0.1 to 1.0. Recently/ there have been reports 
where the exact analytical solutions of the secular equations 
for spin I = 5/2/ 7/2 and 9/2 are obtained [l3]. Analytical 
expressions for eigenvalues in the case of spin I > 3/2 in the 
presence of a weak Zeeman field are also available in the lite- 
rature [l4,15]. Using a spherical tensor technique Krishnan 
and Sanctuary [I 6 ] have obtained secular equations in frequen- 
cies of spin 1 = 1 / 3/2 and 5/2 with non- axial field gradients, 
and solved exactly in terms of quadrupole coupling constant and 
asymmetry parameter. 

Once the eigenvalues are known the eigenfunctions may be 
obtained by the usual procedure.. The transitions between the 
doubly degenerate energy levels and induced by 

the application of a linearly polarized r.f. field which inter- 
acts with magnetic dipole moment of the nucleus. The frequency 
of the transition is given by 



(m m' ) 



.. (1.9) 


For the case of axial symmetry/ transitions with Am = +1 are 
only allowed. But in the case of finite r), transitions corres- 
ponding to Am ^+1 are also allowed but they are generally weak. 
In NQR only allowed transitions are studied in order to obtain 
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the field gradients at the nucleus. 

The relative intensities of NQR spettral lines aris- 

ing from transitions between initial and final Quadrupole 
states and , are x^roportional to |<(i| 

is the perturbing Hamiltonian due to the action of the 
r.f. field on the magnetic dipole moment of the nucleus. The 
intensities of the pure quadrupole lines for various values of 
I are discussed in the literature [9,12]. 

Like in other branches of magnetic resonance spectro- 
scopy, in NQR also the relaxation behaviour of the nuclear spin 
ensemble is characterised by the spin-lattice relaxation time 
T^ and the spin-spin relaxation time T^. These relaxation pro- 
cesses play an important role in determining NQR line widths. 
Experimentally observed NQR lines possess finite widths which 

can be related to the apparent spin-spin relaxation time T 2 * . 

-2 -5 

The relative line widths, vary between 10 and 10 of the 
transition frequency. We now consider, briefly, the factors 
that contribute to the line widths in NQR. 

In general, width, shape, and structure of NQR signals 
can have contributions from both 1) static and 2) dynamic 
effects [s]. Among static effects one has to consider both 
i) magnetic and ii) electrical effects. 

i) Magnetic Effects: 

Two types of magnetic interactions between the resonant 
nucleus and the other nuclei of the same or different soecies 
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are possible. These are, a) direct magnetic dipole-dipole 
interactions [l7j. b) Indirect spin-spin interactions 
(J-interaction) [l7]. 

(a) Direct dipole-dipole interaction between the reso- 
nant nucleus and close neighbouring nuclei leads to the fine 
splitting of quadrupolar spectra, whereas the magnetic dipole- 
dipole interaction between the resonant nucleus and other 
distant neighbouring nuclei, other than the very close neighr- 
bour, broadens the fine structure produced by the inter- 
action of the close neighbour. 

(b) The indirect spin-spin interaction between two 
nuclei, originates from the magnetic interactions of two 
nuclei v/ith the electrons that form the chemical bond between 
them [ 17 ]. This interaction can lead to fine splitting of the 
quadrupole spectra if its short range effect is sufficiently 
strong to overcome the dipole-dipole interaction. Broaden- 
ing [is] due to J-interaction can take place only if it occurs 
with equivalent strength between a resonant nucleus and a 
number of neighbours. 

ii) Electrical Effects: 

The main sources of electrical effects are dislocations 
and strains in the crystal or powder grains. These are found 
to broaden the NQR spectra [l9-2l]. Qualitatively, the broaden- 
ing due to strains and dislocations can be understood as arising 



from random distortions of intermolecular interactions, and a 
corresponding random distribution of field gradients at nuclear 
sites. Presence of impurities also causes a random distribu- 
tion in environments around different nuclei and hence leads to 
a broadening of quadrupolar lines. 


2 ) Dynamic Effects: 

Among the dynamic effects the motions arising due to 
torsional oscillations of molecules are important. In NQR, 
torsional motions contribute most to the spin lattice relaxa- 
tion. Since the spin lattice relaxation limits the life time 
of the excited state, it contributes to a line broadening 
known as the "life- time broadening". Hence the contribution 
to line widths from dynamic effects are related to T^. 

Thus, in general, a composite line shape g(co) of the NQR 
line can be expressed as 

00 oo 

=_/_J 9ayn (U.U') (u'' du» 

.. ( 1 . 10 ) 


where ^dip ^ele line shape functions corres- 

ponding to the dynamic effects, static magnetic dipolar inter- 
actions and static electric effects, respectively and Wq is the 
the NQR frequency. Now an approximation can be made that the 
net half width of the .NQR signal gCw) is given as a sum of the 
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half-widths of the individual line shape functions ^dip 

and Qgjg* The net half-v/idth 6ci) can be represented by 

6co = ^ .. (1. 11) 

dyn dip ele ' 

normally the net line shape in NQR can be described by a 
Gaussian. Hence T 2 * can be defined as 

^ 2 ( 7T log 2 ) 

T = 2 _ (1.12) 

^ 6co 

So far we have discussed the basic theoretical aspects 
of NQR. In the next section we present some of the applica- 
tions of NQR spectroscopy. 


1*B SOME APPLICATIONS OF NQR 

NQR spectroscopic investigations yield information about 

2 

the quadrupole coupling constant, e qQ, and asymmetry para- 
meter, T) , Of the efg. These parameters can be used to get 
considerable information regarding the electronic environment 
around a quadrupolar nucleus which helps in understanding the 
chemical bond and other related properties of the molecules or 
ions in which the nucleus is present [8,22-24]. However, pure 
NQR study is possible only in the solid state. This is because 
of the fact that the efg tensor like the dipolar 'interaction 
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tensor is traceless and the rapid motion present in liquids 
and gases averages it out to zero. This does not imply that 
quadrupole interactions do not exist in liquids and gases. 

The information about quadrupole interactions in liquids and 
gases can be obtained by other spectroscopic methods like NMR, 
EPR/ microwave and molecular beam techniques. Thus, the 
information from NQR complements the' information from other 
techniques. The chemical applications of NQR have been review- 
ed extensively in the literature [8,23-30J. In what follows 
we present briefly typical applications of NQR spectroscopy. 

It must be emphasized here that no attempt for a thorough 
review of the literature is made. 

2 

The quadrupole coupling constant, e qQ, is a product of 
nuclear quadrupole moment (eQ), a nuclear property, and (eq) , 
the principal Z-component of the efg tensor, which is a mole- 
cular property. Hence, if we know one of these, the other 
quantity can be determined using the experimental value of QCC. 
The quantity 'q' could be estimated if the charge distribution 
over the molecule is known. This requires a knowledge of the 
wave functions for the electrons in the molecules. 

Both ^ initio and serai- empirical wave functions have 
been used in the theoretical calculation of the efg's. Ab 
initio level calculations are expected to yield accurate efg's. 
Results of ab Initio calculations of efg's for free molecules 
have been reviewed extensively [31-34]. Recently, Palmer [35j 
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and Brown et al . [36] have reported initio calculations of 
14 

N quadrupole coupling constants using different kinds of 
basis sets. However, in the case of efg calculations on large 
molecules, as the initio calculations become cumbersome and 
require large computer time, one has to resort to the approxi- 
mation methods like . empirical and serni-erapiirical methods. The 
two approximation methods, the Townes-Dailey method [37] and 
the Cotton-Harris method [38] are very commonly used for the 
calculation of efg's. Molecular orbital approach has been very 

popular in this regard. Bray et al. [39] have analyzed the 
14 

results of N NQR data in solids using Townes-Dailey approach. 

In the Townes-Dailey approach, Sternheimer effects [40-42] are 
not considered. These effects play an important role in the 
evaluation of efg's when deformation of electron core takes 
place either by efg due to outer orbitals or by nuclear quadru- 
pole moment. Results of quadrupole coupling constants obtained 
using several semi- empirical methods have been reviewed [43]. 

In covalent molecules a knowledge of efg's can yield 
information on the ionic character of the chemical bond invol- 
ving the quadrupolar nucleus and the hybridization of the atomic 
orbitals involved in the bond [8,22,23]. Asymmetry paramter 77 
does not involve any properties of the nucleus and it can be 
interpreted directly in terms of the electron distribution in 
molecules. NQR spectral parameters can also be correlated with 
other chemical- and physicochemical constants such as Hammet- 

Toft 0 parameters, electronegativities and dissociation constants 
[8,22,23] . 
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NQR can be used to distinguish between "Chemically inequi- 
valent" and "physically inequivalent" sites of the same quadru- 
polar nucleus in a crystal [s], A lattice site is said to be 
chemically inequivalent to another lattice site if the com- 
ponents of principal efg tensor at the respective sites are 
different. The physically inequivalent lattice sites have at 
the respective sites the same magnitude for the principal efg 
components but they are oriented in different directions. Hence, 
chemically inequivalent sites have different NQR frequencies, 
whereas, physically inequivalent sites generally, have the same 
frequency. 

As the efg is insensitive to isotopic substitution, the 
ratio of quadrupole moments of different isotopes can be obtain- 
ed' from the ratio of quadrupole coupling constants of the same 
compound with different isotopes [44-46J. 

The efg and hence the NQR frequencies are quite sensitive 
to temperature changes. Systematic studies of the temperature 
dependence of the NQR frequencies yield useful information 
regarding molecular motions in solids [8,28]. As precise fre- 
quency measurements are possible, temperature dependence of NQR 
frequency has been used for accurate measurement of tempera- 
ture [47,48]. Normally the temperature dependence of the NQR 
frequency do^^/dT is negative. Anomolous temperature depend- 
ence of NQR frequencies in coordinate compounds have been 
reviewed in the literature [49]. 
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NQR studies of pressure dependence can reveal intermole- 
cular interactions present in the crystal [50]. Several studies 
on the pressure dependence of NQR frequencies have been made in 
the literature [28,50-55], 

Temperature and pressure changes on a crystal can induce 
phase transitions in solids. Since NQR frequencies and inten- 
sities undergo changes at these phase transition points, NQR 
can be used for the study of phase transitions [8,27,28,56-60], 

As mentioned in the previous section relaxation processes 
play an important role in determining line widths. is found 

to be more sensitive to the dynamic structure of the solid than 
[6l]» T 2 can have contributions from both static and 
dynamic effects. In quadrupolar spin systems, T^ can be quanti- 
tatively related to the nature and frequency of the molecular 
and lattice motions in the solid phase [61-63]. Since these 
relaxation times are temperature and pressure dependent, accu- 
rate measurements of NQR parameters as a function of tempera- 
ture and pressure can yield information about* structural and 
dynamical effects in .solid state [64], 

NQR frequency changes due to the transition from the con- 
ducting to superconducting state have been observed [65]. Care- 
ful investigation of NQR parameters is expected to lead to a 
better understanding of the superconducting state of the mate- 
rial. Recently, several workers have investigated the high 
temperature superconducting oxides by NQR spectroscopy [66, 67], 
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Using NQR spectroscopy the existence of the quanfum mecha- 
nically predicted Berry's phase has been demonstrated experi- 
mentally [68]. 

Nuclear spin with I >2 possess/ in addition to the quadru- 
pole moments, electric hexadecapole moments [20]. Usually hexa- 
decapole coupling constant values are small, which lie between 
few hundred Hz to few KHz, hence the direct detection of this 
interaction is rather difficult. Several workers have used 
accurate measurements of NQR frequencies to obtain hexadeca- 
pole coupling constant [20,69-74], 

From the above discussion it is clear that NQR is one of 
the most important tools for the investigation of static and 
dynamic phenomena in solids. The research work in NQR is ref- 
lected in the nine bi-annual "International Symposia" that were 
held in the recent past. Proceedings of the first two sympo- 
sia have been published in the form of books [75,76]. The 
papers presented in the subsequent symposia have been published 
in scientific journals namely, "Journal of Magnetic Resonance" 
"Journal of Molecular Structure" and "Zeitschrift fur Natur- 
forschung" . The papers presented at the recent ninth sympo- 
sium (Jan. 1988) will be published in the "Journal of molecu- 
lar structure". 

in the next section of this chapter we outline, briefly, 
the methods of detection of NQR and also present a brief 
description of the pulsed NQR spectrometer employed in the 
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experimental investigations in the present thesis. 

I.C DETECTION OF NUCLEAR QUADRUPOLE RESONANCE 

The information about nuclear quadrupole interactions 
can be obtained from a variety of the spectroscopic techniques. 
For example/ to name a few/ in the gaseous state, the quadru— 
pole interactions can be measured by molecular microwave spec- 
troscopy [ 77 ] and NMR spectroscopy can be used to measure the 
quadrupole interactions in liquid state [78]. Quadrupole 
interactions in excited molecular states in solids can be 
studied by a novel technique, known as optically detected 
magnetic resonance (o.d.m.r.) [79]. For a detailed discussion 
of various techniques for measuring quadrupole interactions the 
reader is referred to excellent articles by Smith [30, 80 ]. In 
this section we deal only with the methods of pure nuclear qua- 
drupole resonance in solid state. 

Like in NMR, the direct detection of Nqr can be made 
using continuous wave (c.w.) or pulsed and Fourier transform 
(ft) techniques. In what follows, we present a brief account 
of these methods and no attempt has been made for a rigorous 
survey of literature. 

I.C.l Continuous Wave Techniques ; 

In continuous wave spectrometers the sample is subjected 
to an r.f. field, by placing it in a coil which is part of an 
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r.f. oscillator. The same oscillator usually detects the NQR 
signal induced in the coil. The crossed-coil method [Sl] of 
detection which is commonly used in NMR spectroscopy is not 
possible in pure NQR as there is no net component of induced 
signal at right angles to the transmitter coil' because of the 
opposite senses of quadrupole precession associated with the 
+ m^ degeneracies of the quadrupole state. However, if a 
small magnetic field is applied then this degeneracy is liftdd, 
and the crossed-coil method can be used. Since a component of 
nuclear induction then appears normal to the axis of the r.f. 
coil. A few NQR studies have been done v/itiri crossed coil 
arrangements [82-84 j. Since the magnitude of the splittings 
of nuclear energy levels by the interaction of a nuclear qua- 
drupole with electric field gradients in solids can range from 
zero to several thousand MHz, NQR spectrometers are required to 
operate over a large range of frequencies. Because of this 
reason c.w. NQR s^rjectrometers are based on frequency- swept 
oscillators, in which the frequency can be changed by slowly 
varying the tuning capacitor in the tank circuit of the oscilla- 
tor. 

Most of the c.w. NQR spectrometers are based on one or 
other of four main types of oscillators, namely, 

i) the marginal oscillator [85] 

ii) the regenerative oscillator [86] 

iii) the superregenerative oscillator (SRO) [87,88] 
and iv) Limited oscillator [89], 
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Among these^ SRO's have the highest sensitivity and hence they 
have been used, extensively/ for NQR studies. However, these 
oscillators are associated with frequency instabilities and it is 
very difficult to get pure absorption line shapes using simple 
SRO's. Injection- and phase-locking of an SRO to a stable 
external crystal oscillator can improve its frequency stability 
and line shape response [90-93]. In such spectrometers signal 
averaging can be used to improve s ignal-to- noise ratio of weak 
NQR signals [93] . For a detailed discussion of c.w. NQR spec- 
trometers the reader is referred to [90-94]. 

I.C.2 Pulsed Techniques : 

For the detection of NQR, pulse techniques, offer more 
flexibility over c.w. methods. In principle, a modern pulsed 
NMR spectrometer (with multinuclear capability) which is used 
for the study of NMR of solids can be employed for the study of 
NQR by simply removing the magnet. 

There are distinct advantages in pulsed spin echo methods 
of detection of NQR over the c.w. methods; For example, in c.w. 
methods one often faces the problem of modulation broadening 
and inhomogeneous broadening of NQR line, whereas spin echo 
signals are free from these broadening effects. By following 
the spin echo envelope modulations it is possible to detect 
the small splittings in NQR lines due to dipolar or indirect 
spin-spin interactions [S]. Pulse techniques are better 



20 


suited for the measurement of various types of relaxation 
times. The fast signal averaging facilities in pulsed spec- 
trometers can be used to improve the signal- to-noise ratio in 
the case of wealc resonances. Multiple pulse technic[ues, which 
are commonly used in NMR of solids [95/96] can be used to remove 
the contributions of- interactions which broaden the NQR line 
and thus achieve a narrowed line. The resolution enhancement 
in NQR is / in principle/ facilitates the study of weak inter- 
actions. Using pulse techniques it should be possible to study 
the multiple quantum coherences in NQR also. 

Like in NMR, in NQR also FT techniques have been imple- 
mented [ 97 - 99 ]. FT NQR method has been used for the study of 
spin-spin interactions [lOO] and for the stuc^y of fine struc- 
ture in the NQR spectrum [99]. Details of pulsed and FT-NQR 
spectrometers are available in the literature [8,97-106], 

In NQR spectroscopy the direct detection of quadrupolar 
nuclei with low abundance and/or low quadrupole coupling cons- 
tant, (termed as rare nuclei) is rather difficult. This problem 
can be surmounted to some -extent by using the signal averaging 
facilities available on pulsed NQR spectrometers. However, the 
double resonance techniques offer better sensitivity for the 
detection of the rare nuclei. All double resonance techniques 
take advantage of the large nuclear resonance signal charac- 
teristic -of an abundant species of nuclei 'C (observe nuclei) 
in the sample containing the rare nuclear spiecies 'R' the two 
species being strongly coupled. 
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If the observe and the rare nuclear species are irradia- 
ted with suitable r.f. fields near their resonance frequencies, 
the dipole-dipole and exchange interactions between 'O' and 'R' 
spins can facilitate the observation of the resonance signal 
from '0* with enhanced sensitivity. Different kinds of double 
resonance techniques- have been developed and implemented in the 
literature [107-115]. 

Double resonance techniques can be classified as 

i) Pure quadrupole double resonance (PCDR) techniques, 

ii) NMR-NQR double resonance techniques. 

In PQDR the observe nucleus 'O' and the nucleus to be detected 
'R' are both quadrupolar. The NQR signal of 'R' spins can be 
detected by monitoring the NQR signal of 'O' spins. In NMR 
NQR double resonance experiments the NMR signal of observe 
nucleus 'O' (usually protons) is monitored while the 'R' spins 
are irradiated in the absence of magnetic field. Various double 
resonance techniques that are employed in the literature are 
reviewed [ 116-120]. Here we present a brief discussion of PQDR 
techniques only. There are three types of PQDR techniques that 
are employed in the literature. These are: 

i) spin-echo double resonance (SEDOR) [l07, 108] 
ii) double resonance in the rotating frame (DRRF) [109,110] 

iii) steady-state double resonance (SSDR) [ill]. 

Recently, Ramachandran and Narasimhan [l2l] have developed a 
new DRRF technique based on Off-resonant pulse methods. The 
high sensitivity and low transmitter power requirements of this 
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technique have been demonstrated by detecing Na in NaClO^. 

Employing all the above mentioned double resonance tech- 

V 22339 

niques the resonances of rare nuclei such as, H, Na, K, 
058Vl33l7 2*7 

' Rb, Cs, 0, Al and other nuclei have been detected 
with enhanced sensitivity in many chemically interesting com- 
pounds. In most of these cases the chlorine and bromine nuclei 
have been used as observe nuclei [113,121-124], 

There have been several reports describing pulsed NQR 
double resonance spectrometer in the literature [106,110,117, 
125], Computer and/or microprocessor-controlled pulsed spec- 
trometers have been described in the literature [126,127], 
Recently, Marsimha Reddy et al . [128,129] developed a micropro- 
cessor-controlled pulsed NQR spectrometer. The microprocessor, 
employed in this spectrometer, with proper software is not only 
capable of generating the various kinds of multiple- pulse 
sequences but also controls the entire spectrometer. This 
spectrometer can be effectively used for automatic acquisition, 
averaging and prdcessing of NQR signal data. 

In the next section we present the description of the 
pulsed NQR spectrometer system which has been employed for the 
experimental investigations in the present thesis. 

I.C.3 Description of the Pulsed NQR Spectrometer Employed 

for the Experimental Investigations in the Present Thesis; 

The block diagram of the home-made pulsed NQR spectro- 
meter used in the present investigations is shovm in Fig. I.l. 




Fig.M BLOCK DIAGRAM OF THE MICROPROCESSOR CONTROLLED PULSED 

NQR SPECTROMETER 
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A photograph of the spectrometer system is given in Fig. 1.2. 
Here we present a brief description of various modules of the 
spectrometer. For complete details of the spectrometer system 
the reader is referred to [l06, 128-130 j . The spectrometer 
used in this work can be conveniently operated over a range of 
15 MHz to 40 MHz. In this spectrometer the pulse sequence 
generation is done by an 8085-based eight-bit microprocessor 
system (Microfriend- I, Dynalog Microsystems, Bombay, India). 

The software required for the various kinds of pulse sequence 
generation (see ref, 129) can be stored in the EPROM and it can 
be called on to the on-board user RAM and executed for the pulse 
sequence generation. The pulse sequence is output through one 
of the I/O ports of the microprocessor system and are available 
at the BUG connectors on the side panel of the card cage of the 
microprocessor system. The pulses available at these BNC 
connectors can be used for r.f. gating and triggering of peri- 
pherals. Apart from the ■ microprocessor system an hardware pulse 
programmer [l30] capable of generating a maximum of three pulses 
is also available, in this spectrometer assembly. 

The r.f. signal source is a signal generator (Model 8640B 
from Hewlett-Packard, USA). The r.f. is gated by an r.f. gate 
based on double balanced mixers. The gated r.f. pulses are 
amplified by a home-made low noise tuned power amplifier and 
coupled to the sample circuitry through an impedance matching 
and tuning network. This home-made power amplifier has consi- 
derable power droop and it has poor long-time pov/er stability 
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especially in the experiments involving long train of r.f. 
pulses where large average power is drawn. Hence, in all the 
multiple pulse experiments a commercial broadband ENl (Model 
A-300, Electronic Navigation Industries, Inc., New York, USA) 
power amplifier is used to drive the final tuned stage of the 
home-made power amplifier. This power amplifier combination 
gives a better long time stability of the pulse power and it 
gives an approximate r.f. pulse power output of 800 W. Since 
the r.f. amplitude level, at the output of the r.f. gate is 
unable to drive the ENI power amplifier it is amplified by a 
broadband amplifier (Model Hp-461A from Hewlett Packard, USA) . 
This power amplifier combination can be operated conveniently 
in the range 15 MHz to 35 MHz. A low power commercial broad- 
band amplifier ENI (Model 440LA, Electronic Navigation Indus- 
tries, Inc., New York, USA) is also available for the use at 
higher frequencies. 

The NQR signal induced in the sample coil, by the appli- 
cation of r.f. pulse is coupled to a low noise tuned pre- 
amplifier based on 3N200 FET transistors. The cross diodes 
at the input of the preamplifier protect it from the high power 
r.f. pulses. The output signal at the preamplifier is further 
amplified and detected by a broadband (2-200 MHz) receiver 
(Model 625, from Matec, USA). Although this receiver has the 
provisions for both diode and phase sensitive detections we 
have used, in all the experiments, only the phase sensitive 
detection. In the phase sensitive detection the reference 
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signal is obtained by passing the output of HP 8640B signal 
generator through a delay line (Model 11 - 16 , 0.4 sec, 10 turn, 
variable delay line of. ESC, Electronic Corporation, U.S.A.). As 
the different delay line settings changes the level of the r.f . 
signal, the output of the delay line . is passed through a tuned 
r.f. amplifier before' connecting to the reference input of the 
receiver. The overall recovery time of the preamplifier and 
receiver combination is 10-15 ^^s. In this spectrometer system 
receiver is also gated by a DBM (Model HP 10534A) to avoid the 
overload of the receiver from the r.f. pulse feed-through during 
the 'on' time of the transmitter. 

The detected signal at the output of the receiver can 
either be monitored on an oscilloscope or it can be acquired into 
a signal analyzer (Model SM 2100B, from Iwatsu Electric Company 
Ltd., Tokyo, Japan) . The signal analyzer has two input channels 
through which time domain signals can be acquired. It has, a 12 
bit ADC whose maximum clock frequency is 400 KHz. It can be 
operated in "Auto", "Program" and "BASIC" modes. The acquired 
time domain signals can be stored in any one of the eight blocks 
of the memory and these can be processed using the inbuilt rout- 
ines of the system. The inbuilt routines like average, FFT , IFFT, 
auto-power spectrum, cross power spectrum, autocorrelation and 
cross correlation etc. can be executed by the ‘Function Keys' avai 
lable on the front panel of the signal analyser. In BASIC mode of 
operation, the BASIC Key board is used to write the software 
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program for processing the signals stored in the block memory 
of the signal analyzer. In the Basic program one can use the 
systems routines whenever necessary. ■‘The signals stored in 
the block memory can be displayed on the CRT display available 
in the signal analyzer. This CRT display is also useful for 
BASIC programing. The processed signals stored in the memory 
blocks can be transferred to a mini floppy disk for the perma- 
nent storage. These signals can be read into the block memory 
of the signal analyzer and can be recorded on to a plotter 
(Model SR 6602, from Iwatsu Electric Company Ltd., Tokyo, Japan) 
which has been interfaced to the signal analyzer through the 
GPIB interface facility available on the signal analyzer. In ■ 
this spectrometer boxcar averager (Model CW-1, from Princeton 
Applied Research, USA) facility also available for experiments 
where single point averaging is required. 

Facility for performing double resonance experiments also 
exists in this spectrometer. In the double resonance channel 
a frequency synthesizer (Model HP 3325A from Hewlett Packard, 
USA) is used as the signal source. The output of the synthe- 
sizer is gated by an analog switch (NS-AHOOl 4CD from National 
Semi-Conductor Corpn., USA) for lower double resonance channel 
frequencies. The gated r.f. is amplified by a suitable r.f. 
amplifier depending upon -the r.f. power requirements of the 
double resonance experiment, and coupled to the double reso- 
nance probe through a tuning and matching network. 
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So far we have discussed the principles, applications 
and methods of detection of KQR. In the previous section the 
advantages of the pulsed methods were brought out. In the 
present thesis we shall be mainly concerned with the theoreti- 
cal and experimental studies of NQR signals of spin 1=3/2 and 
5/2 nuclei by pulse -techniques. 

In order to understand the various kinds of pulse experi- 
mental results one needs to have suitable theoretical tools 
for the description of such experiments. In the next section 
we present a brief outline of basic theoretical aspects for the 
description of pulse experiments. 


THEORETICAL DESCRIPTION OF PULSED NQR 

The first experimental observation of pulse responses in 
pure NQR was made by Bloom and Norberg [l31j and Kohn and Herzog 
[l32j. Although the origins of NQR and NMR are entirely diffe- 
rent the transient signals observed in NQR are similar to those 
in NMR. In NMR case, in thermal equilibrium there exists a net 
magnetization and the main interaction Hamiltonian is linear and 
p>roportional to a vector and hence it is possible to explain the 
transient signal formation in terms of simple vector models 
[l33j in several cases. In -the case of NQR in zero field there 
is no net magnetization [134] and the main interaction Hamil- 
tonian is a bilinear operator which is proportional to a second 
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rank tensor. Hence/ it is always not possible to visualize the 
formation of free induction and spin echo signals, in terms of 
simple models. Using a semiclassic al model Bloom et al. [l34] 
have attempted to give a geometrical explanation for free 
induction and spin echo signal formation in the special case 
of axially symmetric, quadrupole interaction in spin I = 3/2 in 
zero field. However, ' no simple geometrical picture is possible 
for explaining the transient signal formation in NQR in the 
case of asymmetric quadrupole interaction and in the presence 
of a weak Zeeman field. However, the formation of transient 
signals in NQR can be studied by considering the quantim mecha- 
nical description of the evolution of the spin system. 


Basic Theory ; 


In this sub-section we present a brief discussion of the 
procedure involved in the description of pulse experiments in 
NQR, in general terms. For the sake of simplicity we consider 
an ensemble of non-interacting nuclear spins with I = 3/2 in 
the absence of a Zeeman field. 


The total interaction Hamiltonian for this case can be 
written as 
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where (^q = *• (1*14) 

a' = 1_ [31^ - + 7] (I^ - ij) J .. (1.15) 

The first term is the quadrupolar Hamiltonian and the second 
terms “H^Ct) defines the interaction between the spins and 
the applied radio frequency (r.f.) field. The response of 
nuclear spin system to r.f. pulse excitations can be calcula- 
ted using the equation of motion of the density matrix. From 
the time- dependent Schrodinger equation, one can arrive at the 
equation of motion for the density matrix [l35]. 


d cr(t) 
dt 


-i [i((t). £r(t)] 


.. (1.16) 


This equation is known as Liouville-von Neumann equation or 
density operator equation. Its formal solution can be written 
• as 

cr(t) = U(t) «r(0) U(t)"^ .. (1.17) 

t 

where U(t) = T exp (- i J ^(t')dt') .. (1.18) 

■ o 

Dyson time-ordering operator T defines a prescription for eva- 
luating the exponential functions in cases where the Hamiltonians 
at different times do not commute i.e. [^^(t')/ ?^(t")] 0. 
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In the high temperature approximation, the thermal equili- 
brium desntiy matri;-: is given by [l36j: 

a(t<0) = {Tr [exp ( exp (- ■^) 

1 

- [Tr(i)J ^ (1 - -^) 

- [Tr(i)j”^ (1 — '^Q) .. (1.19) 

'/6 kT 

where ^ = — = — is the quadrupolar Hamiltonian. In the 

^ ^ kT “ 

right hand side of this equation, the first term being a cons- 
tant it is not affected by any evolution of the spin system. Sq we 
need to follow only the evolution of the second term, proportional 
to reduced density matrix P, which is defined by setting 

o-(t<0) - [Tr(i)_]"^ (1 ^ — P(0)) .. (1.20) 

'fe kT 

where P(0) = "Kq (1.21) 

In the quadruijolar interaction representation defined by 
the transformation operator 

u= exp (i -“c t) . .. (1.22) 

v6 ^ 

the total Hamiltonian can be written as [l36J: 


. . (1.23) 
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where Aw = (w^ -co) 


0 ).. 


^ 'W’ = u la « u 
/6 “ x/i 


-1 


• ( t) 

* ^ dt 


and ^ U 'X^(t) u"^ 


. . (1.24) 


represents the r.f. Hamiltonian, in the interaction repre- 
sentation, obtained after truncating the high frequency terms. 

In the interaction repres&ntation the reduced density 
matrix satisfies the equation of motion 


-1^ ' - - = i [P(t),^] .. (1.25) 

Since, now being independent of time P(t) is given by 

^ r~* 

P(t) = exp (-i'K t) P(0) exp(i^ t) .. (1.26) 

An intense r.f. pulse of width which satisfies the con- 

dition -yH^ Aco/ is referred to as an "hard" pulse. When 
such a pulse is applied then during the pulse the spin system 

f 

evolves under the Hamiltonian only. Folla>/ing the removal 

of the pulse the spin system evolves under the influence of 
the Hamiltonian ^ ‘Hq- Thus/ the Eqn. (1.24) can be used to 
obtain the density matrix of the spin system after a pulse or 


a sequence of pulses. 
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The expression for the signal induced in the coil after 
a pulse or a sequence of pulses can be obtained by taking the 
expectation value of the operator, M, corresponding to the 
magnetization, i.e.. 


<M> = Tr [ ^ m] 


.. ( 1 . 27 ) 


I.D.2 Theoretical Approaches for the Description of Pu.lse 
Responses in NQR. A review ; 

In this section we present a short review of various 
theoretical approaches available for describing the spin dyna- 
mics in NQR. The first quantum mechanical description of the 
formation of the transient signals in NQR was given by Bloom 
et al. [ 134 ]. They used the time- dependent perturbation theory. 
Their method of calculation essentially involves solution of 
the time-dependent Schrodinger equations for the spin system in 
the presence and absence of the r.f. pulses. Using this 
approach they have calculated the single and two pulse responses 
in the case of spin I = 3/2 with axial field gradients in zero 
field as well as in the presence of a weak Zeeman field making 
an angle 'Q' with respect to the symmetry axis of the crystal. 
Following this. Das and Saha [137,138] developed the density 
matrix method for calculating the free induction decay (FID) 
and spin echo signals in NQR. This method involves the matrix 
transformations of the density matrix in .tlie presence of r.f. 
pulses and during the evolution periods. Using this method 



they have calculated the FID and spin echo signals amplitudes 
for nuclei of any general spin I in single crystals with 
axially symmetric field gradients. Pratt [l39j used inter- 
action -representation approach for the description of vhrious 
pulse experiments in the case of spin I = 3/2 nuclei. In this 
method the spin dynamics are follov/ed by matrix transforma- 
tions of density matrix, in the interaction representation, in 
the presence of r.f. pulses and during evolu-tion periods. In 
the above mentioned approaches in obtaining the solution to the 
equation of motion of the density matrix -the explicit matrix 
representations of the operators are used. However, in many 
cases, it is not necessary to consider the entire density matrix 
either because certain elements are irrelevant for the selected 
observable, or because the number of the actual degrees of 
freedom is smaller than the total number of matrix elements. 

In such cases, it is convenient to do the calculations by 
expanding the density operator in terms of a suitably chosen 
set of base operators and to derive equations for their time- 
dependent coefficients [l40]. 

The fictitious spin- 1/2 operator formalism due to Vega 
and Fines [l4l] has been extensively used in the descrption 
of pulse experiments in pure NQR of spin 1=1 systems. Since 
the finctitious spin- 1/2 operators possess simple cyclic commu- 
tation relationships it is relatively easy to evaluate the 
time dependence of the density operator in the presence of r.f. 
pulses and during the evolution periods. Cantor and VJ’augh [l42] 
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have employed this formalism for the analysis of pulsed spin 
locking experiments in the case of nuclei with spin 1=1. 

This formalism has also been employed by other workers in 
the analysis of spin-locking and multiple pulse experiments 
in pure NQR of spin 1=1 [143-146]. Although the fictitious 
spin- 1/2 operator formalism was generalized to higher spins by 
Vega [l47] and by Wokaun and Ernst [l48] there is no report 
where this formalism is employed in pure NQR of spin I> 1. As 
mentioned in Section I, Krishnan and Sanctuary [l6] have used 
the spherical tensor method to obtain the analytical expressions 
for the eigenvalues of spin 1, 3/2 and 5/2 in terms of quadru- 
pole coupling constant and asymmetry parameter. In this method 
they solve the appropriate Liouville equation to obtain time 
dependent polarizations / . They have also given the time evol- 
ution of quadrupolar tensor analytically. Recently/ Bowden 
et al . [149-153] extended the tensor operator formalism/ earlier 
applied to electron paramagnetic resonance [l54] for the study 
of multiple quantum NMR of spin I 1 . They have discussed the 
advantages of this formalism over the fictitious spin- 1/2 
operator formalism and also the differences between their 
approach and the approach used by Sanctuary in the description 
of multiple quantum NMR experiments. 

I.E SCOPE OF THE PRESENT WORK 

This thesis deals with theoretical and experimental 
investigations in pulsed NQR spectroscopy. The tensor operator 
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formalism due to Bowden et al. [149-153] which has been employed 
earlier in NMR has been adapted for the present theoretical 
studies in NQR. The responses of quadrupolar nuclear spin sys- 
tem, to single, double and to a variety of miul tipi e- pulse 
sequences have been experimentally investigated. The tensor 
operator formalism has been employed to analyze these experi- 
mental results. 

In Chapter II of this thesis, we have employed the tensor 
operator formalism due to Bowden et al . [149-153] in the calcu- 
lation of pulse responses in NQR of single crystals. Using 
this formalism, calculation of single and double r.f. pulse 
responses in the case of nuclei with spin 1=1, 3/2 and 5/2 
in single crystals but in the absence of a magnetic field are 
presented. In the case of spin 1=3/2 nuclei, the pulse 
response calculations in the presence of a weak Zeeman field 
are also given. Although we have presented only single- and 
two- pulse responses, especially in the case of spin I = 3/2, the 
results for the evolution of tensor operators given in this 
chapter can be used to calculate the response of any kind of 
r.f- pulse sequence in the absence of relaxation effects. 

Chapter III illustrates the application of tensor opera- 
tor formalism in the calculation of single and two pulse res- 
poses in polycrystalline samples containing spin I = 3/2 nuclei 
in the presence of a weak Zeeman field.. Assuming a configura- 
tion in which r.f. and Zeeman fields are oriented parallel to 
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each other and which in turn are oriented in a random direction 
with respect to the principal axes of the efg we have calcula- 
ted single and two pulse responses and the results are present- 
ed. It is shown that the two-pulse response contains a term 
which is responsible for the spin echo envelope modulations in 
the presence of a weak Zeeman field/ and yields information 

about the asymmetry parameter (r)) of efg. Typical experimental 

3 5 

spectra of Zeeman spin echo envelope modulations of Cl in 
AgClO^ are also given. 

Chapter IV presents the studies on multiple- pulse experi- 
ments in pure NQR. The experimental results of various multi- 
ple pulse sequences in NQR of single crystal of spin I = 3/2 
35 

( Cl) nuclei are given and compared with those obtained from 

23 35 

polycrystalline samples. Results of Na- Cl double resonance 

experiment in NaClO^ powder carried out under the multiple 

pulse spin locking condition are also presented. The experimen- 

T 87 

tal results of various multiple pulse sequences studied on Re 
(spin I = 5/2) nuclei in polycrystalline sample of KReO^ are 
given. Finally, using tensor operator formalism, some of the 
multiple pulse sequences are analyzed in pure NQR of spin I = 3/2 
nuclei and the results are presented. 

Chapter V deals with the theoretical investigation of 
double quantum coherence in pure NQR of spin I = 5/2 nuclei. 
Using tensor operator formalism it has been shown that the 
creation and detection of double quantum coherence in the case 
of spin I = 5/2 nuclei in single crystal can be done by the 
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application of pulse scheme similar to the one employed by 
Hatanaka et al . [155/156] in NMR of three level system. 

SUMMARY 

In this chapter a discussion of basic principles, and 
some applications of NQR spectroscopy are given. A brief out- 
line of NQR detection techniques and the description of pulsed 
NQR spectrometer system employed in the experimental investi- 
gations in the present thesis are also given. The basic theore- 
tical background for spin dynamics calculations and a short 
review of various approaches available for the desorption of 
pulse experiments in NQR are also presented. 
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In pulsed NQR, one is concerned with the study of respons- 
es of guadrupolar nuclei (1^1) to an r.f. pulse or to a seque- 
nce of r.f. pulses. In order to understand the origins of such 
responses of nuclear spins one needs to have a suitable theore- 
tical tool. As pointed out in Chapter I, several theoretical 
approaches are available to describe pulse responses in NQR. 
However, the most convenient method for spin dynamic calcula- 
tions in NMR or NQR is the one based on Liouville-von Neumann 
equation. In this method of calculation one often needs to eva- 
luate the operator expression 

P(t) = exp i-iiit) P(0) exp (ilit) .. (II. 1) 

for different initial conditions and for different Hamiltonians. 
One way of evaluating the above expression is by expressing the 
operators in their explicit matrix representations and carrying 
out the matrix multiplication, which may often lead to lengthy 
expressions. But this problem can be considerably simplified 
when the density matrix is represented in terms of a complete 
orthogonal basis set of operators [l]. Various sets of base 
operators have been used in the description of spin dynamics in 
NMR spectroscopy [l]. 

The fictitious spin- 1/2 operator formalism [2J, has been 
used extensively for the description of spin dynamics in NQR 
in the case of spin 1=1 [3-5], Although these operators have 
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len generalized to higher spins [6/7], there is no report where 
lese have been employed in the spin dynamics calculation in the 
3R of spins 1^3/2. 

Since the quadrupolar nuclei may be present in axial or 
on-axial efg's the study of pulse responses of various spins 
I >l) in both these cases is of importance. Especially, in the 
ase of spin I = 3/2, in order to obtain both the quadrupolar 
nteraction parameters it is necessary to study the pulse res- 
onses in the presence of a small Zeeman field. Keeping this in 
lew, in this chapter we have employed the tensor operator forma- 
ism [S-ll] to calculate the pulse responses in NQR of nuclei 
ith spin I = 1, 3/2 and 5/2 in single cri'stals. Our emphasis 
ill be mostly on spin I = 3/2 case. The results presented in 
his chapter will be used in the subsequent chapters of this 
.hesis for the analysis of pulse responses in a variety of situ- 
(.tions . 

This chapter consists of four sections. In Section II. A 
i brief account of tensor operator formalism is given. Section 
;i.B deals with the calculation of pulse responses of nuclei with 
jpin I = 3/2 and 5/2 with axially symmetric field gradients. 
>ection II. c discusses the pulse response calculation for the 
:ase of spin I = 3/2 nuclei 'with axial field gradients but in 
che presence of a Zeeman field. Pulse response calculation 
:or nuclear spins 1=1, 3/2 with non-axial field gradients 
Ls presented in Section II. D. The last section of this chapter 
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presents calculations for spin I = 3/2 case with non-axial 
field gradient in the presence of a weak Zeeman field. This 
chapter concludes with a siommary . 


II. A TENSOR OPERATOR FORMALISM 

In this section vre present a brief account of the tensor 

operator formalism as used by Bowden et al . [s]. Irreducible 

tensor operators 1" ^ of rank n and order q have been defined in 

4 

the literature [l2j. Irreducible tensor operators for n^ 5 are 
given in Table II. 1. 


Matrix elements ; 

From the Wigner-Eckart theorem the matrix elements of the 
tensor operator T of rank n, order q, are given by [l2J 

Si 


<(lm |t 

i. 


(- 1 ) 


I-m 


( 


Ini 
-m q m 


> <^li 


(II. 2) 


T *T 

where ( _ , ). is the 3i coefficient 

-m q m' 

and 

I II T^ II l]) is the reduced matrix element and is given by 





n : 


.n 


n ! ( 2 1 t n + 1 ) 

(2n) ‘ (2I-n) 1 




1/2 


. . (11.3) 
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Table II. 1. Tensor Operators for after Buckmaster 

et al. (1972) 


T° = 1 ; [3 +1)] ; 

So = '• 1+1 = ^ 5 [^z- t3+'- 


+ I 


^2 1 ^2 
■+ ' J+2 = I ^+ 


[5 ll - {3 Id +1)-1}I^] 

* I I t‘S - Ki+l) - I !, 

(|)^^d[v 4 ]+ 

- 1 ,1.1/2 3 

+ I 

I [35 1(1+1) -25} 1^ + 3 1^(1 +1)^-6I(I+1)] 

+ [n I^[3 Kl+l) +1] I^}, 


= I [n I 2 - Id +1) -5}, I^] 


- ,1.1/2 1 r 3. 

+ ( 7 ) 7 LI 2 ' I+J, 

i I^ 

4 ^+ 


= I (—) [63 I^ + 35{ 3 - 51(1 + 1) } I^ + { 12 - 50 1(1+1) 


+ 15 I^d +1)^} I^] 

z 


^ i ^ ld+l)J ± 

14 {3-1 (1+1) J I, +{12-8 Id +1) +I^'(I +1)^}] 

I (|)^^^ ij [3 I^ +9 I^ +{12-I(I +1) JI ^+{6 -Id +1) 3 ] 

^ ll [5 ± 27 I^ + {24 - I (1+1)}] 

^ (f)^/^ 4 [T^ +2] 




Unit irreducible tensor operators ; 


These are defined as 





Im')> 


r (2n+l) 

LT2I?iy-i 




(II. 4) 


Relationship between unit and non-unit tensor operators ; 


m ^ = JL r (2n-H) (2I-n) ! 2’^(2n)! n 

^ q n ! ^ ( 2 1+n+ 1)1 w q 


Ortho normality ; 

The unit tensor operators obey the following orthonorma- 
lity relationship; 


Tr[ (t J^)' (T 


n 


-q 


q' 




( 11 . 6 ) 


where the adjoint 


(T = (-1)"^ T ^ 




.. (II. 7) 


Symmetric and antisymmetric combinations of the tensors ; 
These are defined as 


T ’^( s) 

^q. 




. 1 


-T [T " - T "] 

/2 -‘J --<J 


. ( 11 . 8 ) 


with q: 0 
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Commuta-bion relationship ; 

The commutation relationship v/hich exists between the 
tensor operators is given by [l3] 


[t^^, tJJ!] = S a(n,n',N,l) {l-(-l) 


n+n ' -N 




] 


NQ 


<(nq| n,q,n' ,q' ^ T 


N 


Q 


. . (II.9) 


where 


i) a(n,n',N,l) = (2N+1) ' (-1) 


1/2 


,N+2I f n. n*N, 11 11 II 


f AJ., AA AN A 

^ I I 


<i|| 

ii) I n, q, n ' , q is Clebsch-Gordon coefficient. 


i> 


iii) { j j is a VJigner 6j coefficient. 


Since the unit irreducible tensor operators form a 'basis set' 
any operator can be expressed in terms of irreducible tensor 
operators. For example, the density matrix P can be expressed 
in the form 


p _ p ^ T 

n,g ^ 


. . (II.IO) 


where the coefficient P'q^/ is referred to as a Pane statisti- 
cal tensor and it is given by 


Pg" = Tr[(Tg")'''p] 


(II. 11) 
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The tensor operators in non-unit form are orthogonal but not 
ortho normal . 

The tensor operators in non-unit form hold for all I. In 
all the calculations in this thesis we use • tensor operators in 
non-unit form. 


II. B PULSE RESPONSES IN THE ABSENCE . OF A ZEEMAH FIELD FOR THE 
CASE OF AXIALLY SYMMETRIC FIELD GRADIENTS 

II.B.l General 

The total laboratory- frame Hamiltonian for an ensemble of 
non-interacting nuclear, spins experiencing nuclear electric qua- 
drupole couplings can be written as 

^ = ifrffciT [31^ - d] + V,U) . . CII.12) 

The first term is the quadrupolar Hamiltonian for the case of 
axially syrnrnetric field gradients, where the symbols have their 
usual meaning [l4] and the second term '^^(t) defines the inter- 
action between the spins and, the applied (r.f.) field. In all 
the calculations in this chapter we assume that relaxation 
effects are not significant. Also we set ti = 1, and express 
energy values in radians/sec. The r.f. pulses are assumed to 
be applied along the x-axis of the crystal (which coincides with 
the direction of the principal field gradient axis V . If the 


r.f. field is applied along the x-axis of the crystal with 
amplitude and frequency oo then 




-2tH^ cos (cot) I = -2 CO cos (cot) 


.. (11.13) 


where co^ = and 

Nov/ ‘^( t) can 


1 is the gyromagnetic ratio of the nucleus, 
be written as 


'K(t) = -^ "3^^ - 2co. cos(cot) I. 

U 1 wX 


.. (11.14) 


where ICq = -p [31 I 


. (11.15) 


Following the convention of Bowden et al . [8],3^(t) can be 
written in terms of tensor operators as: 


^ (t) = T 2oo^ cos(cot) (a) 
^ o 1 


.. (11.16) 


In the quadrupolar interaction representation defined by the 
transformation operator 



.. (11.17) 


The total Hamiltonian can be \/ritten as 
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where = (coq - to) 

^ = U ^ U~^ -1 U — - c :^— ■■■ .. (11.19) 

^ \J6^ 

and 

^1 “ ^ '^l(t) u” = 2to^ cos(( 0 t) U T^(a) U~ .. (II. 20) 

IN# 

is referred to as the r.f. Hamiltonian in the interaction 
representation. 

In the evaluation of the pulse responses we often need to 
know the time evolution of density matrix (thus the time evolu- . 
tion of tensor operators) under the quadrupolar^ r .f . and Zeeman 
Hamiltonians in the interaction representation. We have employed 
the method of “nested commutation relationships” [s] to obtain 
time evolution of tensor operators under the influence of qua- 
drupolar and r.f. interactions. In evaluating the time evolution 
of tensor operators under the influence of Zeeman Hamiltonian, as 
the method of "nested commutation relationships" becomes tedioxos, 
we have adopted the "harmonics of the motion method" [13]. 


II. B. 2 Spin I 


For this case the total Hamiltonian in the interaction 
representation, defined by the transformation operator given in 
Eqn. (II. 17) can be written as 

.. ( 11 . 21 ) 

'16 ° 

where ■Aw= Wq - CO, with Wg = — ^ ■ 9- .. (ii.22) 
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f\J 

can be obtained as follows. Making use of the evolution of 
the tensor operators under quadrupolar interaction (see Table 
II.2), we get 

(t) = 2co^ cos (cot) [t^ (a) cos(cot) + +i T^(s) sin (cot) 

+ T^Ca) -Ir (cos(cot) -1)1 .. (11.23) 

' JTs ■ 

where we have dropped the terms with coherent time dependence as 
non-secular terms. The evolution of the tensor operators under the 
influence of the, Hamiltonian is set out in Table II. 3. 

II.B.2(i) Single-pulse response: 

Let us consider the effect 
of a single r.f. pulse (as shown 
in Pig. II. 1) on an ensemble of 
non-interacting nuclear spins I = 

3/2 with axial field gradients. 

The thermal equilibrium reduced 
density matrix is given by (see 
Chapter I, Sec. I.D.l): 

p- (0) = Tq 


(TT/2) 

X 



Pig. II. 1. Single-Pulse 
Sequence 


Since it commutes with the quadrupolar Hamiltonian we have, the 
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Table II. 2. The Evolution of the Tensor Operators T. (s,a) Under 

oj Q 2 ^ 

a Quadrupole Interaction = — T‘^ for I = 3/2 after 

“ ^ ° 

Bovjden et al . (1986)* 


UT2(s)u“^ 

fV 

= Ths) 

rv -L 

[■| cos (^Qt) + 1 ] - i T^Ca) ^ sin (coQt) ^ 



+ -|r- [cos (o) t)-l] 

x/l 5 ^ 

UT?(s)u“^ 

= -i T^ 

3 v /9 2 

(a) [ ^ sin (co^t) ] + T^Cs) cos (w^t) 



-i T^Ca) ^ — sin (w^t) 

.UT2(s)u"^ 

= T^Cs) 

3 2 

cos (0 t) - i T„ (a) — sin (w.^t) 

W rs/ z w 

UTi (s)u"^ 

= T^(s) 

r*/ i. 

1 /I" [ cos (W^t)-!) -i T^Ca)/^ sinCco^t) 



+ T^Cs) [I cosCcj^t) +-|] 

UT^(S)U“^' 

2 3 

= -i ~ ^' 2 ^®^ cos (w^t) 


U 


= exp (-i 


For the antisimiinetric combinations T^^Ca), replace s by a 

2 123 

and a by s in the above equations, co,. = e'qQ/2. Th T , T_(s) 

3 

and To (a) all commute with the Hamiltonian. 
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able II • 3 . 


The Evolution of the Tensor Operators T^(s,a] 

Si 

the Influence of the Hamiltonian = + 0 )^ [-j 

+ -pr- T^(a)] for I = 3/2 
s/15 


1 -1 
rT U 
o 


'IT O #0 

= T (C0S£ +4) + T^ ■:t ■— (cOS£ - 1) 

-.o t v/s 


v/I 


+ i T, (s) sinE+ i T^ (s) 


3J5 


sin £ 


Jt5;(s)u’ 


T^(s) i [4 C 0 S£ t 1) + T?(s) — ^ (cos£ - 1) 

^ ^ JTs 

- T^(s) 4 (cos£ -1) + i T^ ^ sin£ +i T^ 4^ 

w ^ O O ^ O j— 

-JlT^Cs) ^ Sin£ 


0Tj(a)u“^ = T^(a) (cos£ +4) - T^Ca) i (cos£ -1) 

-T?(a) — ~ (cos£ -1)' -i Ti(a) sin£ 
/Ts ^3 3 


UT^U”^ = T^ COS £ +i t4(s) sin,£ 

IS» O ' rJ JL 


UT?(s)u"^ = T?(s) cos£ + i T^ sin £ 


UT?(s)u“^ = T?(s) ~ (7 coS£+ 3) + T^Cs) — (cos£ - 1] 

10 -.1 

+ To(s) (cos£ -1) + i T^ — ^ sin £ 

2j5 5J5 

+ i T^ sin£ + i Tp(s) sin£ 

5 jYb 


Under 

T^a) 


sin£ 
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Table II.3 (contd.) 

UTj(a)u"^ = T^(a) ~ (3 cos£+7) + T 3 ( a) ^ (cos £ - 1) 

rJ ^ 2 J 5 

- T^Ca) — \zr (cos£ - 1 ). + i T^Ca) sin£ 

^3 JIO 


UT2(s)U ^ 


3 3 1 

To(s) cos£ + i To(s) sin£ 

^2 ^3 ^ 


+ i Tt (s) sin£ - 1 (s) sin£. 


UT 2 (a)u'’^ = T 2 (a) cos£ + i ^ 3 ^^^ ~h sin£. 


+ i T:f(a) sin£ -i T^Ca) -—r sin£ 


3 -1 
UT^U 


Ut2(s)u“^ 


= ^ (1 +4 COS £ ) + T^ ~ (cos £ -1) 

wO 5 .0 

+ i T? (s ) ■ •^-~ sin£ - i (s ) ^^^sin£ 

= To(s) i (l+cos£ ) +T?(s) — (cos£-l) 

2JT5 

- T^(s) (gos£-1 ) +i T^Cs) — sin£ 


UT2(a)u"^ 


T 3 (a) ( 1 + cos ,£ ) + (a) (cos £ - 1 ) 


(a) ~ (cos E - 1 ) + i ( a) — ^ sin £ 
-mI 10 J 2 


where U = exp (-i t) 

^ = 'Is 't 
2 2 2 

TT(a), Tt(s) and T^(a) all commute with the Hamiltonian 
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reduced density matrix in the interaction representation also, 
as 

P(0) = tJ . . (11.25) 


Immediately following the application of an r.f. pulse of width 
‘tyj' the density matrix of the spin system is given by (see 
Table II. 3): 


^ 2 2 
P (t,,,) = cosd + i sinf T. (s) 

w O r\; i. 


where £ = /s u ^ t uj 


. . (11.26) 
.. (11.27) 


Following the removal of the r.f. pulse the spin system evolves 
under the influence of 
by (refer to Table I 1.2): 


^ T^ and the density matrix is given 
v/i ° 


P(t') 


u 


Q 


P(tJ U 


-1 

Q 


where = 


P( t' ) 


f . Aw „2v 
exp(-i — T^) 

v6 


= cos£ T^ + i sin£ [-1 T^(a) sin (Awf) 

O ^ rJ D 

+ T?(s) cos (Awf) - i T^ (a) ' sin (Awt')l 

.. ( 11 . 2 , 8 ) 


v/here f = 

The observed signal is proportional to the expectation 
value of the I^ operator, and it is given by 
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<I (f )> aTr[P(t') (T^ (a)"^)] .. (II.29) 

Making use of Table 11.2 and the orthogonality properties 
of tensor operators we have for the signal in the laboratory 
frame, 

asinf s in (Aco t ' ) ["I cos (cot') + -1] Tr[ (T^(a))(T^(a))’^] 

+ sinCAcot' ) (cos (cot ' )-l) Tr[ (T^ ( a) ) (t^ ( a) )'^] 

+ — ^ cos (Au t ' ) sin(cot ' ) Tr[ (T^ ( s) ) ( ( s ) )'^] 1 

r^l 

.. (11.30) 

which simplifies to give 

xTo ^ 

<(I^(t')^= ^ sinS [sin(Aco +co) t' ] .. (II. 31) 

Using Eqn. (11.28) and Table II. 2 it can be easily shown that 

< I^> = 0, < I^\ =0 . . (11.32) 

These results are in agreement with those reported earlier in 
the literature [14-16]. 

Like in NMR, in NQR also one is generally interested in a 
phase- sensitive detected signal. This can be obtained by consi- 
dering the expectation value of the magnetisation operator in 
the interaction representation [3], The magnetization operators 
(at the resonant frequency) in phase and 90 out of phase with 
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the applied r.f. pulses' in the laboratory frame are given res- 
pectively by 

m' = -cos(cot) T^(a) (in units of -y '^) .. (11.33) 

X N 1 

= -sin(cot) T^a) (in units of -yh) .. (11.34) 

X rO 1 

In the interaction representation these are given by 


M 


cos(wt) U Tha) U"^ 

r>J ± 




(11.35) 


rO 

M 


sin(a)t) U T^(a) u”^ 


i — ^ T^(s) 
2^2 ^ 


(11.36) 


In obtaining Eqns. (II.35) and (11.36) we have made use of 
Table II. 2 and dropped the high frequency terms. 


Now using Eqns. 11.28, 11.35, 11.36 and also mabing use 
of the orthogonality properties of the tensor operators, we have 
the expression for the in-phase signal as 


/3co„ - 

" ~rfe 2 • sinf [sin (A^t')] 


and the 90 out-of-phase signal as 

1 

N ” 4 ' 2 ' sine [cos(Aa)t ')] 


t-t . 

H 


. . (11.37) 


. . (11.38) 


where f 
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In all the subsequent calculations In this chapter we consider the ^ 

35 

expectation value of M only. Typical off-resonant FID's of Cl 
obtained experimentally with our spectrometer from a single carystal 
of NaClO^ are shown in Fig. II. 3. The observed behaviour can be 
explained by Eqn. 11.38. However, the decay of the signal is due to 
relaxation effects which are not included in the theoretical 
discussion. 


II.B.2(ii) Two-pulse response: 

Consider the response 
of an assemble of nuclear 
spins of I = 3/2 to the two- 
pulse sequence shown in Fig. 

II. 2. Making use of Tables 
II. 2 and II. 3 it can be easily 
shown that the density matrix 
after the second r.f. pulse is 
given by 


CV2)^ (^/2) 



Fig. II. 2. Two-Pulse Spin-Echo 
Sequence 

i'iCa)){cos^ t/2 sin(Aot) 

2 

-sin t/2 sin(Aa) ( t-T ) ) } 
cos (&o)t) +cos£ cos { Aw( t-t) ) 

£/2 cos (AW(t-2t))3] 

^ cos( ) 


P(t-T) = sin£ ['/2(T]^(a) ^ + 


ilS 

+ cos£ sin(Aa) (t-T) ) 


2 2 
+ iT^ (s) { cos £:/2 


. 2 
- Sin 


, m2 / 2j, .2 

+ Tq (cos t -sin 


. . ( 11 . 39 ) 
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(b) with resonance off-set AiJ = -8 KHz 
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NOW/ the expectation value of M is given by 

CO 

[sinE cos^ £-/2 cos (Acot) 

+ sin£ cosE cos (Ao)(t-t) ) 

, - sin£ sin^ E/2 cos ( Aco( t-2T) ) j .. (11.40) 

where % » t^ 

t > (t + 2t ) 

UL 

In the above expression, the first term corresponds to the 
remnants of the FID due to the first pulse and the second terms 
is due to the FID of the second pulse while the third term is 
responsible for the echo. 


II. B. 3 Spin I = 5/2 

In the case of spin I = 5/2 nuclei, there are two possi- 
ble transitions, one i^^ith the frequency coq and other with a fre- 
quency 2 Wq . First we consider the transition with frequency. 

II.B.3(i) Irradiation near frequency: 

The total Hamiltonian, in the interaction representation, 
for this case is given by 


= 


A(o 



.. ( 11 . 41 ) 
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where Aco = Wq - w 


CO 


Q 20 


3 2^ 

e qQ 


.. (11.42) 


and K (t) = 2co, cos (oot) U Tj’(a)u'*^. 

1 X /V X 


The evolution of tensor operators under the influence of pure 
quadrupole Hamiltonian ^ = -H ^ is given in Table 12 of ' 
ref. [ 9 ]. Using Table 12 of [9],^. can be written as 


3^^ = 2o)^ cos(wt) T^(a) (9 +16 cosCwt) + 10 cos(2(i)t)) 

+ i ~ — T^(s) (4 sin(a)t) + 5 sin(2cLit)) 
14v/2 

T)f(a) (3 +2 cos (cot) -5 cos(2a)t)) 


OxiTs 

^ t5(s) (2 sin(cot) -sin(2(i)t)) 


S'/? 


+ — - — (3-4 cos(wt) + cos('2cot))] 


3 v/2T0 


.. (11.43) 


'K. - wi [- i| Tj(a) + -4rr T^(a) + — - — T^(a)] .. (11.44) 

1 1 3b X a r 1 /o -1 X 


9 yTs ^ 


3v/210 


This equation is obtained after truncating the high frequency 
terms.. A few selected commutation relations for spin I = 5/2 
are set out in Table II. 4. The evolution of a few selected tensor 
operators under the influence of the Hamiltonian is given in 

Table II. 5. 

Single pulse response: 

Nov7, the thermal equilibrium reduced density matrix in the 
interaction representation is given by 



73 


Table II 


[T^(a) , 
[T^a) , 

[Tha), 
[Tha) , 

/V 

[T^Ca) , 

[Ti(a), 

[Ti(a), 

[T^Ca) , 

^ -L 

[1^(3), 

[Ti(a), 

[Ti(a), 
[ri(a), 
[Tha) , 

[Tj(a), 
[T?(a) , 


.4. A Few Selected Commutation Relations for Spin 
I = 5/2 


T^] 


Ss tUs) 


tJ(s)] = 

TjCs)] = 
= 

T^(a)] = 
T=(a)] = 
Tq] 

;?<=)] = 
T^(S)] = 

T=(a)] = 


T|(a)] 

/] 

tJCs)] 

/(s)] 


- T2(s) - 

- Ths) 

f^j± 


-il 

-Jf -71 ;?(-> 

-Ji 

-\£o tJ(s) 

— p (s) - v/To tJ 

J 2 rJ^ 






-ll T|(a) 


- 2 JI T5(a) -;7 if (a) 
/(a) -J6 T|(a) 

81 . - rr ^4 , . 


7/5 


81 

70 


v/Ts T^(a) - §i ;s - 5JI + 
= Ip /(s) + 75 ^/T5 T.i(s) -v 5 ^ /(a) 


15 


2 So5 


. . .contd. 


K) 



Table I I. 4 (contd.) 


[Ti(a), TjCa)] - 


Sol 


[iJXa), Ti(s)] 


" T^a) 

+ — T^(a) 

+ ---■- T^(a) 

27 



7 yio 

2 vf21 

~ T^(a) 

. 17 „3- . 

•f To(a) 

- -1 T^(a) 


re 

JTo 

/2 



- t 2 (s) - f/| Tt( = ) 


135 

7 


tI 

14 pjo 


[THa) , 

[Tj(a), 

[T^(a) , 
[T^a), 
[T?(a), 

rv -b 


[T^Ca) 


[t? ( a) , 
[T^(a) , 


tJcs)] 

if (a)] 
T|(a)] 

sl‘®>i 

Tq] 

/N/ O 

tJ(s)] 

T^s)] 

T^(a)] 

To(a) ] 


4-, 

To] 




7 7 30 14 


- I/I 

7 V 2 /vo 
18 


^ T 2 (a) + I T 2 (a) 


2 v/n 

25 ^3, > 

r~ ^3^^^ ” 
3\/2 

si " 


T.(a) + 

sJlO ^ ''15 




1 T^(a) 


3J2 


T.(a) - -— X.. 


-/5 T^(s) 

JL 

4/1 T*(s) - 5 if T 


- 1^/2 + 1/1 st‘=> 

- § - 2 ''^ S 3 <^) + if /I 

3 /| T^Ca) _ II /35 T|(a) - § \/7 T^Ca) 

- ^ /e T^(s) + 15 /1 T^(s) 


[T^Ca), Ti(s)] 


75 JTS _ 4/|j T^(s) 


7 


T ^ 2 ^^^ " "T 


^v/ 5 +T^ 


■ T 

io 7 




. . .contd 



Table II. 4 (contd.) 


[T=(a), tJ(s)] 


15 5 4, , , M 15 2 

^ V 3 ^ 


JI| T?(3) -Tt(s) i| 


[T=(a), T^a)] 


[T^(a) , T^Ca)] 


r? T^(a) + ^ if Tf(a) + T^^Ca) 




5 x/S 3!(^) - ^ 'Tl T^Ca) + Ifil Tf(a) 



Table II. 5. The Evolution, of a few Selected Tensor Operators 
Under the Influence of the Hamiltonian 


U T^ U"^ ' = (25 -f 3 cos(£) ) + T^^' ^ (cos(£) - 1) 

^o 28 ^o ^O 

+ i sin(E) (-^ T^Cs) - T^Cs)) 

U T?(s} U“^ = T?(s) i (2 cos(£) + 5 cos(E/2)) 

W -1* W -b / 


- Tts) — — (cos(£) - cos(£/2) 

3 vr 4 

vTs 2 5 4 

+ i sin(£) T^ - T^) 

7^2 " 3'JTo ■" 


+ i sin(£/'2) ( ^2^®^ T 2 ( s ) ) 


U Jo = Jo ^ - Jo M <=“<£> -1) 


+ • i sin(£) (- — T^(s) + 

7 Nfyo 


U Tf(s) U~^ = - T?(s) — (cos(£) -cos(£/2)) + T^Cs) i (5 cos(£) 
-*1 7 714 17 


+ 2 cos (E/2)) + i'sin(E) (- ^ T^ + T^) 

14 77 1475 

+ i sin(£/2) T 2 (s) + I tJCs)) 


'i? = 0). (M T^a) - -4=^ ^T^a) ^ T^a)) 

i ^35-^1 9 JT 5 ^ prrzr>>i 


— i- - I 
\/210 


U = exp- (- i / £ = 780). t^ 




77 


^ O 

P(0) = T 
wo 

Immediately following the application of an r.f. pulse of width 
the density matrix is given by (see Table II. 5): 


^(t ) = + 3 cos£) + — (cosf-1) 

u; 28 wo wo 

- i T?(s) sin£ + i ~ — T^Cs) sin £ .. (11.45) 

7\/8 ■ 2j21 


where H = \/a t^ 


Following the removal of the pulse the spin system evolves under 
the influence of and the density matrix is given by 

^le 


P(t') = — [25 +3 cosEI + — 5 — (cosE-1) 

28 4 J 1 F 5 

+ [- sin£ sin(Aa)t')] 

+ (a) [sinS — ^ sin(Aa)t ' ) J 

6v/5 

+ T^(a) [ s in £ sin(Awt')] 

7v/35 

+ T^(s) [-i sin£ (-^ cos(Acot'))] 

+ Tf(s) [i sin£ ( — — cos(Acot'))] .. (11.46) 

2^/21 


The magnetization operator (in units of^ti) ( 90 '^ out of phase 
with the applied r.f. pulse) in the interaction representation. 
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for this case, is given by 

M = -sin(£ot) UT?-(a)u“^ - f.[ — T^(s) + T^(s)] 

.. (11.47) 

Making use of the orthogonality relationship of tensor operator, 

from Eqns . 11.46 and 11.47 we have for the expectation value of 

■ 

the magnetization operator M 

<M > = 1 sin£ cos(Awt') .. (11.48) 

^ 6 kT 

where, t' = t-t^j 

II.B.3(ii) Irradiation near 2 w frequency: 

u 

The r.f. Hamiltonian for this case is given by 

(t) = 2w, cos(2(0t) T^a) .. (11.49) 

X 1 ^x 

In the interaction representation it is given by 

11 

'K^(t) = cos(2wt) T^(a){9+16 cos(wt) + 10 cos(2cot)J 

+ ■ — — T?(s){4 sin(Qt) + 5 sin(2o)t)} 

14{2 

- — T?(a){3+2 cos(wt) - 5 cos (2 cot) 3 
9 VTS ^ 

- — ^ (s) {2 sin(a)t) -sin(2cot)} 

+ — — — T^(a) { 3-4 cos ( o)t) t cos ( 2tot) 3 ] 

3\[2To ■ 


.. (11.50) 
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After dropping the high frequency terms we have 



f-" 7 


_5 

9v/r5 


•T J ( a) 
X 


1 

3\/2r0 



(11.51) 


The evolution of a few selected tensor operators under the influ- 
ence of is given- in Table II. 6. 

II.B.3(ii)a Single pulse response: 

The thermal equilibrium reduced density matrix in the 
interaction representation is given by 

P(0) = T^ . . (11.52) 


Immediately following the application of an r.f. pulse of dura- 
tion 't^^* the density matrix is given by (see Table II. 6): 


P(t, ) = 4 (4+3 cos£ ) + ^ T^ (cos£ -l) 

w 7 ^o _ . ~o 

- i T^(s) sin£ - i Tt(s) sin £ .. (11.53) 

7 


where £ = \/5 t^ 


Following the removal of the pulse/ the spin system evolves 

under Xw (see Table 12 of ref. [9]) and the density matrix 
/6 

is given by 
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Table II. 6. The Evolution of a few Selected Tensor Operators 
Under the Influence of the Hamiltonian "U. 


U T];(a) U"^ 


U T^Ca) U"^ 


= T5;(a) ^ (19 +16 cos(£/2)) 

- T?(a) (cos(f/2) - 1) 

1 9'JTS 

- T^Ca) -4=: (cos(£/2) - 1) 

^ 3V210 

- i sin(£/2) Tp(a) + — — T^(a)) 

= T^Ca) — (1 - cos(£/2)) + 

1 SSnTIS 

+ T^(a) (14 + cos(£/2)) 

+ T^(a) (cos(£/2) - 1) 

1 srn; 


+ i sin(£/2) (T^(a) + — ~r T^(a)) 

2 3^2 5m ^ 


U T^(a) U"^ 


2 -1 
U T^ U 
o 


- T^a) (cos(£/2) - 1 ) 

1 7n/2T0 

+ T?(a) (cos(£/2) - 1 ) 

1 sm 

+ T^(a) (11 + 10 cos(£/2))) 

+ i sin(£/2) (^ iha) + Tf(a)) 

■i T^ (4+3 cos(£)) + -~r- (cos(£) 

7 o Vl^5 o 


+ i Sin(£) ('^ tJ(s) + T^(s)) 

' ^ 'TTTo ^ 


U T^(s) u"^ 


= T^(s) rj (2 cos(£/2) Hf- 5 cos(£))4- 

+ Tt(s) — ~ (cos(£) - gos(£/2)) 
3^4 

^ i (^15 ^2 ^ ^ 

- i s in ( £/2 ) (~ tI ( s ) + ^ T2 ( s ) ) 


. . .cont( 



Table II. 6 (contd.) 


4 -1 

U U 

o 


U tJ(s) u“^ 


- 1) + tJ i (3 + 4 cos(£)) 

+ i sin(2.) (—2^ T?(s) + T^ i 

7V7 ^ 7 1 

2 


^ T’(s)) 

T^(s) (cosCe) - cos(£/ 2 )) 

+ T^(s) (5 cos(£,/ 2 ) + 2 cos(£)) 

+ i sin(£)(-i^ ^ tJ) 

7 \A 210 o 7 o 

+ i sin(J/2) T^(s) + ^ tJ(s)) 


U T^Cs) u“^ = cos(£/ 2 ) 


+ i sin(£./ 2 ) (- T^(s) + T^(s) - ^ 


U-T 2 (a) U”^ = T^Ca) cos(£/ 2 ) 


+ i sin (£/ 2 ) ( ^ Tha) + ^ 

7 'JJO -‘■ 


+ ^ T^(a) - -4- T^(a)) 

\ri 1 S'/e 3 


Tj(a) - ^ T 3 
Sn/" 2 ^ 2^6 3 

2 „5 




U T^Cs) U"^ 


U T^Ca) U“^ 


= T^Cs) cos(£/ 2 ) 


4i sln(£/2) .(-^T^s) * = ) 

= T^Ca) cos('£/2) 

4isln(a/2) 

* ffi 




u = exp (-i V 


— 4=:r T^a}) 

3'^210 


S t„ 



02 


P(t‘) = [sin£ sin(2A6jt') {- ^ 

- -1— T^(a)} 

15^/7 


+ i sin£ cos(2Awt') {- T^(s) - T^(s))] 

7 -1 


.. (11.54) 


where t‘ = t.-t^J. 


The laboratory frame magnetization operator, for this case is 
given by 


M = -sin(2Qt) T^(a) (in units of 
In the interaction representation it is given by 
M = -sin(2ut) T^(a) 

c; - X[t^(s) ^ — T^(s)] .. (11.55) 

^ 28^2 12 '/?^^ 


The signal in the interaction representation is given by the 

W 

expectation value of the magnetization operator 


<M^>a Tr [P(t') 


Making use of the orthogonality property of the tensor operators 
we have 


rsj 

<M > 


6-fe — 


(2Acjt ' ) 


. . (11.56) 
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II.B.3(ii)b Two-pulse response: 

We consider here the response of nuclear spins with spin 
I = 5/2 to the two- pulse sequence shown in the Pig. II. 2. MaJcing 
use of Tables II. 6 and 11.12 of [9j and proceeding as in the 
single pulse case we obtain for the signal following the second 
pulse as 


2 ^ . - 

[sin£ cos^ S/2 cos(2A(ot) 

+ sin£ cos £ cos ( 2Aci)( t-T ) ) 

2 

- sin£ sin £/2 cos ( 2Aa)( t-2T) ) J .. (11.57) 

where £ = 'fs t^^ 

T t 

UJ 

t > (t + 2 t^^) 

In Eqn. 11.57 the first term corresponds to the remnants of the 
FID due to the first pulse and the second term is due to the FID 
of the second pulse while the third term is responsible for the 
echo. These results are in agreement with those available in 
the literature [l4]. 
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II *C PULSE RESPONSES IN THE PRESENCE OF A WEAK ZEEMAN FIELD FOR 

NUCLEI WITH SPIN I = 3/2 IN AXIALLY SYMMETRIC FIELD GRADIENTS 


II.C.l Transitions for a Nuclear Spin I = 3/2 in the Presence 
of a Weak Zeeman Field 

If we apply a constant magnetic field in a direction with 
polar coordinates 0, -0 in the 'efg' principal axes system then the 
Hamiltonian corresponding to the interaction of Zeeman field with 
nuclear spin system is given by 

^ = - I H f cos© sin0 + I sine sin©' + I cos©) 

z o ‘ 

.. (11.58) 

The presence of the magnetic field (for ^ g g;Q) removes the 
degeneracy between states m = +3/2 and first order perturbation 
theory gives [is] 


E+m = A [3 m^-l(I+l)]+ m cos© 


.. (11.59) 


corresponding to the states 'i' and ? / respectively, with 


A = 


2 ^ 
e qQ 


41(21-1) 


.. (11.60) 


The states with m = +l/2 have non-zero off-diagonal elements 
c;^l/2j Hence, diagonalization of sub- 

matrix for m = +1/2 leads to tlae zero order mixing of the states 
'^+1/2 ^ l/2 yields two new states 
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’*^ 4 - “ ^+ 1/2 ^- 1/2 


where 




'^_l/2 sin a - '^^1/2 ^ 


tan a - L (|_|-5 J 

r 2 .1/2 

f = [ 1 + 4 tan. q] 


.. (11.61) 


.. (11.62) 


The energies of these states E_j_ are given by 


E 


_^ = A [3/4 - 1(1+1)] + f/2 co^ cosQ 


where 


.. (11.63) 

.. (11.64) 

Now, four transitions with Am = +1 are possible between 
the mixed states and with states '^^ 2/2 frequencies: 


0 ) ^ = 
a 


CO. 


00 / “ 
a 


= 6a -- 


6A - 


6A + 


%' ~ 


= 6A + 


3-f 

2 

“0 

cos© 

3+f 

2 

“0 

cos© 

3-f 

00 

cos© 

2 

0 


3+f 

2 

^0 

cos© 


.. (11.65) 


Thus, the original line is split into- a symmetrical quadruplet, 
the outer pair of lines is referred to as the (X doublet and the 
inner pair as the 8 doublet. It should be noted from Eqns. 

11.63 and 11.65 that the Zeeman pattern is independent of azi- 
muthal angle /. The energy level diagram for the above situation 
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is given in the Fig. II. 4. 

I I . C . 2 Evaluation of Pulse Responses 

In this section, we assTome in our analysis that a weak 
Zeeman field is applied in a general (G, jZ5) direction in the 
principal axis system,' and that the r.f. field is applied para- 
llel to the principal x-a;<is. Since, for the case of axially 
symmetric field gradients (r? = 0) the Zeeman pattern is indepen- 
dent of ^ (see Eqns . 11.63 and 11.65) we may set the azimuthal 
angle j6 equal to zero without loss of generality. 

Now, the Zeeman Hamiltonian, in the laboratory frame, 
corresponding to this situation is given by 

"K = - [i sinG ' + I cosG] .. (11.66) 


v;here being the strength of the magnetic field. In 

the interaction representation it is given by 


= u u' 

2 Z 


^ [sinG {2/5 T];(a) 


tJ (a) } -cosG tJ] 
-1 -O'* 


(11.67) 


In arriving at Eqn. 11,67 we have made use of Table 11,2 and 
truncated the high frequency terms. 

In the presence of weak magnetic field we consider the 
response of nuclear spin system to resonant r.f. pulses (i.e.. 



4“ — ; — + 3/2 COq Cos 0 



nergy level diagram for nuclear spin 1 = 3/z in the presence 
weak Zeeman Field Hq* 
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Aid =0). So in the interaction representation the total Hamil- 
tonian for this case can be v/ritten as, 

'ii = '5 + if .. (11.68) 

J- 2 

The evolution of a few selected tensor operators under the influs' 
ence of the Zeeman Hamiltonian is set out in, Table II. 7. Table 
II. 8 gives the commutation relationships of the Hamiltonian with 
the tensor operators. In all the calculations', we neglect the 
effect of the Zeeman field in the time interval (O^t^tjjj) 
during which the pulse is on, assuming that the magnetic field 
causes negligible precession of the spin system during applica- 
tion of r.f. pulse. 

II.C.2(i) Single pulse response: 

The thermal equilibrium reduced density, matrix, in the 
laboratory frame for this case is given by (see I.D.l), 

p(o) = T^ + 2__ [T::(a) sine - t;: cose] .. (ii.69) 

However, if the magnetic, field satisfies the condition ' 

then, to a good approximation, the initial density matrix can be 
taken as : 

P(0) = T^^ 

w O ■ . .. , , 
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Table I I. 7 


Uz5i<=>'^z^ 






The Evolxition of a few Selected Tensor Operators 

Under the Influence of the Hamiltonian for I = 3/2 

2 


T?(s) [(J^) cosA + (f~) COSB] 


+ i T^(a) [ (|^) sinA + sinS] 

2 f -F ^ _ 1 ^ 

if [cosA -cosB] 

2 ( f 

+ i T 2 (s) [ (sinA-sinB)] 


T^Ca) cosA + (-1^) cosB] 

+ i T^(s) [ (|^) sinA + (f^) sinB] 


2 1/2 

+ T 2 (s) [-^ (cosA - cosB) ] 


2 (f^-1)^^^ 

+ i T 2 (a) ■:^ (sinA -sinB) 


= T 


2 ( -F^ 1 ^ ^/^ 

^(s) 2 f [cosA -cosb] 


2 f -F^ 1/2 

+ i ^i(a) (sinA - sinB) 


^2 t cosA + (f^) cosb] 


tiT^Cs) [(|^) SinA + (^) SinB] 


. . .contd. 



Table I I. 7 (contd.) 


90 




2 _ 1/2 

Ti(a) [cosA -cosb] 

2 ( f ^ — 1 ') 

+ i T. (s) rrr- (sinA -sinB) 

i. . z r 

+ T^Cs) [(^1^) cosA + (^§|) cosB] 

+ iT 2 (a) [ (^^) sinA + sinB] 


2 -1 


where 


exp ( -i ) 

oa [sine { I T^a) -T^(a) > - cose tM 

o'- 5 ^ 

[1+4 tan^e]^"^^, A = (•^^) cose Ci) t 

^ O 


cose t, w = YH 
2 o o o 



✓ ' 


Table I I. 8. Commutation Relationships of Irreducible Tensor 


Operators T!i(s,a) for I = 3/2, v/ith the Hamiltonian 

^^'7 = [j3{|. T|(a) -T^Ca) ] - aT^] = oF 

Z O ^ ^ 5 rul rip- O 
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Immediately following- the application of an r.f. pulse of width 
'tm' the density matrix is given by the Eqn. 11.26. Since we 
are considering' the pulse response/ on resonance (Aw= O), the 
density matrix following the 'removal of the pulse evolves under 
the influence of the 'Zeeman Hamiltonian and it is given by 

(see Table II. 7) : 


P(t') = cosH, + i sin£ cosA + (^'^) cosB} 

+ i Ti(a){(-^) sinA + (^) sinB} 
2 ( f 1 ) 

+ T2(a) (cosA-cosB) J 

2 C f 1 ^ 

+ ijr 2 (s){-^ 2 ^— — (sinA-sinB) ) 


1 . (11.70) 


where A = 3+f/2 cos© 

B = 3-f/2 cos© (o^f 


- .. (11.71) 


Now the signal in the interaction representation is given by 


<M^> a Tr[P.M^] 


J3 


0) 


Q 1 


sin£ [ ('^^) cos A + 


4 kT 


cosB] 


.. (11.72) 


Similarly/ in the laboratory frame the operator corresponding to 
the Y-component of the magnetization (in quadrature with the 
applied r.f. pulses) is given by 

My = sin(cat) ly = i sin(6ot) T^(s) (in units of Tti) 

. . (11.73) 
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In the interaction representation it is given by; 




. (11.74) 


<V 


f3 


0) 


^ - «in£ [(^) sinA + (^) sinB] 


4 kT 2 


. (11.75) 


These results agree with those presented earlier in the litera- 
ture [l5/16j. Also, it can be easily shown that 



These results indicate that in the presence of a weak Zeeman 
field, a crossed-coil experiment is possible in NQR. Now, the 
FID signal is modulated with the frequencies which correspond to 
the splitting of the steady-state pure quadrupole line due to the 
magnetic field. For example, at a particular orientation of 
magnetic field i.e., with Q = 54*^44', the steady-state splittings 
of the aoc' lines coalesce into a single pair at the original 
pure quadrupole frequency [l4j. For this orientation of the 
magnetic field we have 

f = (1 + 4 tan^ Q) =3 

From the Eqn. 11.72, for this orientation of magnetic field it 
can be seen that the FID gets modulated by only one frequency 
corresponding to and is given by 

3 0)^ cos ( 54*^44 ' ) 
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. ' 35 

The experimental results of the •modulation of cl FID signal in 
the single crystal of NaClO^ when the Zeeman field orientation 
corresponds to e = 54*^4 4' are shown in Fig. Ii.5. From this 
figure it is clear that the FID is being modulated by only one 
frequency/ which is equal to 3 cos(54°44'). 

II.C.2(ii) Tv;o pulse response: 

Again considering the two- pulse sequence given in Fig. II. 2, 
the density matrix after the second pulse can be written as 

P(t-T) = P(0) u-^(x) 

.. (11.76) 

where ^ “ exp(-i tu^) .. ( 11 . 77 ) 

n (x) = exp(-i^^'c) .. (11.78) 

z ^ 

r-J ^ 

'K and 'i^ry are given respectively by Eqns . 11.24 and 
1 4 

11.67. 

NOW/ making use of Tables II. 3 and II. 7 we obtain, after some 
trigonometric manipulations/ for the expectation value of the 

magnetization operator 

“ 4^ T [sinE cos E /2 i ( 2 ^“^ cos (■^^ cosG m^t) + 

(^) cos (^ cos© oo^t)} . . (i) 

. . . .contd. . 
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(b) In the presence of a Zeeman Field. HgCi 15-9 Gauss. 
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+ sine cose { (^) cos cosQ u^Ct-T)) ' 

+ cose('^^ cosS co^(t-T;})3 ..(ii) 

2 'F— 1 2 '^4»'F 

- sine sin £ /2{ cos cose (ji^{t-2x)) 

+ cos cose to^(t-2T) ) J - . (iii) 

2 

- sine sin^ e/2 [ (^-^) { cos (•! cose co (t-2T)). 

2f f 

cos (“ cosG CO t) 
z o 

f 3 

+ cos (-r cose (0 (t-2T)) cos (t: cose CO t) 

2 0 2 o 

- cos (ry cose w t) cos (■?■ cose ^ ■t)}]1.,(iv) 

.. (11.79) 

Our results are in agreement with those reported earlier in the 
literature [l5,16]. The first term corresponds to the remnants 
of FID due to the first pulse, the second terms gives the FID due 
to the second pulse while the third term contributes to the echo 
at 2 1 and the fourth term -is responsible for the echo envelope 
modulation as a function of t . 

Experimental recording of two-pulse response, in the pre- 

35 

sence of a weak Zeeman field, obtained from Cl in single crystal 
of NaClO^ is shown in Pig. II. 6. 
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II. D PULSE RESPONSES IN THE ABSENCE OF A ZEExMAM FIELD FOR THE 
CASE OF NOH-AXIALLY SYMMETRIC FIELD GRADIENTS 

This section deals with the calculation of single- and 
two-pulse responses of nuclear spins 1=1 and 3/2 with non- 
axial field gradients, we first consider tlie 1=1 case. 

II.D.l Spin 1=1 

Since the asymmetry parameter (V) mixes the basis functions 
whose m values differ by +2, the eigen functions of are 

linear combinations of these basis functions [is]. The eigen- 
functions, energy levels and the transition frequencies for spin 
1=1, nucleus with non-axial field gradients are given in 
Fig. II. 7. Now, as there are three transitions possible we 
discuss the observation of signals corresponding to each of these 
transitions separately. 

The evolution of tensor operators under the quadrupolar 
Hamiltonian = 

In the interaction representation defined by the transformation 
operator 

U = exp (t i — (tJ + r)/\f3 T^(s))) til. 80) 

^ .,2 

the total Hamiltonian is given by 

X = ^ (T^ tLs)) 


x/3 




Q 


pT t Ti/\f3 Tf(s)] are given in Table II.9. 


(11.81) 





6 

IL 
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Table II. 9. The Evolution of Irreducible Tensor Operators Under 
the Influence of Quadrupolar Hamiltonian 

= — w,, [t^ + Ti/\f3 T^Cs)] for sp-in 1=1* 

UT^U ^ cos (f.t) + i T?(a) \f2 sin (f.t) 

r^O -rviO Jl 1 

UT^(s)u'"^ = T^(s) cos (f 2 t) - i T^(a) \/2 sin (f 2 t) 

UTT(a)u"’^ = TiCa) cos (f^t) - i T?(s) 'f2 sin (f^t) 

O rv#J. O 

2-12 1 1 
UTi(s)U = T- (s) cos (fot) - i T-, (a) — sin Cf^t) 

^1 3 v/2 ^ 

UTT(a)u”^ = T?(a) cos (f-t) - i t}(s) — sin (f«t) 

wl ~i ^ ^ z 

UT^(a)U ^ = To(a) cos (f.t) + i T^ ~ sin (f^t) 

~2 1 ^o ^ 1 

where U = exp (-i ’^^gt) 

2 2 

T and T^Cs) commute with the Hamiltonian 

^ O r>J Z 

fl = Irjcog , f2 = Wg (1 -T)/3), f^ = Og (1 + 7?/3) 

3 2 

Wq ~ ^ 

*These results can also be obtained from Table 1 of ref. [lo] 
by setting 1i = 1, & = 0. 





(11.82) 


where Aco = “q > 

3 2 

“q = 4 ^ 

II.D.l(i) Irradiation near frequency; 

The r.f. Hamiltonian for this case can be v/ritten as 

= 2to. cos (f^t) U T,(a) U” .. (11.83) 

Using Table II. 9 and dropping the high frequency terms 

we get 

" “l •* 

The evolution of tensor operators under = w^Tj^(a) can be 

obtained from Table 3 of ref., [s]. 

II.D.l(i)a Single-pulse response: 

The thermal equilibrium reduced density matrix is given 
by 

P(0) = (T^ + ^/v/3 T^(s)) .. (11.85) 

Immediately follov/ing the application of an r.f. pulse of width 
'tyj' the density matrix is given by 



p(\,) = [ !■ (1-77) + (1 + 77/3) cos£: J 


3 

4 

i T^(s) sin£ 


J 3 


(1 +77/3) 


+ T^Cs) [- ^ (1 -77 ) + ^ (1 +77/3) cos£ ] 


.. ( 11 . 86 ) 


where E = 20 )^ 


Following the removal of the r.f. pulse the spin system evolves 


under the influence of 
matrix is given by 




r3 


(T -h r )//3 T^(s)) and the density 


P(t') = ( T (1-77) + I (1 + 77/3) cos£ ) 


J3 


/3 


+ ^2(3) [- — (1-77) + — (1 + 7 ?/ 3 ) cosH:] 

+ T^(s) [“ i — (1+77/3) sin£ cos (Af^t')] 
+ T]^(a) [- (1 + 77/3) sin£ sin (Af3t')] 


2 <2 


. ( 11 . 87 ) 


with t' = t-t 


UJ 


The magnetization operator for this case is given by 


= 

-sin(f3t) T^(a) (in units of rii) 


4 -^ 

M = 

-sin(f^t) U T?-(a) U"^ 


X 

0 rsil 



-i I/V2 T?(s) 

.. (11.88) 


Now, the expectation value of the magnetization operator is 
given by 



I — tn/3) . sin£ cos (Af^f) .. (11.89) 


AlsO/ making use of the orthogonality property of the tensor 
operators, from Eqn. 11.87 it can be easily shown that 


<M > 
y ^ 

1 

<T^(a)> = 0 



1 

II 

= 0 

.. (11.90) 


II.D.l(i)b Tv 70- pulse response; 


Let us consider the effect of two pulse sequence shown 
in Fig. II. 2 on the nuclear spin 1=1 system. The density 
matrix after the second pulse can be written as 


P(t-T) = Uj^tt^) Uq(x) 0(0) URhtJU'Cx) 






rs/ 

where Up^(tuj) = exp (-i tjj^^) 

Uq(t) = exp (-1 ~ Au)^) 

rs/ 

'^1 is given in Eqn. 11.84. 


.. (11.91) 


.. (11.92) 


Making use of Table II. 9 of this thesis and Table 3 of ref. 
[8] we obtain for the expectation value of the magnetization 
operator, 



I UH 


A [sin£ cos^ E/2 cos (Af 3 t) 

+ sin£ cos £ cos ( Af ^ ( t-T) ) 

2 

- sin£ sin E/2 cos (Af ^ (t-2T ) ) J 

.. (11.93) 


where t ;$> t^J 

t > ( 'c + 2 t^j) and £ = 2a)^t^, 


II.D.l(ii) Irradiation near ±2 frequency: 

if we wish to observe the signal corresponds to frequency 
±2 we need to apply the r.f. pulse, with this frequency about the 
Y-axis of the crystal. The r.f. Hamiltonian for this case is 
given by 

= -i cosCf^t) 

- -icj, T^s) . . (11.94) 

The evolution of tensor operators under the influence of the 
Hamiltonian = -iw.T;r(s) can be obtained from Table 4 of 

ref. [s] . Making use of the Table II. 8 of this thesis and 
Table 4 of ref. [Sj the expectation values of the magnetization 
operators in the case of single and tvro pulses can be written as, 

II.D. l(ii)a Single-pulse response: 

<( My> ~ § sinf cos(Af 2 t') .. (11.95) 

= 0 . 


where t' = t-tyj ; an.d = ^M^)- 
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II.D.l(ii)b Two-pulse response: 

= f It [sin£ cos^ £/2 cos(Af2t) 

+ sinE cosE cos (A f 2 ( t-r) ) 

2 

- sin£ sin £/2 cos ( Af 2 ( t-2i:) ) ] 

.. (11.96) 

<\> = <^2) =0 

where t 

t > (T +2t^^) 

II.D.l(iii) Irradiation near frequency: 

In order to detect the signal corresponding to frequency 
f^ we need to apply the r.f. pulse with this frequency about the 
Z-axis of the crystal. 

The r.f. Hamiltonian corresponding to this situation is 

.. (11.97) 

In the interaction representation it is given by 

= -2co, cos(f-t) tJ; ^ .. (11.98) 

JL J, J.wOX'^O 

r-* 

The evolution of tensor operators under the influence of 
can be written as (see ref. [s]). 



I UU 


exp Tq(s) exp (i 

= tJ^(s) cos(q.co-t ) + iTjJ(a) sinCq.co^t ) .. (11.99) 

^ JL UJ UJ 

T^i(s)'s will also have a similar form. Using these relationships 
wq 

and Table II. 9 the expectation value of the magnetization opera- 
tor after single and , two-pulses (see Pig. II. 1) can be obtained. 

II.D.l(iii)a Single-pulse response: 

u 

<M^'> = I ^ (| n) § sinE cos(Af^t') .. (II.TOO) 

Also <M^> = <My> = 0 ; t' = t-tjjj. 

II.D.l(iii)b Two-pulse response: 

CO ^ 

1 [sin£ cos E/2 cos(Af^t) 

+ sin£ cos E cos ( Af ^ ( t-t) ) 

- sin.E sin^ E/2 cos (Af ^ ( t-2T ) ) ] .. (II. 101) 

<M^> =<My> = 0 

where ; t>(T + 2t^j) . 

II. D. 2 Spin I = 3/2 

In spin I = 3/2 case the asymmetry of the efg, does not 
remove the degeneracies of the energy levels. The basis lm^)> 
functions are no longer eigen functions. The asymmetric compo- 
nent of the efg mixes the 1-3/2)* with the \+l/2)> function and 
the l-l/2)>’ with the )+3/2y function [15]. Due to the Kramer's 
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degeneracy even with the asymmetric efg's only one transition 
is possible. The transition frequency between the two degene- 
rate levels, is now given by [iS]: 


COg = Oq Cl + 71 /3) 


.. (11.102) 


where 


The evolution of tensor operators under the Hamiltonian 

(jD ^ O 

^.= — ^ (T'^ + 7)/V3 T^Cs)) is given in Table II. 10. The r.f. 

Q vf6 

Hamiltonian in this case is 


'^1 = 2w^ cos(cOgt) T^(a) 


.. (11.103) 


where co' = w(l 


In the quadrupolar interaction representation, the total Hamil- 
tonian for this case is given by 


^ ^ + n/J3 T^(s)] + 'S. 

v/6 --2 1 


.. (11.104) 


where 


= 


2w; cos(w* t) UTj'(a) u"^ 
i q r-f 1 


t CO ^1 ^i3|i2l_ ^ ( 2 z i ; ^ ) ^3^,) 

1 ^ 3^2 5 ~1 ^ ..1 


- - T^(a)3] . . (II. 105) 



Table II. lO. The Evolution of Tensor Operators T^ Under a 

Cl} 0 O 

Quadrupole Interaction [t + T^(s)] 




for spin 1 = 3/2 


v/6 


-^o 


1 -1 
UT U 
w O 




15r± . 4n- . icos .. 2J2 


15 


3n/3 


- JTq 7 )^ (cos (wt)-l) 4- i T^ (a) sin(at) ] 


UT^ (s)U ^=^[T^(s) 33 { ( 3 - 71 ) ^ cos (wt) +2(3 + 271 ^') J 

+ T?(s) — ^ { (2 + 71 ) ( 3 - 71 ) (cos (wt)-l) } 

3\/l5 

- 13 ( 3 ) i 71(3 - 71 ) (cos (cot) -1) 

-iT?(a) — i— (3 - 7i) sin(wt) ] 

~ 3\/2 

UT^(a)u“^ = -3 [T];(a) ~ {3,/ cos(cot) +2 ( 3 + 271 ^- 3 ti) } 


+ T?(a) (cos (cot) -1) 

sv/Ts 


- T^Ca) ^ (3 + 7 l) (cos (cot) - 1) 


- 1 


its) 

3/2 


2 -1 
UT-U 

rjO 


t ^O 1 ^ +71^(cOS (cot) -1)} 


- T^(s) (cos(cot) -1) + i T^(a) sin(cot)] 


nTs 



able II. 10 (contd«) 


Ti(s)u"^ 


[ T? (s) f cos ((j>t) - 1 T^(a) (3+'/]) sin (cat) 

^ rJl D 


n/2 


iT, (a) (2-77) sin(a)t) +i T,(a) 7}sin(a)t)j 


2 -1 
Ti(a)U 


F t cos (cot) f - iT^(s) ( 3 - 77 ) sin(ci)t) 


3, , '/2 




- iT^ (s) - (2 4-7)) sin(a)t) + iTo(s) -^77sin(cat)] 


2 -1 
To(s)U 

fj Z. 


[ T2^s){ + (cos (o) t) -1) ) - T^ (cos(wt)-l) 


iTo(a) — - sin (cot)] 


2 -1 
T,(a)U 

/Nf Z 


1 r „2 


2 . 


— r T^(a) f cos (cot) -i T^(s) — sin (cot) 

jj L ^2 ^ -3 


+ i 4 - 7 ^ sin(cot) - i ^ sin(cot)] 

r^o J 2 3v^ 


fT? (s)u“^ 

/xT Z 


^ [ tJ(s) ^ { 3(4 +477 + 77^) cos (cot) 4- (18 +777^-1277)} 


+ ihs) — — { (2+77) (3-77) (cos (cot)-!)} 

lOv/Ts 


_ Tt(s) (cos((0t)-l) 

2\/r5 


- iT:(a) i (2 +77)f Sin(wt) }] 

2 \ttO 


JTj(a)u“^ 


— [ T^(a) { 3 (4-477+77^) cos (cot) +(24+777^-1277)} 


+ Tt (a) - (2-77) (3+77) (cos (cot) -1) 


2^5 


- T4(a) (cos(cot)-l) -i T?(s) (2-77)f sin(cot)] 

2 VTo 


I t u 


Table II. 10 (contd.) 


3 -1 
UT U 

r-J O 


•4 [ ri { 15 / -f r}^cos(o)t)-l) } +T^(s) (cos(cot)-l) 

p-i ID J 


- 4 ■~T=“ (cos(a)t) -1) -i T^(a) sin(wt)] 


UT 2 (s)u^= 4 [T 2 (s) {j>^ + (cos (w t) -1) } - (cos (wt) -1) 


^ ^ (cos(cot) -1) -i T 2 (a) sin(a)t)] 


+ — =- 

s/l5 


UT 2 (a)u'” ~ j I- ^2^^^ cos (cot) - iT 2 (s) -y sin(cot) +i ^ sin((ot)] 


.3, x,-l 1 r _3, . 1 2 . 2 


UTo(s)U = -«■ [ To(s) {6 p +71 (cos((ot)-l) } -T:r(s) il i-. ^t^ .(cos (co t) -1 ) 

b J 2 y/r 5 


- T?'(s) (3-71) (cos (cot) -1) +1 1 ^( 3 ) sin(cot)] 

^x iu ^x 2J2 


UT 3 (a)U-^ 


^ [^ 3 ( 3 ) § is + 71 ^ (cos(cot) - 1 )} 


- Ji (a) (3 +7i) (cos (cot) - 1) 


+ T^(a) (Gos(cot) -1) +i T?(s) sin(cot)l 

^^ 2^15 ^ 2 J 2 


where' U = 


exp (-1 ^Qt) 


(1 + 7 iV 3 )^'^^ 


1 1 1 

In arriving at this equation we have made use of the Table 
II. 10, and, as usual, dropped the high frequency terms. The 
evolution of the tensor operators under the infl ;nce of 
are set out in Table II. 11. Having obtained the time evolution 
of the tensor operators under the influence of 

it is a s traighforward matter 'to obtain • the expression for the 
pulse responses. The single- and two- pulse responses obtained in 
this case are: 


II.D.2(i) Single-pulse response: 




4 kT 




sinjS cos(Aa) 

■ 2^1/2 ^ q 


(3+7? 


t ' ) 

(11.106) 


where /3 = 




(3+7?^) 


(11.107) 


Aw = (w -CO ) (1 +7? V 3 ) 
^ ^ . 

f = t-tj. 


Also <My> = <M^> = 0. 


II.D.2(ii) Tv^70-pulse response: 

(1 




CO. 


4 kT 


2 V 1/2 

7? /3) 


^ ^ ' 2^'l ' /2 cos^ (8/2 cos(AcOqt ) 


(3 +7?‘') 


+ sinjS cos8 cos(Aco (t-r)) 

\ q 

2 

- sinj 0 sin 3/2 cos(Aco (t-2T))) 


. . (11.108) 
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,n 


Table II. 11. The Evolution of the Tensor Operators T Under the 

^ ^ CJ 

Influence of the Hamiltonian for I = 3/2 


UT^U~^ 
/-j o • 


•4 (15 4 + (3+5 77^) (COS3-1) J 




-T^(s) T)(cos|3-l) + T^ — — (cos|3-l) 

3^I6 3yfl0 


+i Ti(s) ( 3 + 77 )? sin^ +i T?(s) sinjS 

5\/3 3v/5 


-1 To(s) f sinjs] 


v/3 


UT ?'( s ) u "^ 




[T^(s) {15 j3 + 2 ( 6+77 -317) (cos3-1) } 


+ T?(s) — ( 3 + 77 - 277 ^) (cos3-l) 

" 3JT5 


- T 3 (s) (COS 18 -I) +i T 


1 .1 


5^/3 


(3+i7)P sin^ 


-i T^(s) ~ f sin0 + i T^ ( 2 - 17 ) f sin3] 


v/l5 


UT^ (a)u“^ 




I [Ti(a) /+ 3(1-77^) (COS3-1)} 


T?(a) ^^4 ^(cos 3-1) -T^(a) —^( 3 + 277 ) ( l-rj) (cos/3-1) 

^ i. 3^ I 


- iTo(a) 4 f (i_yj) sin3] 


2 -1 
UT U 
r^O 


= ^ +(cos 3 -l)J +T^{s) ^(cosj3-l) +i ^ 1 ^^)'^’ sii^/3] 


UtJ(s)u”^ 


— rT?(s)f cosiB + iT^(s) 4 sinB] 

O ^ rs/ i. ^ /s; O 


f 


UT2(s)u“^ 


i [T 2 (s) I {'3 /+ n^(cos/3-l)} +T^ 4 (cos/3-1) 


/3 


+ iT^(s) 4^ sin/3] 


^^3 


. - contd 
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Table II 


3 -1 
UT U 

r>) O 


UTt (s)U 

^ X 


UT? (a)U 


UT?(a)U 


UT2 ( S ) u' 


.11 (contd.) 


P 


\ [Tq {15 (577^ + 12) (cosP-1) } +T 2 (s) — (cos^-1) 

V 15 


4 . yl — 5 — (cosjS-l) + i T^(s) (4 + 3 77)f sirijS 


5n/To 


+ i tJCs) — ^ (2-T7)P Sin3- i t:^(s), sin0] 


5x/30 


r<J O 


r) 


V30 


[t^(s) { 10 j>^+ (7 +2 77^ + 477 ) (cosjB-l) } 

+ T?(s) — - — ( 3-277) (cosj3-l) 

2 slls 

^ (3 +'77- 27}^) (cosj3-l)] +i T^Cs) -^^Psin^ 


+ T^Cs) [ 

. J2 


lQ\fl5 


vTo 


+ ITq (4 +3?7)f sin0 + i T^ ^oTs 


jjV tTi(a) 3 Y {30 +(3 +2r7)^(cos3 -1) 

(3 + 277 ) (cos^ - 1) 


+ T,(a) 

2 Pr 5 


- T^a) — ^ (3 + 277 ) (I- 77 ) (COS0-1) 

^ J- 


10 \.^ 


+ iT?(a) ( 3 + 277 ) sln3 ] 

^ J30 


= i Ft? (a) P cos3 +i T^a) sin3 + 1 T^ (a) (3+2 77)sin0 

J2 ./^n 




/ 3 O 


- IT? (a) — ^ (l-77)sin/3] 

5/2 


P 


■| [T 2 (s) {f%(cos0-l)) -T^ (cos^-l) 

+ T^ (cosj3-l) + 1 T?(s) — sing 

Pis v/2 

+ lT?(s) f sing -i T?(s) — ~ sing ] 


^/^o 


S'/I 


. . .contd. 
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Table II. 11 ( contd . ) 


UT 3 (s)u"^ = ~ [^ 3 ( 3 ) I [e/ + (3 +2rj^) (cos/3-1) 




+ T^Cs) — ^ (3-2n) (cos^-1) -Ths) (cos^-l) 

2^5 10 


^ . ^3 , . P . ^ . „1 2 

+ iTo(s) — sini3 - 1 T„ 

^/2 10 ^ 
UT 3 (a)u"^ = [T^Ca)' | { 2 J>^ + (cos 13 - 1 )} 

+ T^(a) — - — ( 3 + 27}) (cosp 1) 

^ 2 CS 


.^3 _llf_ 

■° I0I3 


77fsin3- iT^ sin^] 


T^^Ca) ~ ( l-n) (cosj3-l) +i T^Ca) — sin/3 ] 


where U = exp (-i ) 




3 


1 w' 

i= o, [-\ [T^i(a) Tha) - 

^ 3 r vi 5 j j 

( 3 + 7 ?) 

1 r3f w 


f = (1 + 



where x t 


t > (x + 2t ) 

UJ 


II. D. 3 Calculation In the 'Diagonal Frame' 

Pulse response calculations can be considerably simplified 
if we work in a representation in which the quadrupolar Hamil- 
tonian is in diagonal form. It can be easily shown that (for 
I = 3 / 2 ) that quadrupolar Hamiltonian can be brought into the 
diagonal form by the transformation operator 


2 . 

V3 


= exp X T^Ca)) 


where x = tan h / Jd 


. . (11.109) 

.. (II. 110) 


The transformation of a few selected tensor operators 

under U, are set out in Table 11.12. The total Hamiltonian 
d 

in the diagonal frame can be written as 


^'^( t ) =^2 


.. (II.lll) 


(0 


where }«g ^ '"o ' ' 


and = 2 cos (m„t) [Tha) { 1 1 4/5 (cosx-l) i 

1 (4 ^ J. 


+ Tj(a) {— (cosx-1) + ^ sinxi. 


+ To(a) {— sinx- i(cosx-l))] 
r^3 ^ -i 


(11.113) 


'd' refers to the diagonal frame. 
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Table 11.12. Transformation of a few Selected Irreducible Tensor 
Operators in the Frame Defined by the Transformation 
Operator = exp (2/\f3 x T^Ca))/ for I = 3/2 


2 -1 
^d^o^d 


2 2 

T^ cosX - sinX To(s) 
^ o z 




To(s) cosX + T^ sinX 


U3Tj(a)u-^ 


?^(a) [l + 4/5 (cosX- 1 )] + T^(a) [-“ (cosX - 1) 


+ “ sinX] 


+ T3(a) [-^ sinX -1/3 (cosX- 1 )] 




Ta ( s ) [ 1 + 4 ( cosX - 1 ) ] + T? ( s ) [-per ( cosX - 1 ) 

j srnx] 

+ T^(s) r 4 (cosX -1) + ■— sinX'l 

L 3 \/3 


1 -1 
^d^o^d 


fo (cosx-l)] - xl (COSX- I) ] 


T2(s) sinx] 


2 - 1/2 

where cosX = (1 Hhr) /3) ; 




sinX 


Now the evolution of tensor operators under the influence of 
can be obtained from Table II. 2 by substituting cOq by co . In 

w C[ 

the interaction representation defined by the transformation 
operator 



.. (11.114) 


the total Hamiltonian given in Eqn. II. ill can be written as 


lid AUrr 2 '^d 

n/6 J- 


.. (11.115) 


where Aco = (co _(*) ) 

q q q 


d 

and 




CO 

1 ‘ 3f 




(11.116) 


In arriving at this equation the high frequency terms have been 
omitted as usual. The evolution of tensor operators under 
^1 can be obtained from Table II. 3 by substituting £ by 0 


« (3+n)t^ 

where |2> = ' .j •• (11.117) 

(3 + V ) ^ 

<3. ^(3. 

Using the evolution of tensor operators under and and 
proceeding as earlier, the single- and two-pulse responses 
obtained are exactly same as those given by Eqns. 11.106 and 


p a 
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11.108, respectively. From here onwards we present all the cal- 
culations in the 'diagonal frame' only, 

II. E PULSE RESPONSES IN THE PRESENCE OF A WEAK ZEEMAN FIELD 

FOR NUCLEAR SPIN I = 3/2 IN NON-J^IALLY SYMMETRIC FIELD GRADIENTS 

In this section we. present the calculation of pulse res- 
ponses in the presence of a weak Zeeman field applied along X-, 

Y-, and Z- principal axes of the efg tensor. In this analysis we 

assume that « Wg and also; neglect the effect of Zeeman field 
during the r.f. pulse assuming that the spins scarcely process 
during the period when r.f- is on. 

lI.E.l Zeeman Field Applied Along X-Axis of the efg 

The Hamiltonian corresponding to the Zeeman interaction 
for this case is given by ■ I 

^zx = •* 

■Using Tables II. 2 and 11.12 

- CO [/i {1+ i(l-U)}T^(a) - (1+1 (l-Tl)}T^(a) 

+ |{1 (l-n) jT^Ca) j .. (11.119) 

where ^ = (1 -fU 


( 11 . 120 ) 


The evolution of a few selected tensor operators under this 
Hamiltonian are tabulated in Table 11.13. The single- and two- 
pulse responses obtained for this case (for Ao = o) are; 


II.E.l(i) Single-pulse response: 


<M > = 7 -p= (1 - siniS cos ( ---p -- t ' ) 


where = 


i^+^) .... t • t' = t-t 
2v 1/2 ^ ^ ^uj 


( 3+r) ) 


.. ( 11 . 121 ) 


II.E.l(ii) Tvj’O-pulse response: 


<M> = (l+nV3)^^^ — i-ltal 


(3+7?^) 


2 s 1/2 


[sin^ cos 3/2 cos(a)Q — t) 

+ sin3 cos3 cos(co^ - -y— (t-t ) ) 

- sin3 sin^3 /2 cos (w^ -- - - y- ( t- 2t ) ) ] 

.. ( 11 . 122 ) 


where t t 


t ^ {%+ 2tJ 


II. E. 2 Zeeman Field Applied Along Y-Axis of efc 


The Zeeman Hamiltonian corresponding to this situation is 




(11.123) 
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Table' 11.13. The Evolution of a few Selected Irreducible Tensor 

^ c3 

Operators Under the Action of the Hamiltonian 

ZjA 

T?(s) cos(a> — t) + i To{s) sin ( co — ■ y— ■ t) 

2 ^ 

Tj^(a) cos(o>^t) + i T^Ca) sin(co^t) 

To(s) cos (w — ' ^ ■■^t) + i T?(s) sin (w i . il- . r . I l. ) -t) 

-2 ° f ° f 

To (a) cos (o.t) + i T?(a) sin (w t) 

^ O -L O 

m2 

where U = exp 

I -?!<=>> 

+ I^a)U tl - i (1-,))] 
j> = (1 +77V3)^'^^ 

“o=^«o ■ 
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^ d 

Now, 






(11.124) 


Using Tables II. 2 and 11.12 and dropping the high frequency terms, 
we get 


- -icoo [| (1+1 (1-7?)) T^s) - (1 +|(l+r7))T^(s) 

f \/l5J 

+ -I (“1 +j (1 +'n))i^is)] 

.. (11.125) 


The evolution of a few selected tensor operators under this 
Hamiltonian are given in Table 11.14. Pulse responses obtained 
in this case are: as under. 


II.E.2(i) Single-pulse response: 


0) 


^ M y = — ^ (1 tuVs)^^^ " ' ^2 ' ^i 72 ' cos(w 

"" 4 kT {3+ri)^^^ ° 


f 


.. (11.126) 


where t' = t-t 


m 


lI.E.2(ii) Tv/o-pulse response: 




4 kT 


^+ , ^ . [sin/3 cos^/2 cos (co ■t) 

+ sinB cos6 cos(co^ . (ll?] .., ^ (t-x) ) 

- sinjS sin‘^/3/2 cos (o)^ , ( ,. A?I3J , ( -t, 2 t ) ) ] 

J 

. . (11.127) 


where t t 


UJ 


t ^ (i:+2ty^) 
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Table 11.14. The Evolution of a few Selected Irreducible Tensor 
Operators Under the Action of the Hamiltonian 


UtJ ( s ) u “^ 


2 ’ 2 
T-, (s) cos (w^t) - T-Ca) sin (co^t) 
^ i . O o 


UTha ) u "^ 

± 


T^Ca) cos (o)^ -T 2 (s) sin (o)^ - Q-t i ?) . t) 


UT?(s)u"^ = t5(s) cos t) + T?(a) sin (o) t ) 

r>j ^ o y X o y 


UT ^( a ) u "^ 

fO Z 


2 2 
To(a) cos (w^t) + T. (s) sin (w t) 

^ O rJ X O 


2 -1 
UTqU 


where U = 


exp (-iX^^t) 


^ d 


-iw© [tJCs) 3(1+ j (1 +h)} 

+ T3(s) I { -1 + ^ (1 +T1)J] 


T?(s) {1+^ (1+^) J 




(1 +n^/3)^^^ 



123 


I I , E . 3 Zeeman Field Applied Along 2- Axis of efq 

The Zeeman Hamiltonian for this case is given by 


= 


rr .1 

O O 


(11.128) 


and 


= 'Jq °d^z '"d" 


- -w [r(i-)-4/f )tJ - — (~ - Dt;;] 

° ^ 3JlO f 


(11.129) 


The evolution of a few selected tensor operators under the 
influence of this Hamiltonian are set out’ in Table 11.15. Pulse 
responses obtained in this case are as follows. 

II.E.3(i) Single-pulse response: 


CO 


<M^> = 


4 kT „2?l/2 cos(co^t‘) 


2 , ^, 1/2 ( 3 + 7 ^) 


ii+rT) 


.. (11.130) 


v/here t' = t-t, 


UJ 


II.E.3(ii) Two- pulse response: 


<( M > = 


00 


n 9 1 /9 (34-7}) 2 y / V 

^ ^ [sin^ cos 0/2 cos(oo t) 


4 kT 




+ sing cosg cos ( 05_ (t-x) ) 


- sin|3 sin^g/2 cos(co^(t-2x))] 

. . (11.131) 


where x » t 


Ui 


t > (T^ + 2t^jy) 
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Table 11.15. 


The Evolution of a few Selected 

• n of 

Tensor Operators Under the Actao 


i-j (ij 

Hamiltonian 

2j 2 


UT?(s)u'”^ = T?(s) cos (to t) - iT?(a) sinCto^t) 

/v/ 


UT? (a)u"^ 

^ X. 


-^(a) cosCto^t) - iT?(s) sin((o_t) 


UT?(s)u“^ 


UT?(a)u"^ 

<*01 ^ 


2 t) 


= To(s) cos(a)_ — t) - iTp(a) sin(a)^ 

^ z ^ P ^ 

= To (a) cos (to — t) - iT^(s) sin(w^ p 
^ Z op rNy J 


where U = expC-i^^ig 


‘^o ['^o I (1 + ^ 


3 4 


yio 


j> = (1 +T]V3)^2'2^ ^ 
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Before concluding this section, it is worth mentioning here that 
since the calculation of the evolution of complete set of tensor 
operators (for spin I =3/2) under the action of Zeeman Hamilto- 
nian in different situations, becomes cumbersome we have present- 
ed here the results of the evolution of only a selected set of 
operators . 


II. F THE EVOLUTION OF TENSOR OPERATORS UMDER' THE ACTION OF r.f . 

HAMILTONIAN; ROLE OF PHASE OF r.f. 

So far we have considered the r.f. pulses, in a two-pulse 
sequence v^ith the same phase. However, if we wish to study the 
response of a nuclear spin system to an r.f. pulse as a function 
of phase then the r.f. Hamiltonian should be modified as follows: 

( cos ■ (w t + /). TjCa) .. (11.132) 


where / is the phase of the r.f. 

We consider the case of spin I = 3/2. In the interaction 
representation it is given by 




-1 „-l 


u. 


. (lIol33) 


(3+r7) 


[(•! T^(a) + T;?(a))cos0' 

v/3(3+ti2)1/2 1 >- 5 ,,1 ^TS 

- i sin/ — T^ (s) 1 
sf2 


.. (11.134) 
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For 0" = 0 this 
For = 90°, 


reduces to Eqn. II 

(3-W) 
n/ 3 ( 3+n^) 


.116. 


_ 1 _ 

v/2 


T^(s) 

rs> J. 


(11.135) 


In arriving at the Eqn. 11.134 we have made use of the Tables 

II. 2 and 11.12 and dropped the high frequency terms. The evolu- 

ol 

tion of tensor operators under the influence of ( ^ is set 

1 0 

out in Table 11.16. Tables 11.16 and II. 2 can be used to cal- 
culate the responses of spin 3/2 nuclei to sequence of pulses 
involving different relative phases. 

It should be quite evident now from the above presenta- 
tion that the tensor operator formalism can be used with rela- 
tive ease to describe the pulse responses in pure NQR. 

SUMMARY 

Tensor operator formalism has been applied to pulsed NQR 
spectroscopy for the first time. The versatility of the forma- 
lism has been illustrated by calculating the single- and two- 
pulse responses in various situations in the case of spin 
1=1, 3/2 and 5/2 nuclei in single crystals. The results 
obtained here are in full agreement with those of others, wherever 
available in the literature.. The evolution of tensor operators 
under the action of the r.f. Hamiltonian including arbitrary 
phase is presented for the spin 3/2 case. For the sahe of 
illustration typical recordings of pulse responses are also 
included. 



'able 11.16. The Evolution of the Tensor Operators T^(a/S) Under 
the Influence of the Hamiltonian ( 'K?) ^ 


1 -1 
IT U 


rT^a) U“^ 


[1 t I (cos3-l)]t (cos^- 1 ) 

+ sir^ [i{ ^ tJ(s) + ■— T^(s) } COS0 + ^ T^^a) sin^] 

T];(a} [1 + I (1 + 3 sinV)(cos3-l)] - T^(a) cos2iZf.(cos0--l> 

JTk 


5 tJ-Cs) + - 2 - T?(sM nnsffi + -i T?(a) sin^] 


jtJ^(s)u“^ 

rJ 


“1 O 

+ i — T^Cs) sin22f(cosj3-l) - ^ T?(a) (cos/3-1) 

2^/2 ■ 3 ^3 

- i sin/3 [-1 { ~ T^(s) +. — tJ 3 siniZf + T^(a) cos/Zf] 

v/i ^/2 , 3-2 

t];(s) [1 +-| (3cos^/2( + 1) (COS0-1)] +-L. T^(s) cos20{cos^- 

^ n/Ts 


- 1 T- 


T^(a) sin20' (cos^-l) T^ (s) (cos/S-l) 


;in0] 


2 -1 
JT_U 


JT^ Ca)u'^ 


JT?(s)U"^ 

^ X 


JT?(a)u'”^ 

/%» ^ 


Jt 5 (s)U-^ 

/vJ 


rslO 


. =inB [ic 1| - 4 T^=, . g 

COSJ3 +3ine [ i T?(s) coBSS +— { xha) + "4* 

-o ^ ■ 5 -1 ^5 , 

T?(a)[l +sin^i2f(cos3-l)] +i[-X_ tJ'(s) +-2- T? (s) ]sin2lZ^(‘^°®^ 

-[3 ^o + 2 o 3 

Ti(s)[l+ cos^0(cos^-l)] — i- Tha) +-1 tJ ( a) 


T?(s)[l+ cos^0(cos^-l)] r 

2/6 ^ 


+ i T^ cos^ sinjs 


T^Ca) [ 1 +sin^0'(cos/3-l)] +i -i To(s) sin 2 j 2 ((cos^-l) 

+ sing [+ 4 ;i‘=> - ■§ 3 'i‘=) S 3 <®d 

O O 'I o 

t:;(s) [ 1 +sin 0 '(cos 0-1)1 +i -i T^(a) sin2i2((cos3-l) 
+ sin0 T];(a) - Tj(a) - T^(a)] sin/Zf 

T^ [1+1 (cos3-l)]+ ^T^ (COS0-1) 


+ sine[i{ 22 ^ Ths) + 2 ^ Tf(s)} cosiZl +-^ sin/Zf Ti(a)] 
5 V.X 5 ^ ^ " 


contd * 


Table 11,16 (contd.) 


UTi (a)u“^ 

Jfsj -*• 


T^(a) [ 1 +T^ (3+4 sin^0) (cos/3 -1)1 + (cos|3-l)[-^§- T^(a) 
J-u j 


"^9 1 1 2 

■ T. (a) COS20' + i sin20' Ti(s)] + i sinjB 

2v/ro 


[i{-~ Tp(s)} singf + T^(a) cos0'] 


UT? (s)u”^ 

4^ X 


T^(s) — cos20(cos3-l) +Tt(s)[ 1+ — (3+4 cos^0’) (cos3-l) 

IG'/IS -^1 10 


+ T^(s) (cosj3-l) - i — ^ — T?(a) sin20' (cos 3 - 1) 

■ 2n/ 5 2-r30 

, r.w 3 „1 , 2\/2 „3 . n/3 3 


+ Sin3 [i{-^ t;^ + ^ t;: + t:;(s)} cos0 — T^(s) sln^f 

PL 5^ .o 5 ^o -2 ^ 2Vro 


UT^(a)u"^ 


O / X 

T2 ( a) [1 + cos 0' (cos3 -l)] -i sin20' (cos^-l) 


+ isinig [{ — (a) + T,^(a) + ~ T^(a) } cos0] 

542 


UT^(s)u‘'^ 


T2 (s) [ 1 +cos^0' (cos3-l)] -i ^2^^^ sin20’(cos3-l) 


+ isin/3 [{ ^ t|(s) + T?(s) + — T^(s)} cosi2f] 

5^2 "" JTo 42 


UT?(a)u“^ 

rj O 


Tj (a) [ 1 + "I (cos3-l) ] + [-^ “ ij' -Jl ^ (cos3-l) 

_ 245 ^ 

+ sln3 [— ^ sin0' T?(s) + ~ T?(a) cosj2(] 

2'/2 \i2 


3 / V — 1 
UT.(s)U 

r-J-J 


T2 (s) [ 1 + j (cos3-l) ] + T^(s) - T^(s) ] (cos3-l) 

2 >^5 


+ sin3 f — ■ " ■ T?(a) sin0' + — T~(s) cos0'l 

2V2 42 -'2 


where ( %\) 


(0. — [(•! Tl’ (a)) COS0’ - i — ^ T? (s)sir^ 


exp (- i ( ^ 1^0 t ) 


(3+77) 

(3+ri)^'^ 


t 



129 


REFERENCES 

[1] R.R. Ernst, G. Eodenhausen and A. Wokaun, "Principles of 
Nuclear Magnetic Resonance in one and two Dimensions", 
Clarendon Press, Oxford (1987). 

[ 2 ] a) S. Vega and A. Pines, J. Chem. Phys . , 66, 5624 (1977). 
b) S. Vega, J. Chem. Phys., 63, 3769 (1975) . 

[ 3 ] R.S. Cantor 'and J.S. Waugh, J. Chem. Phys., 73, 1054 (1980). 

[ 4 ] a) D. Ya. Csokin, Phys. Status Solidi(b), 102, 681 (1980). 

b) D. Ya. Osokin, J. Molec. Struct., 83, 243 (1982). 

c) D. Ya. Osokin, Phys. Status Solidi(b), 1^, K7 (1982). 

[ 5 ] M. Matti Maricq, Phys. Rev., E^, 4501 (1986). 

[6] S. Vega, J. Chem. Phys., 6B, 5518 (1978). 

[ 7 ] " A. Wokaun and R.R. Ernst, J. Chem. Phys., 67, 1752 (1977). 
[s] G.J. Bowden and W.D. Hutchison, J. Magn. Reson., 61 _, 403 

(1986) . 

[ 9 ] G.J. Bowden, W.D. Hutchison and J. Khachan, J. Magn. Reson., 
67, 415 (1986). 

[ 10 ] G.J. Bowden and W.D. Hutchison, J. Magn. Reson., TCi, 361 
(1986). 

[llj G.J. Bowden and W.D. Hutchison, J. Magn. Reson., 7_1, 61 
(1987). 

[ 12 ] a) A.R. Edmonds, "Angular Momentum in Quantum Mechanics'', • 
Princeton Univ. Press, Princeton, N.J. (1957). 

b) M.E. Rose, "Elementary Theory of Angular Momentum", 

John Wiley, New York (1967) . 

c) H.A. Buckmaster, R. Chatterjee and Y.H. Shing, Phys. 
Status Solidi, 13, 9 (1972) . 



130 


[is] W.D. Hutchison, Fh.D. Thesis, University of New South 
Wales, New South Wales, Australia (1987) . 

[ 14 ] A.K. Saha and T.P. Das, "Theory and Applications of Nuclear 
Induction", Saha Institute of Nuclear Physics, Calcutta 
(1957). 

[ 15 ] T.P. Das and E.L. Hahn, "Solid State Physics", Suppl. I, 

Academic Press (1958). • 

M. Bloom, E.L. Hahn and B. Herzog, Phys . Rev., 97_, 1699 
(1955). 


[ 16 ] 



CtlAFXER III 


I GAT ION 
FOI^MAL I SM FOR 
FERTURBED SR IN 


OR TENSOR ORERATOR 

TRE STUDY OR ZEEMAN 

ECHO ENVELORE MODULA 


TIONS IN SRIN I 


3/2 NUCLEI 



132 


This chapter deals with the theoretical investigation of 
Zeeman-perturbed spin echo envelope modulations (ZSEEM) in spin 
I = 3/2 nuclei in polycrystalline specimens . It consists of 
two sections. Section III .A presents , a brief discussion on 
spin echo envelope modulations while section III.B outlines the 
details of the calculation of ZSEEM function for the case of 
spin I = 3/2 using tensor operator formalism. 


III .A ZEEMAN-PERTURBED SPIN ECHO ENVELOPE MODULATIONS 

The parameters that characterize the NQR spectrum are the 

2 

quadrupole coupling constant (e qQ) and the asymmetry parameter 
U of the efg. As has been pointed out in Chapter 1/ in the case 
of the nuclei with spin I = 3/2 in the absence of a Zeeman field 
there is only one transition with a frequency 


0)^ = (1 + UV3) 


1/2 


. .. (III.I) 


Between the pairs of doubly degenerate energy levels. Hence, it 

is not possible to obtain both the quadrupole interaction para- 
2 

meters e qQ and “n by measuring co alone. It is necessary to 
apply a weak Zeeman field to lift the degeneracies of thei energy 
levels and then the resulting four transition frequencies (see 
Pig. III.I) can be used for the determination of both these 
parameters [1,2]. 
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Zeeman NQR spectroscopy has been used to obtain these 
interaction parameters in single crystals [3] as well as in 
powder samples [4-6}. 

The ZSEEM method is ideally suited for obtaining Q using 
pulsed NQR spectroscopy. One can employ both two-pulse and sti- 
mulated echo sequences for this study. However, we confine our- 
seleves, in this chapter, to two-pulse SEEM. Spin echo envelope 
modulation techniques are well known in the area of ESR spectro- 
scopy [7]. The energy level diagram of a Zeeman-perturbed qua- 
drupolar nucleus with spin I = 3/2 is given in Fig. III.l. A 
coherent pulse excitation of all branching transitions would 
lead to modulations of spin echo signals and the echo amplitude 
t{2x) plotted as a function oft, the interval between the two 
pulses will exhibit modulations (see Fig. III. 2). The frequen- 
cies of these modulations correspond to the Zeeman splittings in 
steady-state experiments. The echo envelope modulations can be 
seen even when the modulations on the FID and spin echo signals 
are damped down due to static inhomogeneities. Since the 
energy levels of a Zeeman- perturbed quadrupolar nucleus are 
r)-dependent the Zeeman perturbed spin echo envelope modulations 
(ZSEEM) function £(2 t) is expected to contain information on rj . 
£(2'T) also depends on © and 0 which define the mutual orientation 
of the external fields and the principal axes of the efg tensor, 
owing to the dependence of the various energy levels and the 
transition probabilities on these angles. 
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In order to observe the Zeeman-perturbed spin echo envelope 
Ddulations the following conditions are to be satisfied. i) The 
seman splittings, and (see Pig. III.l) are larger than the 

requency width of the spin packets that make up the inhomogene- 
usly broadened NQR line, ii) The r.f. field should satisfy the 
ondition and "I^he first condition ensures the occur- 

ence of several modulation periods within the echo decay time while 
he latter ensures the excitation of all the branching transitions . 

There have been several reports in the literature regard- 
ng the ZSEEM patterns for spin I = 3/2 nuclei in polycrystal- 
ine samples [8-12], Ramachandran and Narasimhan [ll] have 
iven a more general analysis of ZSEEM patterns. In their 
pproach, they have followed the evolution of the density matrix 
•f the spin system by matrix transformation of the density matrix 
‘ollowing r.f. pulses and during evolution periods and obtained 
in expression for ZSEEM function, £ (2 t). They have computer 
limulated the average value of this function, £(2t), over 
:omplete (©, 0) space, for various values of and T) and comp- 
ired the simulated ZSEEM patterns with those obtained experi- 
mentally for polycrystalline samples. Their study shows that 
E(2t) contains the frequency components co^, co^, and 

[Wa- (see Fig. III.l) contrary to the earlier reports in 

the literature. The sensitivity of £(2 t) to T) is clearly ref- 
lected in the ZSEEM patterns obtained by them. 
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Recently, Narsimha Reddy and Narasimhan [13] have imple- 
mented the extended time excitation technique due to Schweiger 
et al. [14] in pure NQR to obtain ZSEEM pattern in a single short 
experiment. 

Before presenting the theoretical analysis, we illustrate 

here, with some experimental recordings, the nature of ZSEEM 

patterns. In' a two-pulse ZSEEM experiment, the echo amplitude 

is sampled as a function of t , the pulse separation, to obtain 

the echo envelope. Since in the presence of a weak Zeeman field 

the echo amplitude gets modulated the echo envelope also' shows 

modulations (see Fig. II. 2). Figures 111.3(a) and III, 3(b) show 

35 

the experimental ZSEEM patterns obtained from Cl in polycrys- 
talline sample of AgClO^- From these recordings it can be seen 
that the modulation period decreases as the strength of the magne- 
tic field increases. The overall decay of echo envelope modulation 
is due to the spin-spin relaxation processes. If this experiment 
has to be done by manually varying x and sampling the echo ampli- 
tude, mostly by a boxcar, it would take more than an hour. The ZSEEM 
patterns shown in Figs. 111.3(a) & III, 3(b) were obtained using a 
microprocessor-controlled pulsed N.QR spectrometer described earlier 
(see Sec. I.C.3). In the ZSEEM experiment, with this spectrometer, 
the generation of two pulse sequence and incrementation of t was 
done by microprocessor under software control. For each x value 
the complete echo has been acquired into the signal analyzer. 

The echo maximum points have been chosen, off-line, with the 
help of a BASIC program to obtain ZSEEM pattern. Thus, with 
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a magnetic field of 2 9-75 Gauss. 
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this spectrometer a complete ZSEEM pattern can be obtained in 
about 10 minutes. 

The application of tensor operator formalism [15-2o] for 
obtaining the two-pulse ZSEEM pattern for spin I = 3/2 nuclei in 
polycrystalline samples will be discussed in the next section. 

III.B THEORETICAL INVESTIGATION OF ZSEEM USING TENSOR OPERATOR 
FORMALISM 

III.B.l Eigen Values of Spin I = 3/2 Nuclei with Non-axlal 
Field Gradients / in the Presence of a Weak Zeeman 
Field 

The net Hamiltonian in the laboratory frame for a spin 
I = 3/2 nucleus in a non-axially symmetric electric field, grad- 
ient and in the presence of a Zeeman field can be written as 

2 

where Kq = jfxfirTT [3 •• 

is the quadrupolar Hamiltonian and 
- -THh.I 

O 

= sine COS0' + I,, sin© sinfZ^ + I„ cos©] 

. . (III. 4) 



is the Hamiltonian corresponding to the interaction of the 
nucleus with the Zeeman field which is applied in a direction 
with polar coordinates (©^ 20 in the principal axes system of 
efg. I is the nuclear spin, -y is the gyromagnetic ratio of the 
nucleus and is the applied magnetic field strength. If is 
zero, there are two pairs of doubly generate energy levels for 
I = 3/2, which are given by 

2 

El, 2 = f .. (III. 5) 

. ®3,4 = ^ ** 

where f = (1 +ri^ .. (III. 7) 

The transition frequency between these energy levels is 


V 


Q 


= y 


Q 


2 ^ 
6 .30 


s S 


. . (III. 8) 


In the presence of the magnetic field the energy levels are 
split depending upon the strength and the direction of rela- 
tive to the principal axes of the crystalline efg. For arbitrary 
T) values, assuming the magnetic field to be weak, i.e., « e qQ, 

the energy levels can be obtained by treating the Zeeman term as 
a perturbation on the quadrupole levels. For arbitrary values of 
T) the energy levels are given by [l»3] 




The energy levels diagram of a spin I = 3/2 nucleus, in the 
presence of a weak Zeeman field -is shown in Fig. III.l. 

In the following sub-section we present the explicit 
expressions for the interaction Hamiltonians, required for the 
discussion of 2SEEM. 


Relevant Hamiltonians 
III.B.2(i) r.f. Hamiltonian: 

We are interested in the calculation of the pulse respons- 
es of spin I = 3/2 nuclear spin system for the case |[ and 
for an arbitrary orientation of these fields relative to the 
efg principal axes system. The r.f, Hamiltonian for this case 
is 

'^ 2 _(t) = -2w^ii cos(Wg.t) [r^ sin© cos/ + sin© sin/ 

+ cos©] . . (111.12) 
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the resonance frequency in the absence of the Zeeman field. In 
the above eqn. G and 0 are the polar and azimuthal angles of the 
r.f. field vector in the principal axes system and 

“q=“Qf- 

where Y is the gyromagnetic ratio of the nucleus. 

is the amplitude of the r.f. field 
2 

and = e qQ/2 . 

IILB.2(ii) Total Hamiltonian in terms of tensor operators: 

The total -Hamiltonian in the presence of an r.f. field 
is given by , 

■Kct) =Kq +K^(t) .. (III. 13) 

where & K]_(t) are given respectively by Eqns . III. 3, 

III. 4 and III. 12. In terms of tensor operators (see Table II. 1) 
these operators can be written as 

^ tJo ^ .. (III. 14) 

■^5‘o ^3 ^ •• 

and = -2c*)^1i [cos (m^t) ][-T^(a) + T^^ (s) P2'*'’^o ^3^ 

.. (111.16) 

where co = H 
o o 

= sin© COS0, P 2 = i sin© sinjzf/ = cos© .. (III. 17) 
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Total Hamiltonian in the 'Diagonal frame': 

As pointed out in Chapter II, the pulse response calcula- 
tions can be considerably simplified by working in a representa- 
tion in which the quadrupolar Hamiltonian is diagonal. In all 
the subsequent calculations in this chapter we set ti = 1, and 
express energies in units of radians per sec. The transforma- 
tion into the 'diagonal frame' can be effected by the operator 
defined by 

U, = exp (— X T^(a)) .. (III. 18) 

a 

Making use of Table 11.12 (see Chapter II) the total Hamiltonian 
in the 'diagonal frame' can be written as 

^^(t) = = Rq + 'R^(t) 

where Rq = -^ T^ 

= -% [^1 

and R.^(t) = -2w^ cos (Wgt) [j] 

with [j] = -P], { T^(a) + a .2 T^(^) + ^^(a)} 

+P2 ( t];(s) + a^ T^(s) + ag ^3(3)} 

■^^3 < lo * fo * = •• 


(111. 19) 

(111. 20) 

(111. 21) 

(111. 22) 
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and 


^4 = ^7 


(l±± 


) 


a 2 = -V- [3(l-f) +5r)] 

snsf 

■= ij (f+’l-l) 

^6 = ^ -f+'l] 


a 


8 


a 


9 


_4 

3jIof 
2'J2 n 
3f3 f 


[1-f] 


(III. 24) 


Since we will be doing all the spin dynamics calculations by 
following the density matrix of the spin system in the quadru- 
polar interaction representation (QIR) we now transform the total 

J 

Hamiltonian in the 'diagonal f rame ' , ^ , into the QIR. 


III.B.2(iv) Total Hamiltonian, in the quadrupolar 

interaction representation: 

(3. 

The total Hamiltonian in the 'diagonal frame , can be 
transformed into the quadrupolar interaction representation by 
the transformation operator 

Ug = exp (i t ) ... ( 111 . 25 ) 

and is given by 

K'^(t) = Ug^'^(t)u"^ = +ill(b) 


. . (III. 26) 



where 


= -2co. cos(w t) U [j]u"^ 

Si Si Si. 

- “i 3J [Pi‘<3+'»[|3i<®> + 

-P 2 {( 3 -n)[| j^Cs) + ^ tJ(s)]3 

+P 3 C+ ^17 [ 52 ( 5 )])] .. (III. 27) 

and = -o>^ 

-“of ['’ill [f ■'(l-'l)] jila) - ^ [f+ (1-7)] T5(a) 

+ I [j>-(i-r])] T 3 (a) } 


-P 2 i-| [i’+d+r?)] T^(s) - 


1 

v/Ts 


[y+(i-j})] tJ(s) 


+ I + (1+T7)] T3^s) } 

-Poll (y+4) (-1 +f) T^}] .. (III. 28) 

3 5 kO 3 ^q -o 


In arriving at Eqns , III. 27 and III. 28 we have made use of 

Table III. 1 which gives the evolution of tensor operators under 

d. ' 

the influence of the guadrupolar interaction As usual, 

high frequency terms have been dropped since they are non-secular 
terms. 

In order to calculate the pulse responses by studying the 
time evolution of the density matrix in QIR following the r.f. 
pulses and during the evolution under the Zeeman interaction we 
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Table III.l. The Evolution, of the Tensor Operators Tg(s,a) 

d!, 00 cr 2 

Under a Quadrupolar Interaction T 

• Q ^ o 

for I = 3/2 


- i TI + tJ(s) [cos(<o^t)-l] 




iTha) sinCu t):] +T?(s) cos Cw t) 

‘-D Q 


IT? (a) sin(w t) 

JTS ^ 


U T?(s)U^^ = T?(s) cos(w t) -i T?(a) — | sin(co^t) 

q^2 q ^2 q w2 J3 q 


U,t2(s)u/ 


5 V 5 [cos((S t)-l] 


-iT?(a, 4 


i (a) -— — sin(w t) + T. (s) [-r cos(w t) +^] 

1 Jio q -1 L5 q 5J 


U^Tj(s)Ug 


0 1 
- iTo(a) -TT sin(a) t) + t:;(s) cos(w t) 
w z " z q z q 


where = exp(-i Xq t) 

For the antisymmetric combinations T^(a,s) replace s 

q o o 

© O Q T) ■ 1/^2 

by a and a by s iri. the above equations. =-y-^(l +2_) ^ ■ ; 
1 2 3 3 

T , T , To(s) and T^Ca) all commute with the Hamiltonian, 

/V O ^ J o 



need to know the time evolution of the tensor operators under 
the r.f. Hamiltonian and ^ respectively. In what follows 

-L Z 

we present the calculation of the time evolution of tensor 
operators under these interactions. 


III.B.3 The Evolution of Tensor Operators Under the Influence 

^ d. 

of the Hamiltonian 

In order to evaluate the time evolution of tensor operators 
under the influence of any interaction characterized by an Hamil- 
tonian, first of all one has to establish the commutation rela- 
tionships of the tensor operators with that Hamiltonian. One 
may use either the technique of "nested commutation relation- 
ships" [ 19 ] or the "harmonics of the motion method" to arrive at 
the evolution of tensor operators under that Hamiltonian. The 
commutation relationships between the Hamiltonian and the 

tensor operators can be grouped into two distinct sets. It is 
noted that these two sets do not mix, in the sense that commuta- 
tion of any member of a given set with the Hamiltonian pro- 

duces only members of that set. The set which is relevant for 
the present calculations is given in Table III. 2. The time 
evolution of a few selated tensor operators under the action of 

rs.» d 

calculated using the ‘nested commutation relationships' is 
given in Table III. 3. 
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Table III. 2. The Commutation Relationships of a few Selected 
Tensor Operators with the Hamiltonian for 

I = 3/2 


[ = 


“1 [-^1 + Pj T2(a)(3-n) 


- J’3 2”] 


[H^, Ti(s)] 


r_p J n 


[*?' =“l ih [-^2 Io<3-7,)] 




[-P, 2^1] 


[^1' = = 0 


where = 


1 CJ’i(3+l) + ■;|r ijca)] 

-P2(3-„)[ |t1(s) + Jr tJ(s)] 


= cos© sin0 
= i sin© sin0' 


= cos© 






Table III. 3. The Evolution of a few Selected Tensor Operators 

rv/ ^ 'Ar 

Under the Influence of the Hamiltonian 


2 -1 
UT U 
~ o 


= T^ cos£ + i sin£ [c^ ^2 ^3 


UT^(s)u"^ 


tJ(s) t (cos£-l) [gJ tJ(s) + C^C2 


+ C2C3 T2(a)]+ IC^ t'^ sing 


UTj(a)u“^ = Tj(a) + (cos£-l) [-C^C2 ^^(s) - T^Ca) 


- C2C3 T2(a) ] -iC^ sing 


UT2(a)u“^ = T^Ca) + (cos£-l) [c^C^ T^(s) + C2C3 T^(a) 


+ c| T^Ca)] + iC^ T^ sing 


UT2(s)U~^ 


T^Cs) 


where U = exp(-i t^j^^ ) 


C . = 


P^(3+t)) -P2(3-rj) 

' ^2 


2t)P, 


^/ 3 f 


X 




X 


, C3 




X = 


[ — --y [sin^G {(9+ rj^) -6t) cos20'J + 4ti^ cos^©]]^'^^ 
(3+t) ) , 


£ - X 


is defined in Table III. 2. 
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III.B.4 The Evolution of Tensor Operators Under the Influence 

" '■ ■■■'"' fU " ' ' ' ' in-iiii- -T' 1 ir -1,. , , ■ 

of the Hamiltonian ’ 

z 

In this. case also the commutation relationships between 
tensor operators and can be grouped into two distinct sets. 

2j 

The commutation relationships of the tensor operators in one of 
the two sets with are presented in Table III. 4. From this 

table it is clear that the operators in this set cannot be 
divided^ any further into sub-groups such that commutation of 'Kr 
with any linear combination of the first sub-group produces a 
linear combination of the second sub-group and vice-versa. Hence, 
it is not possible to evaluate the time evolution of these opera- 
tors under the action of using the 'harmonics of the motion 

method' . Also, it is rather difficult to get closed-form express- 
ions for the time evolution of these operators using the technique of 
'nested commutation relationships'. However, if : we set J3'= 0, then from 
Table III. 4 it is evident that the commutation relationships can 
be divided into two sub-groups which enable one to use the 'harmo- 
nics of the motion method' to evaluate the time evolution of the 
tensor operators in this particular case. The time evolution of 
this set of tensor operators in this particular case are given I 

in Table III. 5. | 

Since neither of the two methods discussed above could be - ^ 

used for the calculation of time evolution of the tensor opera- | 

fV d "I' 

tors under the action of ^ (see Eqn. III. 28) we have resorted | 
to the matrix multiplication approach for this purpose. 


First 
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Table III. 4. The Commutation Relationships of a few Selected 

fN/ J 

Tensor Operators with the Hamiltonian 

= ..o [-h ^ 

= “o [-h 3'^®> - ^3 3i<=d 

= “o [-^1 = > - ^2 - ^3 j ] 

where cj = 7 H 
o o 

="^0 p T^(a) - ^ [f + (1-79)] T^Ca) 

+ § [f - (1-^)] T3 (a)) 

- ^2 [f+(l-7i)] T^(s) 

^15 

+ I [-f + (1+T))] T^(s)) 

- P 3 {| (f+4) tJ - (-l+f) T^J] 

P^, P 2 , P 3 and ere defined in Table III. 2 
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Table III. 5. The Evolution of a few Selected Tensor Operators 

rj j 

Under the Influence of for I .= 3/2 


T^CsIu”^ = T^ 


^2 

l(s)[- J- oos(Rt) + ^ 


^ cos (St)] 
3 


- T 2 (a) - (cos(Rt) - cos (St))] 

+ i{ T^ ( a) sin(Rt) - sin(St)] 

sin(Rt) - -|^ sin(St)]j 
^ 3 ■ 3 

rr.2, ^r^l^2 

(a)U = T. (a)[— s — cos(Rt) - — T — 

- ^ ^ -L 1 0 X. 0 


-:= — cos (St)] 
^3 


+ T 2 (s)[-|~ (cos(Rt) - 

2 r^^2 

- i{Tf(s)[— ^ sin(Rt) - sin(St) 

^3 ^3 

2 f2h^ f h, 

+ T 2 (a)[— = — sin(Rt) ^ 


— sin(St)]} 
^3 


UzT2(a)u”^ = T2 (s) ^ [cos(Rt) -cos (St)] 


+ T2(a)[^ cos(Rt) - ^ cos (St)] 

3 3 

0 0 

- i{T?(a)[r;i sin(Rt) - — sin(St)] 

' f f , 

+ T2(s) [~ sin(Rt) - ^ sin(St)]j 
"3 ^3 


.'.contd 
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u T=(s)u:^ 


0 0 f 0 

T^(a) ' [cos(Rt.) -cos (St)]- ^ cos(Rt) 

3 3 


ge. 


cos (St) J -i{T-(s) [— — sin(Rt) - 


e^g 


o ■ 3 3 


sin(St)] 


+ T^(a) [— ^— - sin(Rt) - — sin.(St)]} 


2 -1 
U T U 
Z^o z 


= T 


In this table the follov/ing definitions hold: 

= exp(- 13 (^t) 

= “o J [ ‘I ^ [f + d-'rt] T^Ca) 

+ § Cf-d-uHljCa)! 


- P 3 l| [i“ + 4 ] tJ- ^ 1 ] 


= sinG, 


cosG, a^ 




£ 

f 


Ki^ 2 = I [{Pid+a-^) +P3(l+cJ)} + /"{ [?£ ( l+aJ)+P3 ( l+cj) ]^- 

4 (p 5 a^-p|c^)^}] 

Ti = Jk^, ^2 ~ ^2' ^ ~ ^1*^0' ® “ ^2%' “o 

2 2 2 

g = ^ 1^3 2 ~ [^3 2 ^ 

® 1,2 ~ C ^ 3 ^ ■*■ ^ 1 ^ 1 , 2 ^"^ 1,2 

^ 1,2 ~ "”[^ 1 ^ 1 ^ , 2 , 2 ^ ^3 * 


^1^2 "®2^1 
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we have converted the into its equivalent matrix form (see 
Table III. 6 for the matrix representation of tensor operators), 
which is given .by 


a 0 0 b 

0 c d . 0 

0 d* -c 0 
b* 0 0 -.a 

where a = -1/2 ^ 

b = ^ [-Pl{f-(l-77)} + P2{ 

b* = ^ [-P3_ii>-(1-Ti) 5 -P2{(.1 +t]) -$]] 

d = -^ [-P^{ft(l-77) } +P2{ f + (1+^1 ) I ] 
d* = ~ [-PiCf +(1+7 t)} - P2{y + (l+r7)}] 

Now, U - = exp (-i t.) 

= exp (-i[B]coQt) 

= cos ([sjo^t) -i sin([B]w^t) 

Since we find 





(III. 29) 


(111. 30) 

(111. 31) 


. . .contd. 
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Table III 06 . Matrix Representation of Tensor Operators for Spin I = 3/2 




l|> 

!■?) 

l§> 

ll> 










T^Ca) 

^/3 

= - 

0 

1 

0 

0 

tJ(s) =- ^ 

0 

-1 

0 

n 

0 

ml 3 

T = o- 
0 2 

1 

0 

0 

0 



1 

0 

2/\/3 

0 


0 

0 

-2/^/3 

0 


0 

1/3 

0 

0 



0 

2//3 

0 

1 


0 

2//3 

0 

1 


0 

0 

-1/3 

0 



0 

0 

1 

0 


0 

0 

1 

0 


0 

0 

0 

-1 

9 

\[e 

- 




2 , ^ \/e 

_ 



■ _ 

2 , , \l6 

r— 




To = 

2 

1 

0 

0 

0 

Tj^(a) - 

0 

-1 

0 

0 

Ti{s) = 

0 

-1 

0 

0 



0 

-1 

0 

0 


-1 

0 

0 

0 


1 

0 

0 

0 



0 

0 

-1 

0 


0 

0 

0 

1 


0 

0 

0 

1 



0 

0 

0 

1 


0 

0 


0 


0 

L- 

0 

-1 

0 

CM CM 

re 

T 

0 

0 

1 

0 

T^(s) = ^ 

0 

0 

1 

0 

= ' 3 

2^0 

1 

0 

0 

0 



0 

0 

0 

1 


0 

0 

0 

1 


0 

-3 

0 

0 



-1 

0 

0 

0 


1 

0 

0 

0 


0 

0 

3 

0 



0 

-1 

0 

0 


0 

1 

0 

0 


0 

0 

0 

-1 

T^(a) 

-3 

0 

1 

0 

0 

tJ(s) = 

2 Vs 

0 

■ -1 

0 

0 

TjCa) =3^ 

0 

0 

-1 

0 



1 

0 

-./3 

0 


1 

0 

\/3 

0 


0 

0 

0 

1 



0 

-n/3 

0 

1 


0 

^f3 

0 

1 


1 

0 

0 

0 



0 

0 

1 

0 


0 

0 

1 

0 


0 

-1 

0 

0 

T^^s) 

3 

2 \/2 

0 

0 

1 

0 


0 

0 

0 

1 

T^(s) = 1 

0 

0 

0 

-1 



0 

0 

0 

-1 


0 

0 

0 

0 


0 

0 

0 

G 



1 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 



0 

-1 

0 

0 


1 

0 

0 

0 


1 

0 

0 

0 
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[b]3 = b[b]2, [b]" = [[b]2]2 

[B]^ = B[B]^ 

2 2 

where = a + bb* 

2 2 

k 2 = c + dd* 


and 


k^ = [cos^Q (2+p)^ + {(^-1)^ + T}^ + n Gos20'}sin^e]^'^^ 

.. (III. 32) 

k 2 = - 2 ^ [cos^e (f-2)^ + + T}^ -2(y+l)rj cGs2i2r} sin^G] 

. . (III. 33) 


can be written in the matrix form, with the matrix elements 
given as follows: 



= COS 

3 (w^t) - i ^ sin(co^t) 


= 

, (t) ]l3 = 0 

[“z<'^di4 

= -i 

=in(<4At) 


= 0 


[Uz(t)]22 

= cos (Wgt) - i ^ sin(cOgt) 

[^2^^d23 


— sin(a)gt) 

C‘'z‘'^d24 

^ 0 


[“z‘«]31 

=5 0 


[U^(t)]32 

= -i 

(3.^ 

^ sin(a)gt) 


. . .contd 
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= cosCWgt) + i ^ 

sin(cOgt) 




= 0 






= -i 

'K'A' 

— sin(a)^t) 





= 

z‘t)J43 = 0 




[^'■^>144 

= cos (co^t) + i ^ 

sin(w^t) 

.. (Ill 

.34) 

where oo^ 

= k, 

0 

O 

11 




“b 

II 

to 

Wb 

~ T~ 
o ^ 


.. (Ill 

.35) 


Also, it can be seen that U U =1 • Having obtained U in the 
form of matrix, it is now straightforward to arrive at the time 
evolution of tensor operators . In calculating the evolution of 
tensor operators the expressions of the form 

have been evaluated by .multiplying the corresponding matrices . 
The resulting matrix has been converted back into equivalent 
tensor operators. The evolution of a few selected tensor 
operators obtained in this manner are tabulated in Table III. 7. 
Having obtained the evolution of tensor operators under the 
action of and the calculation of pulse responses is now 


s traightf orward . 



Table III. 7 


U T^(s)u:^ 


U T2(a)u;^ 
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. The Evolution of a few Selected Tensor Operators 

c3. 

Under the Influence of the Hamiltonian •J( 

= T^.(s) cos.|((0^-c0j^)t] + A^2 cos [|( (0^+0)^) t]] 

+ T^(a) [A22 { cos[-|(a)_^-u^) t] - cos [■|(a)^+oo^) t ] 3 

A22 sin[|(Q^+co^)t] + A23 sin[i(w^-co^)t]] 

+ T2(a) [A32 { cos[|(a)^-Wj^)t] - cos [|( t] } 

+ A32 1 sin[|(o^+a)j^) t]3 + A33 { sin['|(u^-w^) t])] 

+ T2(s) [A42 i cos[|(co^-(o^)t} -cos[|(w^+co^)t ]3 

+ A42 sin[§((o^+a)j^)t] + A^3 sin[|( t]] 

= T^(s) { cos[-|((o^-(^)t] -cos[|(w^+w^)t ]3 

+ B22 { sin[|(Wg^+(^) t]} +B23 { S'in[|(a)^-a)^)t]3] 
t^21 cos[|(w^-a>j^)t] + B22 cos[|(a)^+6>j^)t]] 

+ T2(a) [B3^ { cos [|(co^-(^)t] -cos [|( 0)^+0)^) t]} 

+ B32 sin[| (Wa+0)j^)t] + B33 sin[| (<o^-(a^)t]] 

: : ^ ®42 ^ ^ ®43 




ro to 
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Table I II . 7 (contd.) 

(a)U 2 ^ = T^(s) [0^2 {cos[^ - cos[i ( co^+Wj^) t] } 

+ D^2 sin [| (o)^-o)^>t]] 

2 1 ’ 1 . 

+ T^Ca) [D2]_ {cos[- 2 (w^-( 05 ^)t] - cos['^ (co^+co^)t]} 

1 1 
+ D 22 sin['j (w^+a)j 3 )t] + D 23 sin[-^ 

2 1 * 1 
+ T2(a) cos['^ "*’ ^32 °°s[^ ( cOg^-Wj^) t]] 

+ ^ 2 ( 3 ) { COs[|- (c 0 g+C 0 b)t] - COs[i (Qg- 03 ^)t]} 

+ D 42 Sin[ I ( 0 )^+w 5 ^)tJ + D 43 sin[| 

= T 2 (s) [ 0^2 {cos[^ “ cos[| (Ug^+w^)t]j 

rsj 

+ G 22 sin[| (w^+Wt,)t] + G ^3 sin[| (a)^-o)^)t]] 

+ T^Ca) [G 2 3 ^ { cos[|- (a 3 _^+w^) t] - cos[i 
*>/ 

+ G 22 sin[| ( 0 )^+ 0 )^) t] + G 23 sin[| (o)^-c^)t]] 

2 1 1 
+ T 2 (a) [G 32 { cos['j (cOg^-cOj^) t] - cos [2 

+ G 32 i sin[| (w^+(^)t]} +G 33 sin[| 

+ 12 ( 3 ) [G 42 i cos[| (w^+(^)t]} +G 42 cos[| (a)^-co^)t]] 



Table III. 7 (contd.) 


where = exp ( -i t) 

d 

The Hamiltonian "K. is defined in Table III. 4. 

z 


^11 

^^1^2 

+ 1 + ac)/2 k^k2 

^12 

= [k.k^ 

- (l+ac)]/2 k^k2 

^21 


= m/2 k^k2 . 

^22 

= ®12 = 

D43 = 0^3 = i(ck^ -ak2)/2 k^k2 

'^23 

^ 1 3 

°42 "" ®32 "" -i(ck^ + ak2)/2 k^k2 

\ — 1 
ro 

^11 

(cs^- aq^)/2 k^k2 

^32 

"" ~°12 "" 

1 ( 32^2 -q2k2)/2 k^k2 

^33 

"" “^13 " 

i(s2k2 -l-q2k^)/2 k^k2 

^41 

= = 

Ccs2 - aq2)/2 k^k2 

^42 

^ ®33 " 

'^12 “ "^23 ^ 

^43 

“®32 ~ 

Gi 3 = -D22 = i (q^k^ +s^k2)/2 k^k 

®21 


+ ac-l)/2 k^k2 

®22 


-(ac-l-))/2 k^k2 

^31 

"°21 "" 

[-(aq2 +CS2) ]/2 k^k2 

®41 

" ^21 " ■ 

-(q^a +s^c)/2 k^k2 

®42 

^ "®22 "" 

-i(q 2 k’^ +32^2)/^ k^k 2 

®43 

^23 

-i(q2k^ -S2k2)/2 k^k2 

^31 

“ ^^1^2 

+ ca + s^q^ + S2q2)/2 k^k2 

°32 

— |^}c2k2 

- (-ca + S2q]_ +S2q2)]/2 k^k2 

°41 


^31 

- (s^q2 

+ S2q]_)/2 k^k2 

^41 

= [k^k2 

- (-CS, tSj^q^^ +S 2 < 32 ^]/^ ^ 1^2 

^ 42 ' 

= [^2^2 

- ca + s^q^ + S2q2]/2 k2^k2 


. . .contd 
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Table III»7 (contd.) 


with 


^2 “ 

1 = 

m = 

■q. =_ 


^ [cos^0(2+f)^ + {(f-1)^ +7)^ +2(?-l)r] cos2jZfJ 
[cos^0(^-2)^ + {(f+1)^ + r}^ -2il+f)n COS20'} 

1 rr .2 2 ,, .2 2 ,, ,2 2-1 

^ [p3_{ -(1-77) } +P 2 { d+r}) -f }] 

p p . 0 — 

^1^2 2f 


1 =-27 

“ ^ ^1 3 

“ ^ ^2 -d+ri )] 

P2/ P3 and j’ are defined in Table III. 2 


^2 

S2 
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1 1 1 . B . 5 Evaluation of Two-Pulse Response; ZSEEM Pattern 

We are interested in calculating the two-pulse response 
of non-interacting nuclei with spin I = 3/2 in polycrystalline 
samples in the presence of a weak Zeeman field. In all the 
calculations the relaxation effects and efg inhomogeneities are 
not considered. 

We assume a configuration in which the Zeeman and r.f. 
fields are oriented parallel to each other and their axis is 
oriented randomly with respect to the principal axes of efg 
tensor. The total Hamiltonian for this situation is given in 
Egn . III. 13. In the quadrupolar interaction representation 
(see Eqn. III. 25) the total Hamiltonian is given by Eqn. III. 26. 

As mentioned in the previous chapters of this thesis the 
two-pulse response in this oase also. can be obtained by follow- 
ing the time evolution of the density matrix under r.f. and 
Zeeman interactions in QIR., As the magnetic field is assumed to 
be weak, i.e., Wq .the thermal equilibrium reduced density 

matrix (see Chapter I, Section I.D.l) for this case also is given 
by 

p(0) = T^ . . (III. 36) 

AJ O 

In QIR then 

P(0) = tJ . . (III. 37) 

In the following calculations we neglect the effect of the 
Zeeman field during the time when pulse is 'on ' , assuming that 
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causes negligible precession of the spins during applica- . 
tion of the pulse. 

Now, consider the response of nuclear spin I = 3/2 system, 
in the presence of a weak Zeeman field, to the two— pulse sequence 
shown in Fig. III. 4. 


( 7 t / 2 )^ {n/2)^ 




/N 





1 ! 

u ^ 

i ! 

t 


t^^ t 

w w 

Fig. III. 4, Tv;o Pulse Spin-Echo Sequence 

Immediately following the application of r.f. pulse the density 
matrix is given by (see Table III. 3) 

P(t^^,) = cos£ + i sinf [c^ ■** *^2 ^3 5*2 ^ 3 

.. (III. 38) 

where £ , C^, C 2 and are defined in Table III. 3. 

Following the removal of the pulse the spin system evolves 
under the influence of Zeeman interaction and the density matrix 
is given by 

P(t) = P(t^j^) U"^ .. (III. 39) 

= exp (-i T ) 


v/here 


. . (III. 40) 
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Similarly, the density matrix after the second pulse can be 
written as 


p(t-T) = U^(t-T) Uj^(t^^) U"^(t^,) U"^(t-T) 

where Uj^(t^j) = exp (-i'K^ t^,) 

■ exp [-i' 3 i^ (t-x)] 


(111. 41) 

(111. 42) 

(111. 43) 


Now, in any experiment, we can only measure the magnetization 
and it can be obtained by considering the expectation value of 
the magnetization operator [2l]. The signal is proportional to 
the expectation value of the magnetization operator. The roagne-^ 
tization operators (at the resonance frequency) in phase and 90*^ 
out of phase w.r.t. the applied r.f. pulses, in the laboratory 
frame can be written, respectively as 


and 


m' = 2 cos(w t) (-Tha) T^s) + P. tS 

X Q rO 1 • X Z ± o r\J O 

(in units of r ti) .. (III. 44) 

= -2 sin(cogt) (-T^(a) Pj + P2 + P3 tJ) 

(in units of fh) .. (III. 45) 


In QIR these are given by 


rs» I 

»x 


= 2 cos(«gt) CUqUj [-Tj(a) P 3 . + t1(s) Pj + P 3 ] 


-1 -1 

U, U } 
d q 


contd . 
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- [-P i I T];(a) + T^Ca)} 

L 1 5 .1 yfS 


+p^ , I tJ(s) h. ^ tJ(s)) 

■^^3 i - > ^^">1 


.. (III. 46) 


and 


nJ 

M 


X 


-2 sin fat) iU [-Ti(a) + tJ(s> P 3 ] 


— 1 —1 

U, U 1 
d q ^ 


^/2 

i 

/2 


- [- P, T^s) 3- P 3 ^ 


(3-n) 


^ T^s) + Cj T^Ca) + T^(a)] 


f T^Ca) -P3 §a T^ta)] 


. (III. 47) 


where = 


^5 = 


^6 = 


.p ■( 3 + 77 ) 
1 3f 

p ( 3 - 77 ) 
^2 3f 

P In 
■ 3 3j> , 


.. (III. 48) 


In arriving at Eqns . III. 46 and III. 47 we have made use of 
Tables Ill.l and 11.12 and dropped the high frequency terms. 
NOW/ making use of Tables III. 3 and III. 7 and Eqns. III.41 and 
III. 47 and considering only the terms that contribute to echo, 
the expectation value of the magnetization operator M can be 

-/C 

written as 

rJ 1 1 

^ " "rtr ^^1 ^Wa+w^)t] + P2^cos[^ (w^-COj^yt] 

+ F 3 sin[-| (w^+coj^y.t] + sin[-| (w^-w^)t] 

. . .contd. 
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1 *1 
+ cos[-| (CO^+O)^) (t-2T) ] + Fg cos[j (t-2T) ] 

1 1 
+ F^ sin['^ (co^ +w^)(t-2T)] + Fg sin[^ (w^-coj^) ( ■t-2T ) ] 

4- Fg COs[| (0)^+W^)T + (t-T)] 

1 1 

+ F^g cos[^ “ 2 ^ ^ 

11 

+ F^^ cos[^ **’2 (w^+w^) ( t-T ) ] 

+ F^2 (Wa+%)(t-T) - I 

+ F ^3 sin[i (w^-u^)(t-T) + I (w^4-w^^)t] 

+ F ^4 sin[| (u.^-(^){t-T) - I (a)^+c^)T] 

1 1 

+ f^ 3 sin[^ (w^+co^) (t-T) + ■^ (co^-Wj^)t]' ■ 

+ F^g sin[-^ (cOg^+cO]^) ( t-x ) - ^ (Wg^-caj^)'^ ] } (III. 49) 


where F's are as given in Table III. 8. Similarly, it can be 

easily shown that the expectation value of the in-phase compo- 

/ 

nent of the magnetization, i.e., M ) is zero. 

The magnetization corresponding to the echo at t = 2T is 
obtained by setting t = 2x in Eqn. III. 49. Thus, 

/aw 

£(2t) Sins {P^ cos [ (oj ] + P^ cos[ ( ] 

+ Pg sin[ (co^+o)^)t] + F^ sin[ 

+ (Fg+F^^) cos(w^x) + cos(Wj^x) 

sinCco^T:) + (F^_^+Fj^g) sin ( cOj^x) + (Pg+Fg) j 

ill . . (111.50) 



Table 111.3. Expressions for F's 


1 68 

(X^-cos^ t./^ (Y^-cos^ £ /2 

+ Cg (Z^-cos ^£/2 

P2 = C4 (X2-cos^ £/2 E^2^ + ^5 (Y2 -cos^£/2 E22) 

+ Cg (Z2 -cos^£/2 E32) 

F3 = (X^-cos^ E ./2 E33) + Gg (Y3 -cos^ £ /2 £33) 

+ Cg (Z 3 -cos 2£/2 E 33 ) 

F4 = C4 (X^-cos^ £/2 E^^) + C3 (Y^-cos^ £/2 

+ Cg (Z^-cos^ £/2 £ 3 ^) 


^5 

— 

C4X5 + 

C5Y5 

+ 

^ 6^5 

^6 

= 

^ 4^6 

^ 5-^6 

4 - 


^7 

= 

C4X7 + 

C5Y7 

4 - 

^6^7 

^8. 

= 

^ 4^8 + 

<^ 5^8 

4 - 

<=6^8 

^9 

= 

C4X5 + 

CsYg 

4 - 

‘^6^9 

^10 

= 

^ 4^10 

+ C5Y 

10 

^6^10 

^11 


^ 4^11 

+ C5Y 

11 

+ CgZ^^ 

^12 

= 

C4X1 2 

+ C5Y 

12 

^6^12 

^13 


^ 4^13 

+ C5Y 

13 : 

>• - 6^13 


contd 
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Table III.B- (contd.) 


^14 = 

C4X14 

+ 

^ 5^14 

+ 

^^ 6^14 

^15 = 

C4X15 

4 * 




^16 = 

^ 4^16 

4 - 

^ 5^16 


^ 6^16 


where = 


Pj^O+T]) P2(3-77) 


-2 77 


E is defined in Table 111.3. 

P 2 , P^ are defined in Table III. 2. 

E's, X's, Y's and Z's are defined, respectively, 
in Tables III.9, 111.10,111.11 and III. 12, 
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Table Expressions for E's and H's 


^ 1^12 “ " ^ 3^11 




*^ 2^12 ^ 3^12 


‘“ 2^13 ^ 3^13 


■C1A21 ■*■ ^2^22 " ^3^21 


C1A21 + C 2 D 2 ^ + C2B22 


C 1 A 22 + C 3 D 22 


'^ 1^23 ^ 3^23 


■^ 1^31 “ '^ 2^31 ^ 3^31 


*^ 1^31 ■*■ *^ 2®31 ^ 3^32 


C1A32 + C2B32 


*^ 1^33 *^ 2®33 


■^ 1^41 " ^ 2®41 ■'■' ^ 3^41 


C1A41 + C2B42 - C3D4, 


^ 1^42 ■*■ ^ 2®42 ^ 3^42 


^ 1^43 ■'■ ^ 2®43 ^ 3^43 


= C^Aj2 - =^22 + =^31 

^1^2^“®11 ■*■^21^ + C^C3(-D2^^ -A3^) + C2C^iT:>2^ -B^^) 


. . .contd. 



Table 1 1 1 , 9 (contd.) 


17 1 

»21 = =1*11 = =^21 + =3°32 

P^O+r]} P2(3-r]) 2'n P 3 

^ '/a f X " ^ \/3fx " Jsj’x 

X = [ — -- 2 {sin^G [( 9 + 77 ^) - 6 77cos 20'][ + 47)^ cos^©}]^'^^ 

(3+7) ) 

where A's, B's, D's are defined in Table III, 7, 

while P^, p 2 , P 3 and f are defined in Table III. 2. 
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Table III.IO. Expressions for Q* s and X's 


= (cOS£-l) 
□ 2 = (cos£-l) 

= (cos£-l) 
= (cos£-l) 
= (cosE-l) 
Qg = (cos£-l) 
= (cos£-l) 
Qg = (cos£-l) 


- "33 

^2 ^ ^6 ” ^44 

X 3 = Q 3 + I 31 

^4 = ^8 + ^42 

X 5 = + I 33 

^ 6 . = ^6 "44 

^7 “ *^3 “ ^31 



* Hi2>] + ^11 

[^ 12(^11 + H^j)] + I 12 

+ Hi2>] + ^13 
[Ei 4CH^2 Hu>] b4 
- » 12 )] + ^21 
[Ei 2 (H 21 - H22)] + I22 
[E 23 (H 2 i - H22)] + I23 

Pi4<H21 - % 2 )] + I 24 


^9 = 

= Q 2 - 

- I 

'34 


= Q 2 

+ 

^34 

Xii 

= % 

- 

^43 

^12 

= % 

+ 

^43 

^13 

= Q 4 


^32 

^14 

= Q 4 

- 

^32 

^15 

Q 7 


^41 

^16 

= ^7 

- 

^41 


contd 
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Table III . 1 o(contdj 


^11 

~ 

^11^12“^21^1l“®31^1l"^4l'^ll 

^12 

= 

® 1 1 2 1 ® 1 1"*"^ 3 1^1 1'^®4 1^ 1 1 

^13 

= 

®21®12‘*'^31°12'^^41*^12 

^14 

= 

^2 1 ® 1 3 3 1^ 1 3 ■*■^4 1*^ 1 3 

^21 

= 

^12®12"^22®11“®32^11"^42^11 

^22 

= 

^ 1 2^^ 1 1'*'^ 2 2 ^ 1 1 3 2^ 1 1‘*'^4 2*^ 1 1 

^23 

SS 

^22®12'*‘®32^12'*'®42°12 

^24 

= 

®22®13'‘‘^32°13’*'®42°13 

^31 

= 

®13^12"®23®11"®33^11"^43'^11 

^32 

= 

® 1 3^ 1 1'*'^ 2 3 ® 1 1"*"® 3 3^1 1'''^4 3^ 1 1 

^33 

= 

®23®12'*’^33°12"^^43^12 

^34 

= 

^23®13'*'®3 3^13'*'^43^13 

^41 

= 

®14^12"^24®11“®34°11“^44^11 

^42 

= 

® 1 4^ 1 1"*"^ 2 4 ® 1 1"*"^ 3 4^1 1'*'® 4 4 '^ 1 1 

^43 

= 

^24®12'*‘^34^12'‘'^44^12 

^44 

= 

®24®13'*'®34°13'''^44^13 


E is defined in Table III. 3. 

A's, B's, D's and G's are defined in Table III. 7. 
C's, E's and H's are defined in Table III. 9 . 
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Table III. 11 . Expressions for R's and V's 

= (cos£-l) [£21^^11 ^12^1 ‘^11 

R2 = (cos£-l) [E22(H^1 + H^2^1 ^12 

R3 = (cos£-l) [E23(H^^ + H^2^] >^13 

R^ = (cos£-l) [E24(H^^ + H^2>j+ ^14 

R5 = (cos£-l) [E2 i(H2;^ - H^2)] + J21 

Rg = (cosa-1) - H22)] + J22 

Rry == (cos^ — 1) t^23^^21. *” ^12^^ 2 3 

Rg = (cos£-l) [E 24 ‘H 21 - “l2>] + ^24 



— 



"^33 

VD 

If 

^ ^2 - 

• J 

34 

^2 

= 

^6 

- 

>"44 

^10 

^2 

+ 

■^34 

^3 

= 

^3 

+ 

'^31 

^11 

= ^5 

- 

^43 

^4 

= 

^8 

+ 

^42 

^12 

= ^5 


^43 

^5 

= 


+ 

*^33 

^13 

= ^ 4 . 

+ 

■^32 

^6 

= 

^6 

+ 

'^44 

^14 

= R4 

- 

■^32 

^7 

= 


- 

'^31 

^15 

= R3 

•f 

^41 

00 

= 

^8 

- 

'^42 

""16 

= R3 

- 

J4I 


Gontd 




175 


Table III. 11 (contd.' 


where J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 


11 

— 

"^ 11^21 

+ 

^ 21®22 

1 

^31^21 

— 

^41*^21 

12 

= 

^ 11^21 

+ 

^ 21 ® 2 l 

4 

^31^21 

4 

^41®21 

13 

= 

^ 11^22 

+ 

^31^22 

4 

^41*^22 



14 

= 

^ll'^23 


®31^23 

4 

^41^23 



21 

= 

“® 12^21 


® 22®22 

} 

^32^21 

- 

^42^21 

22 

= 

^ 12'^21 

+ 

® 22®21 

4 

^32^21 

4 

^42^21 

23 

= 

^ 12^22 

+ 

^ 22^22 

4 

^42*^22 



24 

= 

^12^23 


^32^23 

4 

^42*^23 



31 


^13-^21 

+ 

^23®22 

- 

^33^21 

- 

E 43 G 21 

32 


^13-^21 

4 - 

®23®21 

4 

^33^21 

4 

^43^21 

33 

= 

^13-^22 

4 

^33^22 

4 

®43®22 



34 


^13*23 

4 

®33°23 

4 

®43®23 



41 

= 

®14'^21 

4 

^24®22 

- 

®34°21 



42 


^14^21 

4 

®24®21 

4 

®34^21 

4 

^44^21 

43 

= 

^14^22 

4 

^34^22 

4 

^44^ 2 



44 

= 

^14^23 

4 

^34^23 

4 

^44^23 




£ is defined in Table III. 3. 

A's, B's, D's and G*s are defined in Table III. 7. 
C's, E's and H's are defined in Table I II, 9. 
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Table III. 12. Expressions for S's and Z's 


= (cos£-l) 
S2 = (cosE-1) 
= (cos£-l) 
= (cosf-1) 
Sg = (cosf-1) 
Sg = (cos£-l) 
= (cos£-l) 
Sg = (cosE-1) 


[EbiCHii + H^ 2>] + ^11 
[E32(H^1 + + K^2 

[EgaCH^l + H^2>]+ ^^13 
[E34(Hii f 1 - 1 , 2 )]+ K,4 
[E3i(H2i - H, 2)]+ K21 
[E32(H2i - H, 2)]+ K22 
[E33(H2i - H, 2 )] + K 23 

CE34(H21 - H,2)] + K24 


Z, - S, - K33 

^2 = ^2 - K44 

Z3 = S3 -H K3, 

Z4 = ^8 ^42 

^5 ~ ''' ^33 

^6 = ^6 + ^44 

Z7 = S3 - K3, 

^8 “ ^8 ” 


^2 - ^^34 


'10 


^2 + ^34 


'11 


^5 " ^43 


'12 


Sr + K . -3 
5 43 


'13 


^4 ^32 


'l4 


^4 - ^32 


'15 


S7 + K4, 


'16 


S7 - K4, 


o .contd 
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Table 111,12 (contd.' 


where K 


^11 

— 

"^11^31 

- 

^21®31 ^31^31 

- 

■ ®41°31 

Ki2 

= 

^11^31 

+ 

^2i®31 

4* 

^31^31 

4“ 

■■'41*^' 31 

^13 

= 

^ll'^32 

+ 

®21®32 

4* 

^41*^32 



^14 


^11^33 

4- 

^21®33 

+ 

^41*^33 



K21 

= 

^12^31 

- 

^22®31 

4- 

E32D31 

- 

^42^31 

K22 

= 


4- 

^22®31 

4- 

^32^32 

4- 

^42^31 

^23 

= 

^12^32 

+ 

^22^32 

4- 

^42*^32 



^24 

— 

^12'^33 

4- 

^22®33 

4- 

^42*^33 



K31 

= 

13^31 

- 

^23^31 

4- 

^33^31 

- 

®43°31 

^32 

= 

^13^31 

+ 

^23®31 

4- 

®33^32 

4- 

^43^31 

^33 

= 

^13^32 

4“ 

^23®32 

4- 

^43^32 



^34 

= 

^13^33 

4- 

^23^33 

4“ 

^43'^33 



^41 

= 

-^14^31 - 

®24®31 ^34^31 

- 

■ ^43^31 

K42 


^14^31 

4- 

^24®31 

4- 

^34^32 

+ 

^43^31 

^43 

= 

^14^32 

+ 

^24®32 

4- 

^43^^32 



^44 

= 

^14^33 

4- 

®24®33 

4- 

^43*^33 




£ is defined in Table III. 3. 

A's, 3's, D's and G's are defined in Table III. 7 
C's, E's and H's are defined in Table II I. 9 . 
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From Eqn. III. 50 it is evident that the expression for the Zeeman- 
perturbed spin echo envelope modulations contains all the frequency 
components ( 

Ramachandran and Narasimhan. 

In the following special cases the expression for (2^) 
reduces to those given in [l] (also see Chapter II) . 


(o^+o^b) and (VC^) 


as reported earlier by 


Case 1 : For 'q = 0, = 0, and = 


■+2cj^ cos(co^t) T^(a) 


‘^a “ “b “ 


C (2t) a sin£ sin^£/2 


.. (III. 51) 


Case 2: 


For n = 0, = fn(e) and = 2 0)^ cos ( w^t) Tj(a) 

f = (1+4 tan^e)^^^ 

2 

£(2'^)asin£: sin^ (£/2) [ ( " ■■■ Z^ ) { cos(co t) + cos(co,t) 

2f^ a h 

-COS(cO^T) cos ( 60 ^ 7 ^)} 

,f+1.2 

^' 21 ^ J 


.. ( 111 . 52 ) 


III.B.5(i) Powder averaging: 

Since in the above calculations we have assumed that the 
r.f. and Zeeman fields are randomly oriented with respect to the 
efg principal axes system, the expression for £(2t) is also 
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applicable for this situation. The response of a quadrupolar 
spin ensemble in a polycrystalline specimen, where all values of 
0 and gT occur with equal probability is obtained by considering 
a weighted average overall orientations, 

2nn 

<M^> a J J £(2t) sin© d© d0' .. (III. 53) 

o o 

£(2t) obtained using the tensor operator formalism (Eqn. In.'S©) 
is inserted in Eqn. III. 53 which has to be then evaluated to 
obtain Zeeman-perturbed spin echo envelope modulation pattern 
from polycrystalline samples. 

From the above presentation it is clear that the tensor 
operator formalism can be used for the theoretical analysis of 
ZSEEM in spin I = 3/2 nuclei in polycrys talline specimens. How- 
ever, before concluding this section, the following points may be 
mentioned. As pointed out earlier in this section in the calcu- 
lation of the expression which gives the ZSEEM pattern from poly- 
crystalline specimens one needs to evaluate the time evolution 
of tensor operators under r.f. and Zeeman interaction Hamiltonians 
corresponding to general (©, 0) orientation. The calculation of 
the time evolution of tensor operators under the action of the 
r.f. Hamiltonian is straightforward. However, it is rather 
difficult to obtain closed- form expressions for the time evolu- 
tion of tensor operators under the Zeeman Hamiltonian by any of 
the. two methods viz . , the method of 'nested commutation 
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relationships' or the 'harmonics of the motion method'. Hence^ 
we had to resort to the matrix representation method for evaluat- 
ing the time evolution of the operators and reconvert the result- 
ing matrix into the tensor operators . On this account, the cal- 
culation became rather cumbersome. The tensor operator formalism 
per se does not seem .to offer any additional advantage over the 
other approaches for the analysis of ZS'EEM. 

SUMMARY 

Tensor operator formalism has been applied for the theore- 
tical analysis of Zeeman-perturbed spin echo envelope modulations 
in spin I = 3/2 nuclei in polycrystalline specimens. The present 
results are in agreement with those obtained earlier by Rama- 
chandran and Narasimhan. The tensor operator formalism, however, 
does not offer any special advantage over the other approaches on 
account of the complexity in the calculation of the evolution of 
tensor operators under the action of Zeeman Hamiltonian. 
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The NMR spectrum of a solid sample is usually broady mostly 
due to dipole-dipole interactions in solid-state masking inform- 
ation about chemical shifts and spin couplings. To achieve high 
resolution NMR spectra from solid samples, various experimental 
techniques have been developed, many of which are. based on multi- 
ple r.f. pulses- [ 1 , 2 ]. The first multiple-pulse experiment in 
NMR of solids was due to Ostroff and Waugh [3] and Mansfield and 
Ware [4]. These workers observed that the application of modi- 
fied Carr-Purcell (CP) sequence leads to a lengthening of the 
decay of transverse magnetization. Since then a variety of multi 
pie pulse sequences have been designed to achieve prolonged 
transverse magnetization. These sequences not only achieve 
long transverse magnetization but also give rise to selective 
averaging thereby leaving chemical shift interactions active 
and leading to the high resolution NMR spectrum from solids . 

Multi pie- pulse experimental techniques in NMR of solids are now 
well established [1,2]. 

The interest in the study of multi pie- pulse experiments 
in NQR stems from the following considerations. Multi pie- pulse 
sequences in NQR could also lead to prolonged transverse magne- 
tization and thus to a narrowed line. Information on various 
internal interactions which broaden or split the line can be 
then obtained. This may help in a better understanding of the 
chemical environment of the quadrupolar nucleus in the solid 
state. Relaocation times in rotating reference frame may also 
be obtained from a study of NQR multiple pulse experiments. 
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Multiple pulse experiments can be effectively used iri achieving 
enhanced S/ISI. Under the multiple r.f. pulse irradiation a sensi- 
tive double resonance experiment can be performed to detect rare 
quadrupolar nuclei. Thus, the multiple pulse investigations in 
NQR could yield valuable structural and dynamical information in 
the solid state. 

This chapter consists of four sections. The first section 
gives a brief review of the multiple pulse investigations in NQR. 
Section IV. B deals with the experimental investigations of a 
variety of multiple pulse sequence on spin 3/2 and 5/2 nuclei. 
Section IV. C gives the results of a double resonance experiment 
under multiple pulse spin locking condition. The last section 
outlines theoretical analysis of a few multiple pulse experi- 
ments in NQR of spin 3/2 nuclei using tensor operator formalism. 

IV. A MULTIPLE PULSE INVESTIGATIONS IN PURE NQR - A REVIEW 

The first experimental investigations of multiple pulse 

sequences in NQR were due to Marino and Klainer (MK) [5], They 

14 

employed the Ostrof f-Waugh (OW) sequence to study N (I = 1) in 
a polycrystalline sample of NaN 02 and found effects very similar 
to those observed in NMR experiments, namely, a persistence of 
the transverse magnetization for times long compared to, T 2 / 
spin-spin relaxation time. Following this, Osokin [6-9] inves- 
tigated the effect of phase alternated multiple pulse sequence 
(PAPS), Carr-Furcell-Meiboom-Gill (CPMG) and OW sequence on NUR 
of spin I = 1 nuclei. ■ Ainbinder et al . [lo] have studied the 
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3S 123 

influence of CPMG sequence on Cl (I = 3/2) and Sb (I = 7/2) 
in polycrystalline samples. Recently, Narsimha Reddy and Narasim- 
han [ll] have carried out experimental investigations of the 
influence of CP, CPMG, OW and PAPS sequences in the case of spin 
I = 3/2 nuclei in polycrystalline samples. In all the above men- 
tioned investigations, the persistence of the transverse magne- 
tization for times much longer than T 2 was observed. Recently, 
Ermakov and Osokin [12J studied the influence of the well-known 

pulse sequence, WAHUHA [13^, which is known to suppress the 

14 

dipole-dipole interactions in NMR of solids, on N (I = 1) NQR 
in a single crystal of NaN 02 . 

Apart from the above experimental investigations there are 
some theoretical studies in the area of multiple pulse NQR. Using 
different approaches attempts have been made to give a theoretical 
explanation for the long time and short time behaviour of the 
magnetization in the MK experiment [14-16]. 

Cantor and Waugh [14] employed an approximation method to 
theoretically investigate the short time behaviour and to some 
extent long time behaviour of the magnetization in MK's experi- 
ment. Hitrin et al, [15] employed a canonical transformation 
technique to obtain solution for the density matrix, and obtained 
theoretical agreement with the experimentally observed short and 
long time behaviour of the magnetization. Matti Marick [16] used 
a generalized version of the Average Hamiltonian Theory (AHT) 
based on Floquet solution of the periodic, time-dependent, 
Schrttdinger equation to analyze the MK experiment and compared 
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the results with those obtained in the case of NMR. 

In all the aforesaid cases the theoretical investigations 
were carried out on spin-1 nuclei only. Recently, Ainbinder 
et al . [17,18] have developed a theory which is applicable for 
general spin I. They employed the method of canonical transfor- 
mations and Krylov- Bogolyubov-Mitropolski averaging method [19]. 
They have shown that in the multiple pulse sequence, (^/2) - 
[t- referred hereafter as spin- locking sequence (SLS) , 
spin locking occurs for arbitrary 0 values when 0 = (2k+1)^/2 
where K is an integer. Ainbinder and Furman [17] proposed that 
like in NMR, multiple-pulse experiments in NQR spectroscopy lead 
to two types of averaging of internal interactions which broaden 
the NQR line. These are termed as: (i) general averaging, and 
(ii) selective averaging. A pulse sequence that influences the 
broadening interactions as a whole, regardless of their nature, 
is said to cause general averaging. 

The influence of relaxation processes on the quasi-equili- 
britun state of the magnetization in the spin-locking experiments 
has also been investigated in the literature [lo]. 

The theoretical investigations of PAPS on spin 1=1 
nuclei have also been carried out in the literature [6-8]. In 
these investigations AHT has been employed. -Ermakov and 
Osokin [12] have investigated the influence of PAPS on a system 
of two dipolar coupled equivalent spin 1=1 nuclei and predicted 
the quasi -equilibrium state of the magnetization. Osokin [9] 
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studied the influence of spin-lattice relaxation process on the 
the quasi -equilibrium state of the magnetization in PAPS experi- 
ment. 

It may be mentioned that in a multiple-pulse experiment 
an internal interaction is said to be averaged out if the average 
Hamiltonian of ' that interaction over the cycle of the pulse 

sequence satisfies the condition [20] 

0 or ]= 0 

where is the initial density operator prepared by the prepa- 
ratory pulse. Employing the "fictitious spin- 1/2 operator forma- 
lism" [21] Osokin [6-8] has shown that under the influence of 
PAPS, the average Hamiltonian corresponding to the internal 

interactions due to (i) inhomogeneities in efg's, (ii) torsional 
oscillations of molecules, (iii) heteronuclear dipole-dipole 
interactions commutes with the initial density operator, P^. 

Thus, he concluded that under the action of PAPS the above 
mentioned Internal interactions are averaged out. Osokin [7,8] 
has also shown that the homonuclear dipole-dipole interaction is 
not averaged out under the action of PAPS. Theoretical analysis 
of some of the multiple pulse sequences in NQR have been considered 
earlier by Zueva et al. [22]. Using the fictitious spin- 1/2 
operator formalism Ermakov and Osokin [12] have analyzed the 
influence of WAHUHA sequence on a model system of two dipolar- 
coupled spin-1 nuclei. 


int 



In the next section we present our experimental results 
of some multiple pulse sequences in NQR of spin I = 3/2 and 5/2 
nuclei. 

IV.fi EXPERIMENTAL INVESTIGATIONS OF RESPONSES OF QUADRUPOLAR 
NUCLEI TO MULTIPLE PULSE SEQUENCES 

As has been discussed earlier, under the action of multi- 
ple pulse sequences some or all contributions of the interactions 
which lead to line broadening can be removed and thus one can 
achieve a narrowed line. This kind of resolution enhancement, 
in principle, enables one to study the weaker interactions in 
solid state. 

In this section we present our experimental investigations 

of the influence of a few multiple pulse sequences on quadrupolar 

35 

nuclei. The nuclei investigated are: (a) Cl (I = 3/2) in a 

single crystal of NaClO^, and in polycrystalline samples of 

187 185 

NaClO^/ KCIO^ and SbCl^. (b) Re & Re (both 1 = 5/2) in 
polycrystalline specimen of KReO^. 

The multiple pulse sequences investigated are: 

i) Spin locking sequence (SLS) : (^/2 )^- (t-G^ ^ 

ii) PAPS : 

iii) WAHUHA sequence : (Tt/2) - t (t- (^/ 2 ) -'t-(n/2) -2 t- 

(V2)_y-T(Tt/2)^-'C-)^ 

From the spin locking sequence one can derive well known sequences 
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like CP (G = Tt , 0=0), OW (e = 7t/2, 0 = 90°) , CPMG (G = TT , 

0 = 90°) by choosing appropriate 0 and 0 values. 

II.B.l Experimental 

All the experiments discussed in this chapter have been 
carried out using the pulsed NQR spectrometer described earlier 
(see Section I. C. 3). The DC pulse sequences generated by the micro- 
processor system under the software control have been used to 
gate the r.f. The gating scheme for generating various phases 
in WAPTUHA sequence is shown in Fig. IV. 1. Same gating scheme 
can be used, with slight modification for generating PAPS sequence 
also. Delay line is employed to get arbitrary values of phase 
shifts . In this spectrometer despite the receiver gating the 
transmitter pulse feed-through could not be avoided. In all the 
experiments the responses of multiple pulse sequences have been 
acquired into the signal analyzer. These signals are processed, 
off-line, in the signal analyzer to remove the transmitter pulse 
feed-through. The processed signals can be either stored in a 
minifloppy disk or can be recorded on a plotter interfaced to the 
signal analyzer. 

Since in some of the experiments the multiple pulse res- 
ponses are. studied as a function of the flip angle of the sequ- 
ence pulses it is worth mentioning here the definition of the 
pulse flip angles in the case of the polycrystalline samples. 

We discuss only the case of axially symmetric field gradients 
(t)=0). In polycrystalline samples' a pulse which gives rise to 



(1) pp: Microprocessor 

(2) PSCQ, PSCT are broad band hybrid junctions 
from Mini-circuits, USA 

(3) Double balanced mixers (DBM) are from 
Hewlett-Packard USA (model 10514A) 
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Fig, IV*1 RF Gating scheme for the generation of WAHUHA sequence. 
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a maximTim FID is referred to as a 90° pulse. This 90° pulse 
flip angle is not equal to n /2 but it is 0.66 (for arbitrary 
spin I) . Similarly the flip angle of a 180° refocussing pulse 
is equal to 1.2271 [23]. In all our experiments we have followed 
the above-mentioned definitions in adjusting the pulse widths. 

IV.B.2 Results 

In what follows we present the results of experimental 
investigations of spin-locking, PAPS and WAHUHA sequences. 

IV.B.2 (i) Spin- locking sequence (SLS) : 

We first studied the influence of spin- locking sequence 
35 

on Cl in a single crystal of NaClO^. The crystal was oriented 
such that the r.f. field is perpendicular to the (1,0,0) plane. 
Under the influence of the spin- locking sequence the establish- 
ment of a quasi-equilibrium state of the magnetization could be 
observed. The magnetization has been observed to last for times 

much longer than the spin-spin relaxation time T 2 . In the case 
35 

of cl in a single crystal of NaCl02/ T 2 — 300 bsec and T^^, the 
spin-lattice relaxation time, is - 23 msec [24]. Results 
obtained with OW, CP, CPMG sequences are shown respectively in 
Fig. IV. 2 (a), IV. 2(b) and IV.2(c), while Fig. IV '.2(d) to IV.2(f) 
show the results obtained with arbitrary combinations of © and 0 
values in SLS sequence. From these figures it is evident that 
even in a single crystal the quasi-equilibrium state of the 
magnetization can be achieved with arbitrary 9, gf values for the 
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Pig. IV. 2 (a) . Response of Cl in a single crystal of NaClO^ 

to spin-locking sequence, (^/2) (t-0 / . -t) _ . 

Fig. IV. 2(f) "" 

The plots are given for various combination of 

0, 0 in SLS , In all the cases v= 29.9232 MHz, 

the resonance off-set. Aw --4 KHz and t= 180 

Msec. The relative signal amplitudes in all 

the figures are in the same scale. 




(s>!un ’CIJV) apn;ndujv pubis 











86-61 




o 

o 


00 

ay 

cr> 


(sviun *qJV) 2)pn;nclujV puSis 



200 


sequence pulses. The results of dependence of the decay of the 
quasi-equilibrium magnetization (QEM) on Q, and resonance 

off-set, Ao) / parameters will be presented in subsequent sections. 

In order to ascertain whether the spin-locking phenomenon 

is confined to only specific spin I values or it is general for 

1q7 185 

any spin, we have studied the influence of SLS on Re, Re 

(both with I = 5/2) in the polycrystalline sample of KReO^ . To 

our surprise we observed the spin-locking effect in these cases 

187 

also. In the case of Re the results obtained with OW, CP and 

/ 

CPMG sequences are shown respectively in Fig. IV. 3(a), IV. 3(b) 

and IV. 3(c), while Fig. IV. 3(d) to IV. 3(f) show the results with 

185 

arbitrary combinations of 0, 0" in the SLS. In the case of Re, 

the responses of OW, CP, and CPMG are shown in Fig. IV. 4 (a), 

IV. 4(b) and IV. 4(c). The results obtained with arbitrary 0, 0 

combinations in SLS are shown in Fig. IV. 4(d) -IV. 4(f). The spin- 

187 

spin relaxation time of Re in polycrystalline sample of KReO^ 
is "Sf 90 Asec. Now, from Fig. IV. 3 it can be easily seen that the 
decay of QEM lasts for times much longer than T 2 . These results 
and the results reported earlier in the literature [5,10,11,23] 
indicate that the spin- locking phenomenon is a general one and 
can be observed for any arbitrary spin I. 

IV.B.2(i)a Study of flip- angle '0* dependence: 

We have employed OW sequence for investigating the dpend- 
dence of the decay of QEM, on the flip-angle, '0', of the sequence 
pulses. The resonance off-set is fixed at -4 KHz and x is chosen 



Pig. IV. 3(a) 
to 

Pig. IV. 3(f) 


187 

Response of Re in polycrystalline speci- 
men of KReO^, to spin-locking sequence, 
(n/2)^-(x-0^^_^^j-T)^. The plots are given 
for various combinations of 0, 0" in the SLS . 
In all the cases v = 26.7817 MHz, the reso- 
nance off-set, is about -8 KHz. The 

relative signal amplitudes in all the 
figures are in the same scale. 
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Fig.IV-3(b) CP Sequence: 0=17, 0=0, T = 100psec 
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Fig. IV-3(c) , CPMG Sequence : 9 = 11, 0=90®, T = 100psec 
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Fig.IV-3(f) 0-W Sequence: 0=TT/2, 0=90®, T = 80/Jsec 



Pig. IV. 4(a) 
to 

Pig. IV.4(f) 
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Response of Re in polycrystalline sample 
of KReO. to spin-locking sequence, (7l/2) - 

* !X. 

Fig. 4(a) to 4(f) give the 
results for various combinations of ©, 0" in 
the SLS. In all the cases y= 28.3043 MHz, 
Resonance off-set. Aw, --8 KHz and t = 80 
i/sec. The relative signal amplitudes in 
all the figures are in the same scale. 
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to be 180 psec. For each *G‘ value we have computed the time 
constant for the decay of the QEM, '^2e' fitting the experi- 
mental data to an exponential function. The T 2 ^ values obtain- 

35 

ed for various values of © both in the case of Cl in single 

crystal as well as in polycrystalline sample of NaClO^ are given 
in Table IV. 1. 

Table IV, 1. Dependence of the Quasi-Equilibrium Magnetization 
Decay Time Constant* T 2 ^> on the Flip-Angle of the 
Sequence Pulses in the SLS 

Ao) = -4 KHz; t= 180 Psec 


e 

y 

- - - - - 

^2e 

(in msec) 

For Cl in 

single crystal 
of NaC 103 

' 35 

! For Cl in Powder 

1 sample of NaClO^ 

1 

- 1 

45° 

11.7 

17.6 

90° 

10.5 

16.8 

180° 

9.2 

14.0 


From the above table it is clear that as '©' is decreased 12^ 
increases indicating that the spin- locking effect is more effect- 
ive with smaller flip-angles of the sequence pulses. It is also 
evident from this table that T 2 g, for a given '©' in the case of 
polycrystalline sample is higher than that in for the single 

crystal. Similar kind of behaviour has been observed by Osokin, 
14 r T 

in the case of N in NaM02 L^J* 
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IV,B.2(i)b Study of 0 and T dependence: 

Until recently# there has been no rigorous experimental 
investigations of the dependence of the spin-loching phenomenon 
on 0 in SLS. Only recently Narsimha Reddy and Narasimhaji [11# 2l] 
have studied in detail the dependence of spin- locking effect on 
0 and found that in the case of nuclei with I = 3/2 in polycrys- 
talline specimens spin-locking effect could be observed with 
arbitrary 0 values in SLS. Their results contradict the theore- 
tical predictions due to Ainbinder and Furman [17] . 

187 

We have carried out the spin- locking experiments on Re 
185 

and Re in polycrystalline specimens of KReO^ and observed the 

establishment of QEM for different 0 values. From Figs. IV. 3(a)- 

IV. 3(f) and Figs. IV. 4 ( a) -IV. 4 (f) it can be seen that (for©=Tl/2# 

Tt) with 0=0, ^/2, TC spin-locking can be achieved in the case of 
187 185 

Re and Re also. These results indicate that the observa- 
tion of QEM with arbitrary 0 values in SLS# in the case of poly- 
crystalline samples is not restricted to spin I = 3/2 only but 
it is valid even for higher spins. 

One of the reasons for this observation of the spin-lock- 
ing phenomenon with arbitrary 0 values could be the polycrysta- 
lline nature of the sample where the crystallites are randonnely 
oriented with respect to each other and have their own efp 
principal axes systems which make a distribution of angles with 
respect to the applied r.f. field. 
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In order to verify the above reasoning we carried out 

35 

spin-locking experiments on Cl in a single crystal of NaClO^/ 
with different 0 values. In all these experiments we set 
-4 KHz, T = 180 psec, 9 = ^/2 or 7T . The results obtained with 
© = 45, 135, 225 are .shown respectively in Pigs. IV. 5 (a), 

IV. 5(b) and lV.5(c). The values obtained for 0 = 45, 135 

and 225 are set out in Table IV. 2. Now, from Pig. IV.2(a) - 
IV. 2(f) and Pig. IV. 5(a) -IV. 5(c) it is evident’ that even in the 
case of single crystal QEM has been established for arbitrary 
values of 0. 


Table IV. 2. T 2 ^ Values Obtained from the Study of SLS in the 

35 

case of Cl in Single Crystal of NaClO^ 

Aco = “4 KHz ; 0 = n/2 ; t = 180 psec 


1 

0 (in deg) [ 

-- . - - - - J ... 

T 2 g (msec) 

45 

10.6 

135 

11.9 

225 

13.2 










221 


Thus from this study it may be inferred that the spin-locking 
phenomenon can be observed both in polycrys talline samples as 
well as in single crystal with arbitrary 0 values. 

Marino and Klainer [S] who carried out the first spin- 

14 

locking experiment in NQR of N [I = IJ have investigated the 

— 5 

dependence of the T 2 g on t , and found that r . However, 

Osokin [s] who also studied the spin-locking experiments on 
in NaN 02 shown that the dependence of T 2 g on i is a complex 

one and that there is no simple relation which exists between 
them. Similar results have been obtained by Ainbinder et al. 

[lO] in the case of spin I = 3/2 nuclei in polycrys talline 
sample. 

In all our experiments we have observed spin-locking 
effect for all values of t^T~. However, as the values 

obtained by us are not very accurate owing to the transmitter 
pulse feed-through problem, we could not fit our experimental 
data to any simple relation between '^2q, * 

IV.B.2(i)c Study of resonance off-set, Aco , dependence: 

Although the spin-locking effect has been observed on reso- 
nance as well as off-resonance we have presented all the results 
so far with finite resonance off-set. This is on account of the 
fact that the spin-locked signal amplitude is much smaller at 
resonance than that in the case with finite resonance off-set. 
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Table IV. 3. Dependence of the Quasi -Equilibrium Magnetization 
Decay Time Constant '^2e Flip-Angle of the 

Sequence Pulses in the PAPS, Aco= 0; t = 180 qsec 


e 

-X 

1 

1 

T 2 g (msec) 

1 35 

1 For Cl in 

I crystal of 

1 

single) 
NaCl 03 1 

35 

For Cl in powder 

sample of NaClO^ 

45° 

15.2 


28.4 

90° 

13.5 


27.0 

180° 

11.6 


24.4 


Table IV. 4. Dependence of the Quasi-Equilibrium Magnetization 

35 

Decay Time Constant '^2e (°^'tsined from Cl in 
Powder Sample of SbCl^) on the Flip-Angle of the 
Sequence Pulses in the PAPS 

Aw =0 ; T = 150 jusec 

' - ' - ■' r - 

0 I T,. (msec) 

—X I ze 

> 

45 ° 11.2 

90° 10.0 

160° 8.4 
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35 

In the case of Cl in NaClO^ single crystal, CPMG sequence 
iwith T = 200 jjsec) while in the case of Re in polycrysta- 
lline sample of KReO^ OW sequence (with r = 80 usee) are employed 
to study the dependence of OEM on resonance off-set. A plot of 
Aq versus the amplitude of QEM obtained in these cases are given 
in Fig. IV. 6 and Pig. IV. 7, respectively. Similar kind of beha- 
viour of QEM has been theoretically predicted for the case of 
spin I = 1 [16], 

IV.B.2(ii) Phase alternated pulse sequence (PAPS): 

35 

We have investigated influence of PAPS on Cl in both 
single crystal and powder sample of NaClO^# and in polycrysta- 
lline sample of SbCl^. In this case all the experiments have 
been carried out on resonance. The decay constant of QEM as a 
function of '©' of the sequence pulses has been calculated and 
the results obtained in the case of NaClO^ are given in Table 
IV. 3 and those obtained in the case of SbCl^ are given in 
Table IV. 4. 

Typical recordings of responses of PAPS (with t = 160 usee, 

35 

AWss 0) obtained from Cl in a single crystal and in polycrys- 
talline sample of NaCl02 are shown respectively in Figs. IV.S(a) 

and IV. 8(b). In Figs. IV. 9 (a) and IV. 9(b) we have compared the 

35 

results of PAPS (with same x value and Aco= 0) obtained from Cl 
in single crystal and powder sample of NaClO^. 
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Fig. IV«9(a) Response of in a single crystal of NaCl03 to PAPS • 
9_ =(rr)_xi Ato=0, 't=240/jsec and i.>= 29-9232 MHz, 
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■35 

A typical recording of response of Cl in SbCl^ to PAPS 

(with t;= 120 Psec, Aw = 0) is given in Fig. IV.lO. The spin- 

35 

spin relaxation time, T 2 , of Cl in polycrystalline sample of 

35 

SbCl^ is about 0.5 msec. In the case of Cl in polycrystalline 
sample of NaClO^ T 2 order of a msec and -45 msec [25] 

From the above results again one can infer that under the 

action of PAPS the decay constant of QEM is much larger than the 

spin-spin relaxation time. From Table IV. 2 it may be noted that 

the smaller '©' values lead to longer values, also 12 ^ for a 

given ©,t; values, in the case of single crystal is less than 

that in the case of polycrystalline sample of NaClO^* Response 
187 

of Re in KReO^ to PAPS (with = 80 psec,A(ii = 0) is shown in 
Pig. IV. 11. 

IV.B.2(iii) WAHUA sequence: 

WAHUHA sequence is well known in high resolution NMR of 
solids. It suppresses the dipole-dipole interactions and leaves 
the chemical shift interaction active by selective averaging 
and lead to high resolution NMR spectrum. 

In the case of pure NQR, the only report on the study of 
WAHUHA sequence is that of Ermakov et al . [12] who investigated 

14 ■ 

the influence of this sequence on N (I = 1) in a single crystal 
of NaN02. To the best of our knowledge there has been no litera- 
ture report on the study of WAHUHA sequence in the case of spin 
I > 1. 
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We have experimentally investigated the action of WAHUHA 
35 

sequence on Cl in the polycrystalline samples of NaClO^ and 
187 

KCIO^ and on Re in KReO^. The results obtained in the case 

of NaClO^ and KCIO^ are shown respectively in Fig. IV. 12 and 

187 

Fig. IV. 13 while Fig. IV. 14 gives the results of Re. Prom 
these figures the lengthening of the decay of the magnetization 
is quite evident. This suggests that under the action of WAHUHA 
sequence some of the internal interactions which contribute to 
the line width are getting averaged out. It should be mentioned 
here that one of the important criterion for any multiple pulse 
sequence to be effective in averaging out the internal inter- 

ic 

actions is that the cycle time , should be much less than T 2 , 
i.e./ ^ "^2 L^J* However, because of the pulsewidth require- 

ments the cycle times employed in all the cases studied were more 

ic 

than or equal to T 2 . Despite the large cycle times the lengthen- 
ing of the decay of the magnetization could be observed. The 

decay constants of the magnetization as a function of x obtained 
35 

in the case of Cl in polycrystalline sample are set out in 
Table IV. 5. 

Table IV. 5. T- Values Obtained from a Study of WAHUHA 

35 

Sequence in the Case of Cl in Polycrystalline 
Sample of NaClO^ 

' " " ~ *" I " 

t(/Jsec) j T-^ (msec) 

I ^0 


100 

20.5 

120 

18.0 

150 

17.2 
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Fig. IV- 12 Response of in polycrystalline sample of NaCI 03 WAHUHA 
sequence : Acj=0, T=150psec, r) = 29-9232 MHz . 
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Fig. IV- 13 Response of in polycrystalline sample of KCIO 3 to WAHUHA 
sequence: AcJ =0, X = 300 p sec , 28-0623 MHz . 
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Fig.IV-14 Response of ^®’^Re in polycrystalline sample of KRe04 to WAHUHA 
sequence: A CJ =0 , T = 80jjsec, i) = 26*7817 MHz . 



From the above table it is clear that the smaller x values yield 

larger decay time constants. Similar trend has been observed in 
35 187 

the case of Cl in KCIO 2 and Re in KReO^. 

IV. C DOUBLE RESONANCE UNDER MULTIPLE- PULSE SPIN-LOCKING 

As discussed in Chapter 1/ in NQR the direct detection of 
quadrupolar nuclei having low-natural abundance and/or low qua- 
drupole coupling constant (such nuclei termed as rare nuclei) 
is rather difficult. Double resonance methods offer better sensi- 
tivity in detecting such rare nuclei. There have been several 
double resonance techniques available in the literature (see 
Chapter I, Section I.C.2 ). 

Marino and Klainer [5] who carried out the first spin- 
locking experiment in NQR have proposed that, in principle, the 
sensitivity of double resonance scheme of Emshwiller, Hahn and 

Kaplan [26] could be enhanced under the spin- locking condition. 

23 35 

Encouraged by this proposition, we have carried out Na- Cl 
(in polycrystalline sample of NaCl 02 ) double resonance experi- 
ment under the spin-locked condition. In what follows we present 
the experimental details and the results of the double resonance 
experiment. 

IV.C.l Experimental 

For this double resonance experiment we have chosen CPMG 
sequence. The DC pulse sequence in this case has also been 
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generated by the microprocessor. First the preparatory (Tr/2) 

pulse and the sequence of tt pulses were generated respectively 

at BNC-1 and BNC-2 of the microprocessor system. The Ti pulse 

sequence was passed through a current driver. One of the two 

outputs available at the current driver was used to gate the r.f. 

23 

in Na channel while the other output was used to gate the 90° 

35 

phase shifted r.f. in cl channel and which is subsequently 

35 

used to generate the CPMG sequence in Cl channel. The experi- 

35 

mental set-up in Cl channel is same as described earlier in 

23 

Chapter I. The gated r.f. pulse sequence in the Na channel 
has been amplified by an r.f. amplifier (Model GR 1233-A, 20 Hz 
to 3 MHz from General Radio, USA) which gives a maximum power 
of 15 watts, and coupled to the double resonance probe. Thus, 
basically this experiment is an extension of the Spin Echo Double 
Resonance (SEDOR) experiment [26]. 

IV. C. 2 Results 

Although we have tried various combinations of 9, 0 in the 
spin-locking sequence, CPMG sequence was found to be more sensi- 
tive in this experiment. Since, on resonance, the spin-locked 

signal intensity is rather small we have set an off-set of 

35 23 

-3.5 KHz in Cl channel. Upon irradiation in Na channel at 

23 

the resonance frequency of Na, the amplitude of the spin-echo 
i 35 

i train in Cl channel has been found to be reduced considerably. 

The super imposition of the two echo trains obtained in the presence 
and absence of the double resonance irradiation is shovm in 


?ig. IV. 15. 
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35 

A plot of the Cl spin-echo amplitude, monitored after 

- 5 T 2 time, versus the frequency of "^Na channel is given in 

Fig. IV.16. From this figure it is evident that the resonance 
23 

frequency of Na is - 397.5 KHz at 295 K which is in good agree- 
ment with that obtained earlier by other double resonance methods 
in the literature. 

From the experimental results of multiple-pulse sequences 
presented in Sections IV. B and IV. C the following conclusions can 
be drawn. 

i) Spin-locking phenomenon can be observed with arbitrary 
values of 0 in SLS, in polycrystalline samples as well as 
in single crystals. These results cannot be explained by 
any of the theories currently available in the literature . 
on the spin-locking phenomenon in NQR. A more detailed 
theoretical analysis of spin- locking phenomenon in NQR of 
polycrystalline samples as well as in single crystals is 
therefore warranted. 

ii) It appears that spin-locking effect is not restricted to 
spin I = 1 or 3/2 nuclei only but it can be observed 

in the case of nuclei with any arbitrary spin I. 

iii) Spin-locked signals have higher intensity in the presence 
of a finite resonance off-set. 

iv) Enormous improvement in the S/N can be achieved by co-adding 
all the echoes in the spin locked spin echo train. 




SIGNAL AMPLITUDE (Arb. Units) 




v) Under the Influence of SLS and PAPS the decay time cons- 
tant of QEM increases as the flip-angle of the sequence 
pulses is decreased indicating that spin-locking is more 
effective with short flip angle pulses. 

vi) Under the action of PAPS the lengthening of the transverse 
magnetization could be observed for arbitrary 0 values of 
the sequence pulses for nuclear spins with general spin I. 

vii) WAHUHA sequence leads to the lengthening of the transverse 
magnetization in the case of NQR of spin I = 3/2 and 5/2 
also 

viii ) As pointed out in the beginning of this chapter a sensi- 
tive double resonance experiment under spin-locking condi- 
tion can be carried out to detect rare quadrupolar nuclei. 

To summarize all the aforesaid multiple-pulse sequences lead to 
the lengthening of the decay of the transverse magnetization 
implying a narrowed NQR line in the frequency domain. 

In order to understand the kind of internal interactions 
that are averaged out under a particular sequence one needs a 
theoretical analysis of the influence these sequences . 

We have carried out theoretical analysis of some of the 
multiple pulse sequences, using the tensor operator formalism, 
the details of which will be presented in the next section. 



IV. D THEORETICAL ANALYSIS OF A FEW MULTIPLE PULSE EXPERIMENTS 
IN NQR OF SPIN I = 3/2 NUCLEI 

In this section we present our theoretical analysis of the 
influence of multiple-pulse sequences such as PAPS, WAHUHA [13] 
and Mansf ield-Rhim-Elleman-Vaugham (MREV-8) [27] on spin 3/2 
nuclei. We have employed the Average Hamiltonian Theory [l] and 
tensor operator formalism [28-32] in this investigation. 

In NQR spectroscopy the main internal interactions which 
lead to line broadening are [6,33], 

i) inhomogeneities in the efg's 

ii) torsional oscillations of molecules 

iii) heteronuclear dipole-dipole interactions 
and iv) homo-nuclear dipole-dipole interactions. 

In the present analysis we consider only the first three of the 
above four interactions . The Hamiltonians corresponding to these 
three internal interactions are listed below. 


IV.D.l Hamiltonians of Internal Interactions 

We consider a nuclear spin system with I = 3/2 in a non- 
axial efg. 


rv.D.l(i) Hamiltonian of inhomogeneous efg: 

The Hamiltonian which defines the internal interaction due 
to inhomogeneous efg can be written as 




inhom 


Aco 

v/r 


1^4 - 


+ 7](lx - 




(IV.I) 
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where is the frequency deviation from an average value 

2 . 

® qQ/2 . In terms of tensor operators this Hamiltonian can 
be written as 


^inhom = ^ [^o = •• 

IV.D.l(ii) Hamiltonian corresponding to torsional oscillations; 

Let us assxime that the molecule containing the quadrupolar 
nucleus undergoes torsional motions about the three principal 
axes of the field gradient. The axes fixed in space are denoted 
by primes and the axes fixed in the molecule without primes . Then 
for small rotations and ©^ about the three axes X, Y, 2 it 

may be shown [33] that the primed and unprimed components of the 
field gradient tensor are related by the following equations: 


^X'X' 


V 

XX 


Vyiy« 

= ®z ^XX^ 

(1- 

4-^z > ''yy + ®x V 

^Z'Z' 

= ®y' ^xx + 

®x^ 

V^Y + " ®Y ^ ^ZZ 

V^.Y. 

= -®z ^xx 

^®x 

®Y^2> ^YY-®X®y'^ZZ 

Vy 1 2 • 

= ®Y®Z^XX " 

®x^ 

YY ” ^®Y®Z "®X^ ^22 


^Z'X* “®Y^XX “ ®X®Z^YY ^®X®Z ®Y^ ^ZZ 

-eq/2 H+n) 


where = eq/2 (tj-I); V, 


.. (IV. 3) 


YY 


and ^22 = eq. 
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The quadrupolar Hamiltonian of a nucleus in a molecule which is 
undergoing torsional oscillations can be written as 

^ ° 61(21-1) [I-’'' -I] •• 

where V* is the efg tensor. 

The above Eqn. can be rewritten as 




tors 


. (IV. 5) 


where = — [T^ + rj/Js T^Cs)] 

'/S 


.. (IV. 6) 


is the pure quadrupolar Hamiltonian and Hamiltonian 

which defines the interaction between the spin system and torsi- 
onal oscillations. We are mainly interested in "K^Qj-g* Conside- 
ring only the secular part of 'I^^QPg i*e./ it can be shown 

that 


(4 - “v^v ^4 " 4^ ^ “z^z ^4 ■ 4^ 

tors X X ^ y Y Y ^ 2: ^x z z ^y 


. (IV. 6) 




1 (3+77) _ 1 ( 77 - 3 ) 

X ~ 9 2 ' ‘^y ” 9 2 ' 9 
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e = £ = ©f, and = 0^ 

X X y y z z 


with (^Q == e qQ/2 


In terms of tensor operators can be written as 



"^tors ^2 

\l6 

where = -^ (-co £ +co, £ ) 

1 2 X X Y 


b2 - — (-cOx^x'Vy 
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.. (iy.6a) 

.. (IV. 6b) 

.. (IV. 6c) 


IV.D.l(iii) Hamiltonian of heteronuclear dipole-dipole 
interactions : 

The dipolar interaction between the two magnetic moments 
= Y Ii / "^2 ~ "^2^ ^2 given by the Hamiltonian [34] 


(I, .!„) 3(1, .r) do.r) 


^dip 1'’^ 2^ t 3 “ 5 i 


.. (IV. 7) 


where y gyromagnetic ratios of the nuclear spins 
The above equation can also be written as [34]: 


^dip “ "^2 ^2 ' ^12 


where H ^2 called the local field produced by spin 1 at the 
site of spin 2, 

Thus, the heteronuclear dipole-dipole interaction can be 
treated as the interaction of the resonant nuclear spin with a 
local magnetic field produced by the non-resonant nuclear spin 
at the site of the resonant spin. Thus, 

Z 

= - 7 S j ^ 

m=x m m 


% 


het 


. . (IV. 8) 



where we have set "h = 1, and is the component of the local 
field in 'm* direction. 

In terms of tensor operators written as 

^het = ->'[-Tj(a) + tJ(s) H^] .. (iv.9) 


As mentioned in Section IV. A, under the action of a multi- 
ple-pulse sequence an internal interaction is said to be avera- 
ged out if the following conditions are satisfied [20]: 


^^li[=0or[P,,Jc'«] = 0 


.. (IV. 10) 


where is the average Hamiltonian, of i"*-^ order. 


of an 


internal interaction, and is the density matrix prepared by 
the preparatory pulse of the sequence. 


Since in this section we are mainly interested in the status 
of the various internal interactions under the action of the multi- 
ple-pulse sequences we first calculate the average Hamiltonians 
(we restrict ourselves only to the zeroth order average Hamiltonian) 
of the internal interactions and then evaluate the density matrix 
following the preparatory pulse and subsequently look for the 
commutation relationship of the form given in Eqn. IV. 10. Depend- 
ing upon the result of the commutator between and 
infer whether a particular internal interaction under a given 
sequence averaged out or not. 



In order to present a coherent discussion on the above 
mentioned aspects in the next sub-section we give a brief outline 
of average Hamiltonian theory and in the subsequent sections we 
discuss the evaluation of zeroth order average Hamiltonians of 
the internal interactions, under some of the multiple-pulse 
sequences, and their . commutation with the 

IV. D. 2 Multiple-Pulse NQR and the Average Hamiltonian 

Average Hamiltonian theory'’ is well documented [!}. In this 
sub-section we closely follow the treatment outlined by Haeberlen 
[l]. Prom here onwards in all the calculations in this chapter 
we set IFi = 1. 

Let us consider a nuclear spin system governed by a total 
Hamiltonian ?^(t) , given by 

<K(t) .. (IV. 11) 

where 'Kq is the quadrupolar Hamiltonian, ^^(t) is the Hamil- 
tonian for the interaction between the r.f. field and the 
nuclear spins, and represents the Hamiltonian correspond- 

ing to the internal interactions. 

In the quadrupolar interaction representation (QIR) 


defined by 
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The total Hamiltonian is given by 




'K = %it) +^. 
1 ' ' 1 


int 


In QIR, the equation of motion of the density matrix i 

by 


dp 

dt 




Now let us make a transformation 


P(t) = U^(t) P(t) U^(t)' 


v- ^ ^ 

where (t) = T exp (-i f it. (t) dt ) 

^ o ^ 

T is the Dyson time-ordering operator 


In this new frame the motion of P is governed by 


fl = [Jfint' 




The formal solution of Eqn. (IV. 17) is given by 


P(t) = u.^^(t) 2(0) n-J^(t) 


where U.^^(t) = T exp (-i J ^^^^(t')dt') 


Then from Eqns . IV. 15 and IV. 19/ we have 


. (IV. 13) 

is given 

. (IV. 14) 

. (IV. 15) 

. (IV. 16) 

.. (IV. 17) 
. . (IV. 18) 

. . (IV. 19) 
. . (IV. 20) 
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P(t) = U^(t) U^^^(t) P(0) U^^(t) .. (IV. 21) 

Since the average Hamiltonian theory is based on the cyclicity 
of the r.f. pulse interactions , ( t) , we now briefly discuss the 

cyclicity and periodicity of the multiple pulse sequences . 

A pulse sequence is said to be cyclic if the propagator 
operator Uj^(t^) defined by the sequence at a cycle time t^ is 
unity [35] , i.e., 

Ui(tc) = T exp (-i j^^^(t) dt) =1 .. (IV. 22) 

o 


Using Magnus expansion [36] the above equation can be rewritten 
as 

Ui(tc) = exp { -it[ + ...]} .. (IV. 23) 


where 




K^(t) dt 


(IV. 24) 


= - 2^ J 'it J dt' [ '?^^ ( t) , 'K^( t* ) ] . . (IV. 25) 

^c o o 

From Eqn. IV. 23 it is evident that a sufficient condition 
for designing a cylic sequence is by the requirement that 

= 0 for all i .. (IV. 26) 


A sequence is said to be periodic, i.e., the sequence is repeated 
periodically with a cycle time t^ 
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if %^(t) (t + Nt^) .. (IV. 27) 

For a sequence both cyclic and periodic 

Ui (Nt^) =1 (N = 0, 1, 2, ...) .. (IV. 28) 


Now let us assume that our observation Windows in the 
experiment are set at the integral multiples af the cycle time 
i.e., t = Nt^. Then from Eqn. IV. 21 


5(t) = ?CKtp = P(0) D-l<Nt^) 




. . (IV. 29) 


Since Uj^(t) = U^(Nt^) = 1 for a cyclic sequence 


Since U 2 (t) is periodic i.e.,U^(t) =Uj|^(t+t^) and 
determined through Eqn. IV. 18 


'^int 

hence = [Umt'Vd” ' 

Consequently, 

p(Ntc) = [Ui^t(t^)f P(o) 


. . (IV. 30) 

. . (IV. 31)' 

. . (IV. 32) 


Thus, this equation implies that, in order to describe the state 
of the system, under the action of a cyclic multiple-pulse 



sequence / at any integer multiple of the cycle time t^ it is 
sufficient to know its short-time evolution over one cycle • 

We can now make use of Magnus expansion [36] to evaluate 




with 

^iO) 1 

“^int t 

c 



"^int " 2t i i d0' R. (0') ] 

CO o int ** 


(IV. 33) 


(IV. 34) 


(IV. 35) 


and higher order terms . 

^int referred to as the zeroth-order average Hamiltonian 
while is the first-order average Hamiltonian. 

In order that this series is strongly convergent 


int 


- t w. , 1 

c int 


. (IV. 36) 
•( 0 ) 


In our theoretical analysis we are mainly concerned with 


IV. D. 3 Theoretical Analysis 

In this section we first present the explicit expressions 
for relevant Hamiltonians in QIR and then proceed on to evaluate 
zeroth order average Hamiltonians of the internal interactions. 

Consider a non-interacting nuclear spin system which is 
characterized by a total Hamiltonian, ^(t) (see Eqn. IV. 11) 
given by; 
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•Kct) = +'K,(t) ^ 

Q 1 ' ' int 

where %q is the pure quadrupolar Hamiltonian defined in 
Egn. IV. 5. 

If we assume that the r.f. pulses are applied along x-axis 

of the crystal with a frequency co and amplitude H. then we 

q 1 

have for the r.f. Hamiltonian, 

cos(cOgt + 9 ) T^(a) .. ClV.37) 

with 0 )^ =7 H^; T is the gyromagnetic ratio of the nucleus and 
<p is the phase of the r.f. Considering specific cases of 9 
values we have 

for 9 = 0 3(^(t) = 2 ( 0 ^ cos(( 0 gt) T^(a) 
for 9 = 90° I^^Ct) = -2(j)^ sin(co^t) T^(a) 
for 9 = 180° ^^(t) = - 20 )^ cos(a)gt) T^^Ca) 
and for 9 = 270° '^{^(t) = 2(o^ sin(w^t) T^(a) 

^int '^inhom '^tors ’'‘^het 

where ^i^hom' ^tors ^het given respectively by Eqns . 

IV. 2, IV. 6 (a) and IV. 9. 

As has been mentioned in previous chapters calculations v/lll 
be simplified if we work in a representation in which the quadru- 
polar Hamiltonian is diagonal. In the 'diagonal representation' 
defined by the transformation operator. 


. . (IV. 38) 

.. (IV. 39) 


= exp ( 2 /v^ X T 2 (a)) 


. . (IV. 40) 



the total Hamiltonian is given by (see Table 11.12), 
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where cd 

q 

with ^ 


0 ) 


Q 


(1 


.. (IV. 41) 

. . (IV. 41a) 
.. (IV.41.b) 


For 9=0 ^^(t) = 203^ cos (w^t) [a^T^ ( a) + a 2 T^(a) + a 2 T 2 (a)] 

.. (IV. 42) 

For 9 = 90° ^^(t) = -2(0^ sin(w t)[a^T 2 (a) + a 2 Tj(a) + a 2 T 2 (a)] 

.. (IV. 43) 

^^(t) for 9~ 180° and 9~ 270° can be obtained, respect- 
ively from Eqns . IV. 42 and IV, 43 by changing the sign of 
Hamiltonians . 


% 


where 


•K 




a 

int 

: -t/ 

inhom 

d 

inhom 

Aw 

= -.2 T' 

v/T 

d 

tors 

= A. ■ 

1 r^O 

d 

het 

11 

1 1 

1 


+iH. 


' + H 




tors “^het 


.. (IV. 44) 


.. (IV. 44a) 


.. (IV. 44b) 


r [-H^ [a^Tj(a) + a^Tj(a) + a3T3(a)j 


{a4Tj(s) + a^T^(s) H 

fSf ^ 

z [®7'ro + ^S’^o + 




(IV.44c) 



where = 1 /^ (b^ + ’l/s[3 
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■^2 “ (-‘n/\f3 + b2) 


a, = 


^4 ^7 “ 


SvITsf 


[3(1 -f) + 577 ] 


^3 ~ 3f ^ f +h- 


1 ) 


= 


= 


3 ^^5y 

_1 


[3(1 -f)-577] 


"8 


a„ = - 


3^ [l-f+'l ] 


(i-y) 

n/y 


3^/^o f 
2\f2 

3n/3' 


. . (lV.44d) 


b^, b 2 and f are given respectively by Eqns. IV. 6(b), IV. 6(c) 
and IV. 41(b). Having obtained the total Hamiltonian of the 
system in the diagonal representation, we now transform it into 
the QIR defined by 


U 


q 



. . (IV. 45) 


In the QIR, the total Hamiltonian can be written as (see 
Table III . 1 ) i 


U 


q 




d 




t'K? 


int 


. . (IV. 46) 
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where 


for 9 = 


0 X 1 = 9 = 90° = co^S 


= 180° = -w^R, 9 = 270° 


■■ 1 


^int '^inhom ■*■ "^tors ^het 


Ao), 


q 2 


with 'K • u = T 

^inhom jg- 


-oo^s 


.. (IV. 47) 


.. (IV. 47a) 


(IV. 47b) 


(IV. 48) 


.. (IV. 48a) 


= A.T 

outers 1 o 


.. (IV. 48b) 


*K * 

het 


[-H^ ij (f+l-'DLl T^Ca) ^ Tj(a)]+ ^^(a)) 

+aHy t| (f+l+DCl Ti(s);^ Tj(s)] (-f+l+n) TjCs)) 

.. (IV. 48c) 


U+ii 1 _ 4(1-)’) 3 , 

* i SJ. Jo 3^,j. JeO 


In arriving at the Eqns . IV. 48(a), IV. 48(b) and IV. 48(c) we 

have made use of the evolution of tensor operators under the 

^ a 2 

action of quadrupolar Hamiltonian — ^ T (see Table III.l) and 

1/6 ° 

dropped the high frequency terms as nonsecular terms. 

^ d 

The evolution of tensor operators under the action of 
(for 9 = 0) can be obtained from Table II.3 by substituting £. by 



P - ^ ^ 1 "t . From the same table, the evolution of tensor 

x/J ? ^ ^ 

operators under (for f = 180*^) can be obtained by substituting 

£ by - 3 . The evolution of tensor operators under (for 

^ = 90°) are set out in Table IV. 6 . Prom this table the evolu- 
tion of the operators under the influence of (for 9 = 270°) 

can be obtained by replacing 3 by -3 . 


Having obtained the Hamiltonians corresponding to, the 
internal interactions and r.f. interaction with the spins, in 
the QIR, V 7 e are nov 7 ready to discuss the influence of the pulse 
sequences on the internal interactions. 


We shall now examine whether the pulse sequences under 
investigation are cyclic. We begin with PAPS, i.e.. 


(V2)^-[T-(n)^_2T-(Tl)_^-T] .. (IV. 49) 

Considering only the zeroth order term from Eqn. (IV. 24) we 
see that 

= ~ Jt^(t) dt = (R-R) =0 .. (IV. 50) 

^c o 

where R = [| T^a) + ~ T?(a)] .. (IV.51) 

3^ '-5^1 ^ -• 

Since in PAPS only the operator R is involved in it can be 

easily shown that all higher order terms in the magnus expansion 
— '(i) 

i.e., vanish. Thus, we can establish that PAPS is a cyclic 

jequence. Proceeding in a similar manner it can be shown that 



Table IV, 6. Evolution of Tensor Operators T^(a,s) Under the 
Influence of for Spin I = 3/2 


UTt ( a)u“^ 


UtJ’ (s)u"^ 


UT?(a)U ^ 

/vt X 


UT? (s)u”^ 


UTqU“^ 


3 -1 
ut;:u 


2 -1 
UT U 
^ o 


UT^ (a)U 


2/ 

UT2(s)U 


UT^(a)u”^ 


jj^Ca) [1 +X(cos£-l) ] + T^ (a) (cos£-l) 

\/T 5 


■|(cos£-l) T^(a) 

J ^ J 


— sin£ - — sinE T?(s) 
^ o ^2 


Ti (s) [1 (cosE- 1 ) - — ^ T? (s) (cosE- 1 ) --^(cosE-l) T?(s) 
^1 L 5 ^1 3 ~3 


— i sinE T^(a) 


= t; 


= T 


cosE T?(a) - sinE T^ - sin£ T^ 

b ^o <-0 


T^ (s) 

^[1 +i(cosE-l) +^— T^(cosE-l) +-^ sin£ T?(a) 
= 3\f5 ^ ~ 

+ (cosE-1)[4 tJ + — tJ] + — sinE T^Ca) 
,o L5 ^o -oJ ^ 

T^ cosE + sinE [T^(a) T^(a)] 


T 2 (a) COSE + sinE [-— tJ(s) - tJ(s) - T^Cs)] 

^ Jib \fl ^ 


= T^(s) cosE + sinE [— T^(a) ^ T^(a) - T^(a)] 

^2 ■- 5 n/To Ni 3 


T?(a) + (cosE-l)[T 5 (a) + - - - - - tJ'-' - 

--1 ^--1 10 rn^^^r'l 


- -- - T^Ca) + ^ T^Ca)] 

lOvfTs 2V5 

3 m2 

T -t- — — - -x ^ ' 

v/To 2'JlO 


+ SinE [- -^ t; + --|r- T2 (s)] 


. . .contd . 


Table IV. 6 (contd.) 


UT^ (s)u“^ 

rJ X 


tJ(s) + (cosP-l)[j| tJ(s) 


+ sin£ [-j- Tlia)] 

2 \/Tb 


UT2(a)u~^ 


T^(a) 


^ tJ(s) + — t!(s)] 

10\/r5 2\/5 


UT 5 (s)u"^ = T^(s) 


UTo(a)u“^ 

r-' -J 


UT2(s)u“^ 


= T^Ca) [ (cos£-l) ][i T^Ca) - 

+ ^ sin£ [t^(s)] 

2VZ 

= T3(s)[ 1 4- "I (cos^-l)] + [ 

+ sinjS 

2\/2 


if 3i<s) + ^sj(a)] 


2'^S ^ 


tJ(s) - T^(s) ] (COS3-1) 


where 


« [- a Ths)] 

3f >f? ^ 


= exp(-i iut) 


(3+r)) 
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WAHUHA, i.e.. 


('^/2)^-z(.x-in/2) _^-T - /2 ) 2T - (TT /2 ) _^_x _ ( n/2 ) ^-x - ) ^ 


and MREV-8, i.e.. 


.. (IV. 52) 


(Tr/2) -x(T:-(n/ 2 ) _x-(V2) -2 t- ( V2) -x-(n/ 2 ) -x- 

Jv y 

(71 /2 ) - (71/2 ) 2x - (71 /2 ) y-T - (Tr/2 ) - ) ^ 

.. (IV. 53) 

Pulse sequences are cyclic at least to zeroth order. 


IV.D.3(i) Zeroth-order average Hamiltonians of internal inter- 
actions under the influence of PAPS and their commu- 
tation with case of single crystals: 

In this sub-section we first evaluate the zeroth-order 
average Hamiltonians of internal interactions and seek their 
commutation v/ith P^, the initial density matrix prepared by the 
preparatory pulse. 


The Hamiltonians corresponding to the internal interactions 
are given by Eqns. IV.44(a) -IV. 44(c). Making use of Eqn. iy.34 the 

zeroth-order average Hamiltonian for efg inhomogeneity can be shown 
to be , 


7i 


(0) 


1 


inhom 


r ^ 


inhom 


(x) dx 


A(0 


'f6.2 


^ [t^(1 +COS/3 2 ) +i sin3 


(IV. 54) 
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Similarly, 

( 0 ) _ 


^tors “ [^o (1 +cos) 32) + i sin /3 2 tJ(s)] .. (IV. 55 ) 


and 


?^het - [(l-r])cos 32 + fl [-^(^ T^(a) - Tj(a))} 


f = " ■ ' ^/^5 

jI [-(1-T)) cos ^2 3 ^3^^^ 

\f 2 3 

+ i sin ,8 2 3^ (n-l) T^Ca)} 

+ iHy{[(l+77) +COS 8 2fKj tJ(s)-^ tJcs))] 

+ 3 ^ [(l+7))-cos 02?] ^ 2 ( 3 ) 

-i sin 8 2 ^ T2 (s) J 


•*■ i-sp 


y(cos 82j’ -1) 


3^0 

+ isin^ 2 [-^ T^(s) + 


v/s „ 1 , . . 2 _3 


3\[5 


where 8- 


2 (o). 


(3 +T)) 
N/ 3 f 


t,„) 


vu' 


Ti(s)] 3 ] 


(IV. 56 ) 


.. (IV. 52 ) 


In obtaining Eqns . IV. 54 , IV. 55 , IV. 56 we have made use of the 
Table II . 3 which gives the evolution of tensor operators under 
the action of "Xf (for <p = G) when £ is replaced by 8 = a> t. . 

X ui 

<3 

In obtaining the evolution of the operators under ^ (for 9= 180 ) 
we substituted £ by -8 in Table II. 3 .. 
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The initial density matrix prepared by the preparatory 
pulse is given by 

^ 5'o i3 1 + i sin^^ T^(s) .. (IV. 58) 

where/3 = o) t .. (IV. 59) 

\/3 ? 


Now/ from Eqns . IV. 54, IV. 55, IV. 56 and IV. 58 and commuta- 
tion relationships of tensor operators for spin I = 3/2 (see [29]) 
it can be easily shown that 


[ Kl. ] = ° 

[ 5 . ] = 0 

[ 2. ] = 0 ■ .. (IV. 60) 

when 23 ^ = $ 2 • 


Thus, from the above discussion it is clear that under the 
influence of PAPS the internal interactions defined by ^ 
tors ' ^ het averaged out when the pulse widths of the 

cycle pulses are twice to that of the preparatory pulse. This 
indicates that the PAPS sequence is effecting a kind of general 
averaging leading to the lengthening of the decay of the magne- 
tization. This is in agreement with the experimental results of 
PAPS on single crystal presented in Section IV.B. 
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IV.D.3(ii) Influence of PAPS on spin I = 3/2 nuclei in poly- 
crystalline samples: 

In the proceeding sub-section we dealt with the action of 
PAPS on spin I = 3/2 nuclei in single crystals. In this section 
we examine the influence of PAPS on spin I = 3/2 nuclei in poly- 
crystalline samples.- In the case of polycrystalline samples the 
crystallites are oriented randomly and hence each crystallite 
will make a different angle with respect to the applied r.f. 
field. Hence, the r.f. Hamiltonian for this situation should be 
written as 

^^(t) = -2(0^ cos(Qgt + cp) [l^P^ +IyP2 +12^3] •• (IV. 61 ) 

1 

where P^ = sinG cos^, P 2 = isinG sin/, P^ = cosG, G and / 

are the polar and azimuthal angles of the applied r.f. 
field vector and 9 is the phase of r.f. 

In terms of tensor operators the above equation can be written as 

“K^Ct) = -2(0^ cos( 0 )gt+(p) [-T^Ca) P^ + T^(s) P 2 t F 3 ] 

.. (IV. 62) 

For 9=0 = -2w^ cos(Wgt) [-t]^ ( a) P^ + (s) P 2 + Tq P^J 

.. (IV. 63) 

For 9 = 180° = 2co^ cos(Wgt) [-Tj;(a)P^ +tJ(s)P 2 +J 0 P 3 ] 

.. (IV. 64) 

Transforming ^^(t) first into 'diagonal representation' and 
then into QIR, we have 
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.. (IV. 65) 


for 9=0 Cf co^ [p^( 3 +;j )(3 ^l(a) + 


2 rr,3, 

T^Ca)) 


JTs 


- P 2 ( 3 - 7 )) (| tJ(s) + T^(s)) 

^ JTs 

. P 3 


(IV. 66 ) 


^ Cj[ 

and (for 9 = 180*^) can be obtained from Eqn. IV. 66 by just 

changing the sign of the Hamiltonian, i.e.. 


(Xf) 


f> , o' -‘5.?) 

9 =180 9 = 0 


(IV. 67 ) 


The evolution of tensor operators under the influence of JC. ^ 

(for = 0) is set out in Table III. 3. From Table III. 3 the 
evolution of the operators under (for 9 =180°) can be obtain- 

ed by replacing E hy -E . 

Before evaluating the zeroth order average Hamiltonian, we 
consider the question of cyclicity of PAPS in the case of poly- 
crystalline samples. Again considering only the zeroth order 
term,- from Eqn. IV. 24 and using Table I I I. 3, we see that 


1/t^ J'' ^°(t) dt = W^/4 [T-T] = 0 


^c 


(IV. 68 ) 


where T = 


3 1 


rf hi (3+>l)(f TjCa) + ^ T,(a)) 
- P, (3-7)) (| Tj(s) + Tj(s)) 


P 


3 = 


(IV. 68a) 
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Thus, even in the case of polycrystalline samples PAPS is 
cyclic atleast to the zeroth order. 


Zeroth-order average Hamiltonians of internal Interactions and 

their commutation with p. : 

1 

In the case of • polycirystalline samples we consider the 
internal interactions due to efg inhomogeneities and torsional 
oscillations only. Making use of Eqn. IV. 34 and Table III. 3 the 
zeroth order average Hamiltonian of inhomogeneous efg can be shown 
to be 


— Aw 

^(0) ^ 
inhom 




{P^(3+r)) T^(s) - P 2 ( 3 - 7 )) T^Ca) +P 3 2 t 7 T^Ca)}] 


.. (IV. 69) 


Similarly, 




( 0 ) 

tors 


A- ^ isinS- ^ 

= ^ [t^(1 +COS £ 2 ) + {P.(3+77) T-(s) 


-P2(3-n) T^Ca) + P3 2T}T2(a)}] 

.. (IV. 70) 

The density matrix prepared by the preparatory pulse is given by 
„ isin £ 


P(tuj) = T^ cos £ T + 


jBfx 


- [P3(3 +n) ^3(3) -P2(3 -V) Ti(a) 


+ 2 r) P 3 T 2 (a)] 


(IV. 71) 
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where {sin^e [(9+ri^) - 6n cos2iZf] + 4Ti^ cos^e]]^^^ 

(3+7) ) 

.. (IV. 72) 


^2 = 


From Eqns . IV*69/ IV. 70 and IV. 71 and commutation relationships 
of tensor operators for spin I = 3/2 (see [29^) it can be. easily 
shown that 




. . (IV. 73) 


Thus, even in the case of polycrystalline samples under the 
action of PAPS the internal interactions due to efg inhomogeneit- 
ies and torsional oscillations will be averaged out, when the 
width of the cycle pulses is twice to that of the preparatory 
pulse. Similar kind of behaviour has been observed by Osokin [6[[ 
in the case of spin 1=1. Thus the experimentally observed, 
lengthening of the decay of the magnetization can be explained 
as, partially, due to the averaging effect of ^inhom' ^tors* 
However, experimentally we have observed the lengthening of the 
decay of the magnetization, for arbitrary values of the width of 
the cycle pulses. This observation cannot be explained by the 
theoretical analysis presented above. 

IV.D. 3(III) Influence of WAHUHA and MREV-8 sequences on 
spin I = 3/2 nuclei in single crystals; 

In this section we consider the influence of WAHUH/i and 
MREV-8 sequences on internal interactions due to efg 
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inhomogeneities and torsional oscillations only. 

(a) WAHUHA sequence: 

Using Eqn. IV. 34 and Tables II. 3 and IV. 6 the zeroth— order 
average Hamiltonians of can be shown to be 


inhom 


Aco_ 


[ (l+cos£ +cos 2£) - 2 isint T?(s) 

J6.3 


— ^ sin2£ (| Tha) + T^(a))l 

v2 ^ w-*- jYs 

.. (IV. 74) 


^ tors ~ “T +cos£ + cos^E ) -2 isin£ T^(s) 

p sin2£ [| T^Ca) + T^(a)]] .. (IV. 75) 

■v/2 ^ vi5 


The initial density matrix prepared by the preparatory pulse is 
given by 


P ( t,.) = T^ cosjS + isir^ T?(s) 

^ r«» O ^ JL 


(IV. 76) 


Now/ from Eqns . IV. 74, IV. 75 and IV. 76 and commutation relation- 
ships of tensor operators for spin I = 3/2 (see [29]) it can be 
shown that 


C 

[ p. 


(0) 

inhom 


] / 0 


% 


( 0 ) 

tors 



] ^ 0 


. . (IV. 77) 
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(b) MREV-8 sequence; 

Proceeding as in the case of WAHUHA sequence the zeroth-order 
average Hamiltonions of the internal interactions corresponding 

to efg inhomogeneities and torsional oscillations for this case 
can be shown to be 

^inhom = I .. (IV.78) 

and " "T t5o (1 +COS8 + cos^8 )] .. (iv.79) 

the is given by Eqn. IV. 76. 

Now from Eqns.IV.76/ IV.78 and IV.79 and commutation rela- 
tionshipsof tensor operators for spin I = 3/2 [see 29] it can be 
easily seen that 


[ ^ 
[ ^ 


5}(o) 

^inhom 

^ ( 0 ) 

^tors 


] / 0 
] ^ 0 


. . (IV. 80) 


Thus, the above results indicate that under WAHUHA and MREV-8, 
sequences the internal interactions defined by and 

do not average out. It will be interesting to study the internal 
interactions due to heteronuclear and homonuclear dipole-dipole 
interactions to verify whether there is any kind of selective 
averaging taking place under the influence of these sequences. 
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SUMMARY 

Experimental results of SLS, PAPS and WAHUHA multiple 

pulse sequences on ^^Cl (I = 3/2) ^®’^Re and ^®^Re (both with 

^ ~ 5/2) nuclei are presented. In SLS, the dependence of the 

decay of the quasi-equilibrium magnetization on the parameters 

such as, G, X and Au has been investigated. It is inferred 

that all the multiple pulse sequences investigated here lead to 

the lengthening of the transverse magnetization, implying a 

23 3S 

narrowed NQR line in the frequency domain. Results of Na- Cl 
double resonance experiment carried out under spin- locked condi- 
tion, in a polycrystalline sample of NaClO^ are presented. 

Using average Hamiltonian theory and tensor operator forma- 
lism theoretical analysis of the influence of PAPS, VJAHUHA, and 
and MREV-8 sequences on spin 3/2 nuclei has been carried out. In 
the case of spin 3/2 nuclei in single crystals it has been shown 
that under the action of PAPS the zeroth order average Hamilto- 
nians of the internal interactions ( ^ to 

i) efg inhomogeneities ( '^j_nhom^ 

ii) torsional oscillations ( 

iii) heteronuclear dipole-dipole interaction 

commute with P^, the initial density matrix prepared by the 
preparatory pulse, thereby implying the averaging out of these 
interactions under this sequence. In the case of polycrystalline 
samples under the action of PAPS we considered only two internal 
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interactions namely, '^tors these have been shoxvn 

to be averaged out. Under WAHUHA and MREV-8 sequences we have 

considered "Xinhom ' "^tors It has been shown that, in 

the case of single crystals, these interactions do not average 
out under these sequences . 
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CHAPITER 


V 


TU3DY OF MULTIPLE QUANTUM 
COHERENCES IN RURE NQR 



In pulsed nuclear magnetic resonance (NMR) spectroscopy 
magnetic dipole— allowed transitions between the eigen states 
of the nuclear ensemble can be observed by perturbing the spin 
system with an r.f. pulse which brings the eigen s'tates separa- 
ted by Am = +1 into coherent superposition. Such transitions 
are referred to as the single quantum transitions or single quan- 
tum coherences. Because of the perturbing r.f. pulse the inher- 
ent time-dependence of each state, normally random in thermal 
equilibrium acquires a phase coherence across the ensemble. 

This coherence is known as transverse magnetization and it gives 
rise to the free induction signal. The rigorous definition of 
coherence is based on the elements of the density matrix a [1/2] 
of the spin system. The diagonal element of the density matrix 

0 J.J. = <r| a(t)|r> = lCj^(t)l^ .. (V.l) 

is equal to the probability that the spin system is found , in the 
eigen state | r> , while the off-diagonal element 

= <Cr\ ait)\ sy = C^(t) C*{t) .. (V.2) 

indicates a "coherent superposition" of eigen states C^(t) | r> 

+ CgCt) I s> in the state function ^(t) of the ensemble, in the 
sense that the time dependence and the phase of the various 
members of the ensemble are correlated w.r .t. | r]> and ] s)> . 

Such a coherent superposition is known as "coherence" . The 
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coherence a between two states with Am = +1 leads to the 

J-a j-s _ 

observable transverse magnetization. 

In general a matrix element 0^^ represents n quantum 
coherence which for n +1, does not lead to observable magne- 
tization and can only be detected indirectly. In fact, any 
arbitrary coherence 0~ j-q cSin be made non— zero by preparing the 
spin system with appropriate r.f. pulses. This has been 
elegantly demonstrated in a number of experiments in the area 
of NMR [2] ' 

The first experimental NMR detection of double quantum 

coherence was done by Hatanka et al. [3,4j using a multistep 

27 

process. They carried out NMR experiment on Al (I = 5/2) in a 
single crystal of AI 2 O 2 which was oriented such that the efg axis 
is aligned with the external magnetic field. Out of six energy 
levels that are possible for spin I = 5/2 in high magnetic field 
(see Pig. V.l) they considered only three energy levels in their 
experiment. These levels are labelled as a, b and c in Fig. V.l. 
Hatanaka et al . have used a pulse sequence of the form 
(It /2) (Wg^) (^)j(.(<»>j 3 )-'t-(^ ' Creation and detection of 

double quantum coherence in their experiment can be understood 
as follows. 

A (n/2) pulse applied to the (a, b) transition generates 
single-quantum coherence Iiranediately if one applies a n 

pulse to (b,c) transition, it not only interchanges the popula- 
levels b and c but also swaps the indices b and c. 

Thus, the TT pulse establishes the coherence between a and c. 




CN 




> 

k. 
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which is a double quantum coherence. Since this double 
quantum coherence cannot be observed directly it has to be con- 
verted back into a single quantum coherence by a monitoring pulse 
A second n pulse applied after a time t, to b, c transition, 
converts this double quantum coherence into If the n pulse 

is applied to (a, b) .transition then cr coherence transfer into 

a.c 

^bc * either case, the observable magnetization appears in a 

transition which is not subjected to an r.f. pulse. 

Since the first detection of double quantum coherence in 
NMR by Hatanaka et al. multiple quantum NMR spectroscopy has 
been developed and applied to the solution of problems in mole- 
cular structure and dynamics [2,5-8]. However, to our knowledge 
there has been no study of this phenomenon in the area of pure 
NQR spectroscopy. 

Our interest in the study of multiple quantum coherences 
in NQR stems from the following considerations. As mentioned 
in Chapter I, the nuclei v/ith 2 possess in addition 
to quadrupole moment, electric hexadecapole moment [9]. Since 
the hexadecapole coupling constant (HDCC) is very small of the 
order of few KHz, it is rather difficult to detect it directly. 
There have been several reports in the literature [9-11] where 
the HDCC has been determined by accurate measurements of the 
NQR frequencies. However, it is well known that the NQR fre- 
quencies are rather sensitive to slight temperature changes. 

Hence anv error in the temperature measurement while measuring 
the NQR f regencies could lead to an error in the HDCC . Thus, 
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in order to measure HDCC accurately one needs to detect the two 
NQR frequencies simultaneously, it is in this context we thought 
of the possibility of the multiple quantum experiment in NQR. We 
have therefore undertaken the theoretical study of double quantum 
coherence in NQR of spin I = 5/2 nuclei in single crystals . 

This chapter deals with the theoretical investigations of 
the creation and detection of double quantum coherence in pure 
NQR of nuclei with spin I = 5/2 in a single crystal. Our approach 
is based on the density matrix theory with the tensor operator 
formlism [12-16], In the next section we present the details of, 
the theoretical investigation. 

V.A DOUBLE QUANTUM COHERENCE FOR SPIN I = 5/2 CASE IN PURE NQR 

We consider the case of a non-interacting nuclear spin 
ensemble with axially symmetric field gradient. The energy 
levels for spin I = 5/2, and the pulse scheme considered for 
the double quantum coherence are shown in Fig. V.2. In all the 
calculations in this section we assume that the relaxation effects 
are not significant. We set 1i = 1, and express energy values in 
radians per sec. The r.f. pulses are assumed to be applied along 
the principal x-axis of the efg tensor. 

The net Hamiltonian, in the laboratory frame for a nuclear 
spin ensemble with quadrupole interaction in the presence of 
perturbing r.f. field can be written as 

K(t) =1^ 


CV.3) 



(n)x A I ± 5/2 



coherence is shown on the left. 
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Where X, is the Hamiltonian an. tor the axiaiiy 

symmetric efg case it -ic ^ , r 

is given by [ 17 ] 


'H = - 2 2 

Q 41 (21-1)' [3 i:_ x^l 


(V.4) 


the Hamiltonians which define the inter- 

Bction between the nnr*!.-.-. 

ear spins and the applied r.f. pulses 

with frequencies u,, and 2 ,, 

Q nq 2 Wq, respectively. These r.f. Hamilto- 
nians can be written as 


- - 2(0 cos (CO t1 T / . 

^ Q^-' ix ..(V.5) 

and 

•= -2cou, oos(2Ugt) ..(y.6) 

where ^la ^Ib “ /^^Ib * is the gyromagnetic ratio 

of the nucleus, and^^^ and are the amplitudes of the applied 
r.f. pulses. It may be mentioned here that although in the total 
Hamiltonian we have given r.f. Hamiltonians corresponding to two 
different frequencies, for a given time only one of these Hamil- 
tonians is active (see Fig. V.l). 

Following Bowden et al. [12] the toal Hamiltonian for spin 
I = 5/2 case can be written as 


?o + cos(2o^t) Tj(a) 

.. (V.7) 

3 2 

where co^ = ^e qQ 
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AS mentioned In the previous chapters, the response of 

nuclear spin system to the r.f. pulse excitations can be 

calculated using the Liouv-ino Tr^v, m 

ouviiie-von Neumann equation. Before 

the application of the r.f. pulses the spin system is assumed 

to be at thermal equilibrium under the Influence of the pure 

quadrupolar Hamiltonian. In the high temperature approximation, 

the thermal equilibrium density matrix is given by [18] 


a(t<J0) = {Tr[exp(- ^ )])-! exp(- 


-Hc 


kT 


) 


- [Tr(ll)]“^ (1 - 


\/6 kT ° 




\/6 


.. (V.8) 


The first term on the right hand side of the above equation 
being a constant is not affected by any evolution of the spin 
system. Hence, we need to follow the evolution of the reduced 
density matrix 

P (t40) = T^ .. (V.9) 

only. 

In the quadrupolar interaction representation defined by 
the transformation operator 

U = exp [i t) .. (V.IO) 


the total Hamiltonian is given by,. 
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a = uR(t) u"^ = K 


la " - 'lb 


.. (v.ll) 


where . 2cOj^ cc,s(<o t) u iha) u'^ 


Q 

- “la [ i T?(a) - -i_ -3 

sJTs 


?!(«) 


and 




lb = 2“lb cos (2a)^t) U T^(a) 


""^Ib [f +“|^ T^^(a) + -^T^(a)] .. (V. 


U 




3^210 


. (V.12) 


9\/l5 


3\f210 


13) 


In arriving at these equations we have made use of the evolu- 
tion of the tensor operators under the influence of quadrupolar 
Hamiltonian '3^ (see Table 12 of [13]) and the high frequency 
terms have been dropped as non— secular terms . 

In order to follow the evolution of density matrix, in 

the interaction representation, following the r.f. pulses we 

need to know the evolution of tensor operators under the action 

of the Hamiltonians and . The evolution of a few 

selected tensor operators under the influence of'K^ and X, 

la lb 

can be obtained from Tables. II.5 and II. 6, respectively. 

V.A.l General Expressions for the Density Matrix at Various 
Stages of the Pulse Sequence 

From the NQR energy level diagram (see Fig.V.2) it is 
evident that double quantum cohergnce is possible between the 


states : 
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l+5/2)> « — > |+3/2> 

|-5/2> ^ |-3/2> 

|+3/2> < — ^ I - 1 / 2 ) 

|-3/2> 1+1/2) 

For spin I = 5/2 the operator set consists of 36 tensor opera- 
tors [12^. Out of these, the operators that carry information 
about the double quantum coherence (i.e.l4m| = 2) are 

2 3 4 ^ 

and T 2 ( 9- / a ) 

In what follows, we give general expressions for the density 
matrix of the spin system at various stages of the pulse 
sequence. The thermal equilibrium reduced density matrix in 
the interaction representation is given by 

p(0) = .. (V.14) 

Immediately following the application of the r.f. pulse with 
frequency 2(0^ and width ‘t^j^' the density matrix is given by 
(see Table II. 6); 

= T cos(E 3 _)] + ^ T^ (cos(E^)-l) 

+ i^ tJ(s) sin(E^) + i T^(s) sin(£^) 

where - JE W ]_]3 


(V.15) 
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Immediately after this, if we apply an r.f. pulse with a 
frequency and width 't^^,' then the density matrix just after 
the removal of the pulse is given by 




= Jo «1 + Jo ®2 +^^ = ) 83 + T^(s) 


B, 




.. (V.16) 


In writing Eqn. V.16 we have made use of the results given in 
Table II. 5. The expressions for the B's are given in Table V.l. 
Now, if we apply a (n)^ pulse with a frequency 2co^ and width 
'tm 3 ', then the density matrix immediately following the pulse 
is given by 


P + tyJ2 


+ t^,3) 


= lo =1 Jo ‘=2 = > '=3 = ) =4 

+ T^(s) C5 + T^(s) + T^(3) C, 

. . (V.17) 


Equation V.17 has been obtained in a manner similar to that of 
Eqn. V.15. The expressions for C's are given in Table V.2. 


The magnetization operators at Wq and 2 wq frequences (90° 
out of phase with the r.f. pulses applied) in the laboratory 
frame are given, respectively by (in units of rti) 

^xa = -sin(wQt) T];(a) 

= -sin(2wQt) Tj(a) 


. . (V.18) 



Table V.l. Expressions for B's 


B. 


B, 


B. 




B, 


15 
28 

A xTH' 


•f i 


Ai ig (25 + 3 oosCE^)) - ^2 
+ i sin(E 2 ) (*3 ^ tf 
A. — p= - (cos(£:„) - 1) (3 

i sinCE^) (A 3 ^. 3 . A, ^ 1 

= i sin(£3) (A3 - A2 . 

+ A3 j (2 cos (£3) + 5 cos (£2/2): 

- A^ -4^ (cos(£.) - cos(£ /2)) 
4 7 ^ 7 ^| 2 2 

1 + 5^ , ^ 


14 


= i sin(£2) (- 

- A^ ( 003 (^ 2 ) - cos (£ 2 / 2 ) 

+ A^ j (5 cos (£ 2 ) ■*■ 2 cos (£ 2/2 

= i sin "^4 

= 1 sln(E2/2) CA 3 I A^) 


where 


A 


^ = 1 (4 + 3 cos(£^)), A 2 = 


/105 


(cOS(£. 


A3 = i sin(£^) —!j-' ^4 


”• - i sin{£-) 


■1' vOTo 


(cos(£2) - 1) 


25 cos(£2)) 


) - 1 ) 


£l.= '^^C^b^wl 


£2 y8(j^^t.^2 
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Table V.2. Expressions for C's 


■f i 


Y (4 + 3 cos(£ )) + B (cos(£ ) - 1 ) 

^ 7V105 3 

i sin(£ ) B.) 

/ 3 yVTTo ^ 

2 

TToE - 1 ) + Bj i (3 + 4 oosCEj)) 

— BJ 

'-r A ' 


( 2 ^ . 4^2 

iV? 3 7 “4“' 


+ i sin(£^) (-4- 

i sintEj) (3^ ^ B^) + B 3 i (2 oos (£3 

qn . 

- cos (£ 2 / 2 )) 


"^ "4 — (cos(£ ) 

21 '/T 4 3 

i sin(£V2) Bj. - — - £_ 

3 ' 3 7V/-35 

• 2 

i sin(£_) (B. — 

3' 1 v^Xo 


^ ®4 

4 21'/T4 

B,) 


2 2 \/2 2 

+ + B - 4 = (cos 

37 q 3\/T4 


(£3/2) + 2 cos (£2 


+ i 


(£ 3 )) 

-^B,) 


= i 


+ j (5 cos (£3/ 

i sin(£o/2) (B|- a^) 

3 3\/14 I'flO ° 

= 1 sln(£3/2> (-^33+:^ P5 cosi 

= 1 sin(E3/2) (- ^ 83 4 - 

= 1 sin(E3/2) B3 + Bg) 


B.) + B, CO 

VTTT 4" 6 


£3 = '^5o^,^t^3 


5 cos (£ 3 ) ) 


COS (^ 3 / 2 ) ) 
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In the interaction representation (see Eqn . V.lo) . These are 
given by 


M 


M 




the subscripts a and b refer to the co and 2co frequencies. 

y y 

The signal in the interaction representation is propor- 
tional to the expectation value of the magnetization operator 


“C^xi^ ® '^^[0 * ^xi^ 


.. (V.20) 


V.A.2 Results 

The status of the spin system at various stages of the 
ptilse sequence can be summarized as follows. The first pulse 
with frequency 2w„ is a (^/2) ^ pulse and hence we can set 

W ^ 

t-sj 

^ 2. ~ this case. We can then obtain P using Eqn. 

( V.15). Immediately following the second pulse of vjidth 't.yj 2 ' 
and frequency to Q , we have from Eqns . V.16, V.19 and V.20 and 
using the orthogonality property of the tensor operators 



“q sJq 

7Tt~ 


sin (£ 2 ) 


and 

<«xb> 


to 


Q 


if 

2 


cos (£ 2 / 2 ) 


(V.21) 


6 kT 


. . (V.22) 
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These results indicate that, when the flip angle of the second 
pulse of frequency is different from (jt) then we get free 
induction signals at both and 2 Qq, i.e., a and b channels. 

If we set E ^ n , then from Eqn. V.21 and V.22 it is clear that 
there would not be any signal in either a or b channel. The 
density matrix Eqn. V.16 at this stage contains only the terms 
corresponding to quadrupolar and hexadecapole orders (T^ and T^) 

W.O .vo 

2 4 " 

5*2^^^' 52^®^ multipolar states which carry information about 
double quantum coherence. 


Let us now consider the situation following a (k/ 2) pulse 
of frequency 20 )^ and a second pulse of Wq with a flip angle equal 
to (^ ) . If we now apply a third (tt ) pulse of frequency 2w with- 
out any delay, i.e., = 0 (see Fig. V.2) then immediately follow- 

ing the pulse we have from Eqns. V.17, V.19 and V.20. 


<«xa> 


v/8 

6 kT ^ 


.. (V.23) 


and 



> 


0 


. . (V.24) 


Since the signal which is induced in channel 'a' following the 

2 4 

pulse in channel b is, arising entirely from T 2 (s) and T 2 (s) 
multipolar states it corresponds to the double quantum coherence 
signal. This means that the invisible double quantum coherence 
which has been created by the first two r.f. pulses in the 
sequence has been converted into an obseirvable single quantum 
coherence, in channel 'a', (detection frequency of by the 

second pulse in channel b. 
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We have also investigated the influence of resonance off- 
set on the double quantum coherence, in the particular case 
where the off-set of r.f. pulses in channel b is exactly twice 
to that of r.f. pulses in channel a. In this case using Tables 
II. 5 and II. 6 and Table 12 of [13] we have obtained for the 
signal, immediately following the last pulse in channel b (see 
Fig. V.2), the expressions 

.. (V.25) 

.. (V.26) 

where Aw = (Wq -co). 

It may be noted here that the free precession frequency of the 

double quantum coherence is equal to three times the off-set, 

v/ith respect to the r.f. carrier frequency, w , of channel 

a in this case. If the spin system is allowed to evolve under 

~ T^, during the period “t,, then we get, 

\/6 ° 

,/ jvj \ = — y [cos ('3 Ao)T- +AwT 2 )] •• (V.27) 

6 kT 

and =0 •• (V.28) 

In arriving at the Eqns. V.27 and V.28 again we have made use 

* ■ 

of the Tables II.5 and II. 6 and Table 12 of [13] . 


/ “o \/r 

s ~ "T ( 3 Aw T ) 

6 kT ^ ^ 

and 

^»xb> = 0 ' 
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From the above analysis it is evident that in principle, it 
should be possible to create and detect double quantum coherence 
in NQR of nuclear spin I = 5/2 in single crystals. 

Encouraged by the above results, we have undertaken the 

experimental investigation of double quantim coherence in NQR. 

187 ' 

We have chosen Re (I = 5/2) in KReO^ (r]-0) for this purpose. 
Since the two quadrupolar frequencies of Re are = 26.825 
MHz, V 2 - 53.626 MHz' at 296 K) very widely separated^ a pulsed NQR 
spectrometer equipped with a double resonance probe and two trans- 
mitters working at these frequencies is required to detect double 
quantum coherence. As all the pulses in the pulse scheme (see 
Fig. V.2) for the detection of double quantum coherence have to 
be applied along the x-axis of the crystal it is not possible to 
use a double resonance probe having two orthogonal coils, for 
this experiment. We have constructed a single-coil double reso- 
nance probe, based on the circuit of Kan et al. [19^ end tuned 

187 

it to the frequencies and > 2 ^ have made use of 

V4 cable (at to couple the probe to transmitter and receiver. 

A diagram of experimental set-up, including the double resonance 
probe, is shown in Fig. V.3. Further work is in progress in our 


laboratory. 



= Sample coil 



Fig. Vv3 Experimental set-up for detecting double quantum coherence. 
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SUMMARY 

Using tensor operator formalism theoretical investigation 
of double quantum coherence in NQR of spin I = 5/2 is presented. 
It has been shown that a pulse scheme of the form (n/2) (2ciu)- 
^ Cw ,) -T- (^ )^( 2(0,-,) can be used to create and detect double 

IK. w! DC W 


quantum coherence. 
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CONCLUDING REMARKS 


In this thesis tensor operetor formslisin hss been exten- 
sively used in the description of pulsed NQR experiments in 
different situations. This formalism can be’used with relative 
ease to describe the spin dynamics in the case of pure NQR also. 
However, the calculation of the time evolution of tensor opera- 
tors under the action of the Zeeman Hamiltonian becomes . rather 
cumbersome. Especially, for the case where the Zeeman field is 
oriented in a general (0, 0 ) direction with respect to principal 
axes system of electric field gradient tensor, the calculation 
of time evolution of tensor operators is not possible either 
with the technique of “nested commutation relationship" or with 
"harmonics of the motion method". The tensor operator formalism 
is ideally suited for the study of multiple quantum coherences 


in NOR. 



